y[1][M]=(y[2][M]-c2[i])*y[0][M]/(2.0*(mm[i]+1.0));
scalv[0]=abs(anorm);
scalv[1]=(y[1][M] > scalv[0] ? y[1][M] : scalv[0]);
scalv[2]=(y[2][M] > 1.0 ? y[2][M] : 1.0);
Difeq difeq = new Difeq(mm[i],n[i],mpt,h,c2[i],anorm,x);
new Solvde(itmax,conv,slowc,scalv,indexv,NB,y,difeq);
result=y[2][0]+mm[i]*(mm[i]+1);
System.out.printf(" m = %d n = %d c**2 = %f lamda = %f\n",
mm[i], n[i], c2[i],result);
// System.out.printf(abs(result/expect[i]-1.0));