x[0]=x0[i];
x[1]=x1[i];
x[2]=x2[i];
x[3]=x3[i];
}
MNEWT mnewt =new MNEWT(){
public void funk(double[] x, double[] fvec, double[][] fjac) {
usrfun(x, fvec, fjac);
}
};
mnewt.mnewt(NTRIAL,x,tolx,tolf);
usrfun(x,fvec,fjac);
dy[i]=maxel(fvec);
}
System.out.printf("mnewt: Maximum discrepancy = %f\n", maxel(dy));
localflag = maxel(dy) > sbeps;