pde22_eq.java running
differential equation to solve
d^2u/dx^2 + d^2u/dy^2 = c(x,y)
Uxx(x,y) = 2^16 *
            12*x^2 * (1-x)^4 * y^4 * (1-y)^4 -
            32*x^3 * (1-x)^3 * y^4 * (1-y)^4 +
            12*x^4 * (1-x)^2 * y^4 * (1-y)^4;
Uyy(x,y) = 2^16 *
            12*x^4 * (1-x)^4 * y^2 * (1-y)^4 -
            32*x^4 * (1-x)^4 * y^3 * (1-y)^3 +
            12*x^4 * (1-x)^4 * y^4 * (1-y)^2;
c(x,y) = Uxx(x,y) + Uyy(x,y);
uniform grid on rectangle 0,1 to 0,1
known Solution, for testing method
u(x,y) = 2^16*x^4*(1-x)^4*y^4*(1-y)4
 
xg(0)=0.0
xg(1)=0.1111111111111111
xg(2)=0.2222222222222222
xg(3)=0.3333333333333333
xg(4)=0.4444444444444444
xg(5)=0.5555555555555556
xg(6)=0.6666666666666666
xg(7)=0.7777777777777777
xg(8)=0.8888888888888888
xg(9)=1.0
yg(0)=0.0
yg(1)=0.1111111111111111
yg(2)=0.2222222222222222
yg(3)=0.3333333333333333
yg(4)=0.4444444444444444
yg(5)=0.5555555555555556
yg(6)=0.6666666666666666
yg(7)=0.7777777777777777
yg(8)=0.8888888888888888
yg(9)=1.0
xmin=0.00000, xmax=1.00000, nx=10, hx=0.11111
ymin=0.00000, ymax=1.00000, ny=10, hy=0.11111
u(0.5,0.5)=1.00000
c(0,5,0.5)=-64.00000
internal cells zeroed 
matrix initialized
initial matrix, left side, upper 
time to compute solution = 0.07599997520446777 seconds 
check_soln against PDE
check soln against PDE max error=1.378086533776468E-11
exact solution u, computed us, error
xg[1]=0.11111, yg[1]=0.11111, u=0.00059, us=0.00059, err=-4.72E-12
xg[1]=0.11111, yg[2]=0.22222, u=0.00557, us=0.00557, err=-7.68E-12
xg[1]=0.11111, yg[3]=0.33333, u=0.01521, us=0.01521, err=-1.13E-11
xg[1]=0.11111, yg[4]=0.44444, u=0.02318, us=0.02318, err=-1.33E-11
xg[1]=0.11111, yg[5]=0.55556, u=0.02318, us=0.02318, err=-1.10E-11
xg[1]=0.11111, yg[6]=0.66667, u=0.01521, us=0.01521, err=-4.63E-12
xg[1]=0.11111, yg[7]=0.77778, u=0.00557, us=0.00557, err=2.88E-12
xg[1]=0.11111, yg[8]=0.88889, u=0.00059, us=0.00059, err=7.84E-12
xg[2]=0.22222, yg[1]=0.11111, u=0.00557, us=0.00557, err=-7.68E-12
xg[2]=0.22222, yg[2]=0.22222, u=0.05219, us=0.05219, err=-6.05E-12
xg[2]=0.22222, yg[3]=0.33333, u=0.14263, us=0.14263, err=-4.66E-12
xg[2]=0.22222, yg[4]=0.44444, u=0.21739, us=0.21739, err=-2.15E-12
xg[2]=0.22222, yg[5]=0.55556, u=0.21739, us=0.21739, err=2.28E-12
xg[2]=0.22222, yg[6]=0.66667, u=0.14263, us=0.14263, err=9.20E-12
xg[2]=0.22222, yg[7]=0.77778, u=0.05219, us=0.05219, err=1.68E-11
xg[2]=0.22222, yg[8]=0.88889, u=0.00557, us=0.00557, err=2.53E-11
xg[3]=0.33333, yg[1]=0.11111, u=0.01521, us=0.01521, err=-1.13E-11
xg[3]=0.33333, yg[2]=0.22222, u=0.14263, us=0.14263, err=-4.66E-12
xg[3]=0.33333, yg[3]=0.33333, u=0.38974, us=0.38974, err=1.32E-13
xg[3]=0.33333, yg[4]=0.44444, u=0.59403, us=0.59403, err=5.44E-12
xg[3]=0.33333, yg[5]=0.55556, u=0.59403, us=0.59403, err=1.17E-11
xg[3]=0.33333, yg[6]=0.66667, u=0.38974, us=0.38974, err=2.07E-11
xg[3]=0.33333, yg[7]=0.77778, u=0.14263, us=0.14263, err=3.05E-11
xg[3]=0.33333, yg[8]=0.88889, u=0.01521, us=0.01521, err=4.76E-11
xg[4]=0.44444, yg[1]=0.11111, u=0.02318, us=0.02318, err=-1.33E-11
xg[4]=0.44444, yg[2]=0.22222, u=0.21739, us=0.21739, err=-2.15E-12
xg[4]=0.44444, yg[3]=0.33333, u=0.59403, us=0.59403, err=5.44E-12
xg[4]=0.44444, yg[4]=0.44444, u=0.90540, us=0.90540, err=1.28E-11
xg[4]=0.44444, yg[5]=0.55556, u=0.90540, us=0.90540, err=2.02E-11
xg[4]=0.44444, yg[6]=0.66667, u=0.59403, us=0.59403, err=3.01E-11
xg[4]=0.44444, yg[7]=0.77778, u=0.21739, us=0.21739, err=4.11E-11
xg[4]=0.44444, yg[8]=0.88889, u=0.02318, us=0.02318, err=6.47E-11
xg[5]=0.55556, yg[1]=0.11111, u=0.02318, us=0.02318, err=-1.10E-11
xg[5]=0.55556, yg[2]=0.22222, u=0.21739, us=0.21739, err=2.28E-12
xg[5]=0.55556, yg[3]=0.33333, u=0.59403, us=0.59403, err=1.17E-11
xg[5]=0.55556, yg[4]=0.44444, u=0.90540, us=0.90540, err=2.02E-11
xg[5]=0.55556, yg[5]=0.55556, u=0.90540, us=0.90540, err=2.76E-11
xg[5]=0.55556, yg[6]=0.66667, u=0.59403, us=0.59403, err=3.64E-11
xg[5]=0.55556, yg[7]=0.77778, u=0.21739, us=0.21739, err=4.55E-11
xg[5]=0.55556, yg[8]=0.88889, u=0.02318, us=0.02318, err=6.69E-11
xg[6]=0.66667, yg[1]=0.11111, u=0.01521, us=0.01521, err=-4.62E-12
xg[6]=0.66667, yg[2]=0.22222, u=0.14263, us=0.14263, err=9.20E-12
xg[6]=0.66667, yg[3]=0.33333, u=0.38974, us=0.38974, err=2.07E-11
xg[6]=0.66667, yg[4]=0.44444, u=0.59403, us=0.59403, err=3.01E-11
xg[6]=0.66667, yg[5]=0.55556, u=0.59403, us=0.59403, err=3.64E-11
xg[6]=0.66667, yg[6]=0.66667, u=0.38974, us=0.38974, err=4.12E-11
xg[6]=0.66667, yg[7]=0.77778, u=0.14263, us=0.14263, err=4.43E-11
xg[6]=0.66667, yg[8]=0.88889, u=0.01521, us=0.01521, err=5.43E-11
xg[7]=0.77778, yg[1]=0.11111, u=0.00557, us=0.00557, err=2.87E-12
xg[7]=0.77778, yg[2]=0.22222, u=0.05219, us=0.05219, err=1.68E-11
xg[7]=0.77778, yg[3]=0.33333, u=0.14263, us=0.14263, err=3.05E-11
xg[7]=0.77778, yg[4]=0.44444, u=0.21739, us=0.21739, err=4.11E-11
xg[7]=0.77778, yg[5]=0.55556, u=0.21739, us=0.21739, err=4.55E-11
xg[7]=0.77778, yg[6]=0.66667, u=0.14263, us=0.14263, err=4.43E-11
xg[7]=0.77778, yg[7]=0.77778, u=0.05219, us=0.05219, err=3.96E-11
xg[7]=0.77778, yg[8]=0.88889, u=0.00557, us=0.00557, err=3.58E-11
xg[8]=0.88889, yg[1]=0.11111, u=0.00059, us=0.00059, err=7.82E-12
xg[8]=0.88889, yg[2]=0.22222, u=0.00557, us=0.00557, err=2.53E-11
xg[8]=0.88889, yg[3]=0.33333, u=0.01521, us=0.01521, err=4.76E-11
xg[8]=0.88889, yg[4]=0.44444, u=0.02318, us=0.02318, err=6.47E-11
xg[8]=0.88889, yg[5]=0.55556, u=0.02318, us=0.02318, err=6.69E-11
xg[8]=0.88889, yg[6]=0.66667, u=0.01521, us=0.01521, err=5.43E-11
xg[8]=0.88889, yg[7]=0.77778, u=0.00557, us=0.00557, err=3.58E-11
xg[8]=0.88889, yg[8]=0.88889, u=0.00059, us=0.00059, err=2.04E-11
max_error=6.692762663518259E-11, avg_error=2.3455562794163077E-11
 
total error=1.50E-09, average error=2.35E-11, max error=6.69E-11
writing pde22_eq_java.dat
finished writing pde22_eq_java.dat
total time = 0.15799999237060547 seconds 
finished pde22_eq.java