102103104105106107108109110111112
} // Test method fredin() for (i=0;i<M;i++) { z=(b-a)*i/(M-1); ff[i]=fred.fredin(z); ffexp[i]=sqrt(z); } sbeps = 2.0e-6; System.out.println(maxel(vecsub(ff,ffexp))); localflag = maxel(vecsub(ff,ffexp)) > sbeps;