u[i][j] = 1.0/(1+pow(i-JMAX/2.0,2.0)/4.0);
u[i][j] *= 1.0/(1+pow(j-JMAX/2.0,2.0)/4.0);
expect[i][j]=u[i][j];
}
}
Mgfas mgfas = new Mgfas(u,2);u=mgfas.u;
// Test that solution satisfies difference equations:";
sumsq=0.0;
maxdiff=0.0;
for (i=1;i<JMAX-1;i++) {