fem_bihar3d_la.c running.
solve  uxxxx(x,y,z) + uyyyy(x,y,z) + uzzzz(x,y,z) + 2*uxxyy(x,y,z) + 
       2*uxxzz(x,y,z) + 2*uyyzz(x,y,z) = f(x,y,z)
 
       f(x,y,z) := 0.24  ,
 
      in the cube   xmin<x<xmax and ymin<y<ymax and zmin<z<xmax  ,
 
Dirichlet boundary conditions are computed using:
ub(x,y,z) := 0.01*x*x*x*x + 0.02*y*y*y*y + 0.03*z*z*z*z - 0.04*x*x*y*y -
             0.05*y*y*z*z - 0.06*x*x*z*z + 0.1*x*x*x*y + 
             0.3*y*y*y*z + 0.2*z*z*z*x - x*y - y*z - x*z + 1
 
xmin= 0.000, xmax= 1.000, ymin= 0.000, ymax= 1.000, zmin= 0.000, zmax= 1.000 
nx=6, ny=6, nz=6 
x grid and analytic solution at ymin zmin
i=0, xg[i]=0.000000, ub(xg[i],ymin,zmin)=1.000000
i=1, xg[i]=0.200000, ub(xg[i],ymin,zmin)=1.000016
i=2, xg[i]=0.400000, ub(xg[i],ymin,zmin)=1.000256
i=3, xg[i]=0.600000, ub(xg[i],ymin,zmin)=1.001296
i=4, xg[i]=0.800000, ub(xg[i],ymin,zmin)=1.004096
i=5, xg[i]=1.000000, ub(xg[i],ymin,zmin)=1.010000
 
y grid and analytic solution at xmin,zmin 
ii=0, yg[ii]=0.000000, ub(xmin,yg[ii],zmin)=1.000000 
ii=1, yg[ii]=0.200000, ub(xmin,yg[ii],zmin)=1.000032 
ii=2, yg[ii]=0.400000, ub(xmin,yg[ii],zmin)=1.000512 
ii=3, yg[ii]=0.600000, ub(xmin,yg[ii],zmin)=1.002592 
ii=4, yg[ii]=0.800000, ub(xmin,yg[ii],zmin)=1.008192 
ii=5, yg[ii]=1.000000, ub(xmin,yg[ii],zmin)=1.020000 
 
z grid and analytic solution at xmin,ymin 
iii=0, zg[iii]=0.000000, ub(xmin,ymin,zg[iii])=1.000000 
iii=1, zg[iii]=0.200000, ub(xmin,ymin,zg[iii])=1.000048 
iii=2, zg[iii]=0.400000, ub(xmin,ymin,zg[iii])=1.000768 
iii=3, zg[iii]=0.600000, ub(xmin,ymin,zg[iii])=1.003888 
iii=4, zg[iii]=0.800000, ub(xmin,ymin,zg[iii])=1.012288 
iii=5, zg[iii]=1.000000, ub(xmin,ymin,zg[iii])=1.030000 
 
 
calling gauleg xmin=0.000000, xmax=1.000000, npx=6 
xx[1]=0.033765, xx[2]=0.169395, wx[1]=0.085662, wx[2]=0.180381 
calling gauleg ymin=0.000000, ymax=1.000000, npy=6 
yy[1]=0.033765, yy[2]=0.169395, wy[1]=0.085662, wy[2]=0.180381 
calling gauleg zmin=0.000000, zmax=1.000000, npz=6 
zz[1]=0.033765, zz[2]=0.169395, wz[1]=0.085662, wz[2]=0.180381 
galk(xx[1],yy[1],zz[1])=74896.8
galf(xx[2],yy[2],zz[2])=0.362349
compute stiffness matrix 
k computed stiffness matrix, see above 
 f computed forcing function, see above 
 writing fem_bihar3d_la_c.dat 
finished writing fem_bihar3d_la_c.dat 
check_soln values 
maxerr=2.22738e-12, rmserr=9.87001e-13, avgerr=8.36268e-13 
ug computed Galerkin, Ua analytic, error 
ug[0,0,0]=1.000000, Ua=1.000000, err=0 
ug[0,0,1]=1.000048, Ua=1.000048, err=0 
ug[0,0,2]=1.000768, Ua=1.000768, err=0 
ug[0,0,3]=1.003888, Ua=1.003888, err=0 
ug[0,0,4]=1.012288, Ua=1.012288, err=0 
ug[0,0,5]=1.030000, Ua=1.030000, err=0 
ug[0,1,0]=1.000032, Ua=1.000032, err=0 
ug[0,1,1]=0.959952, Ua=0.959952, err=0 
ug[0,1,2]=0.920384, Ua=0.920384, err=0 
ug[0,1,3]=0.883056, Ua=0.883056, err=-1.11022e-16 
ug[0,1,4]=0.850848, Ua=0.850848, err=0 
ug[0,1,5]=0.827792, Ua=0.827792, err=0 
ug[0,2,0]=1.000512, Ua=1.000512, err=0 
ug[0,2,1]=0.919856, Ua=0.919856, err=0 
ug[0,2,2]=0.839232, Ua=0.839232, err=2.22045e-16 
ug[0,2,3]=0.760368, Ua=0.760368, err=1.11022e-16 
ug[0,2,4]=0.686144, Ua=0.686144, err=1.11022e-16 
ug[0,2,5]=0.620592, Ua=0.620592, err=0 
ug[0,3,0]=1.002592, Ua=1.002592, err=0 
ug[0,3,1]=0.880624, Ua=0.880624, err=0 
ug[0,3,2]=0.757888, Ua=0.757888, err=1.11022e-16 
ug[0,3,3]=0.636112, Ua=0.636112, err=0 
ug[0,3,4]=0.518176, Ua=0.518176, err=0 
ug[0,3,5]=0.408112, Ua=0.408112, err=0 
ug[0,4,0]=1.008192, Ua=1.008192, err=0 
ug[0,4,1]=0.843888, Ua=0.843888, err=0 
ug[0,4,2]=0.677696, Ua=0.677696, err=1.11022e-16 
ug[0,4,3]=0.511344, Ua=0.511344, err=0 
ug[0,4,4]=0.347712, Ua=0.347712, err=0 
ug[0,4,5]=0.190832, Ua=0.190832, err=0 
ug[0,5,0]=1.020000, Ua=1.020000, err=0 
ug[0,5,1]=0.812048, Ua=0.812048, err=0 
ug[0,5,2]=0.600768, Ua=0.600768, err=0 
ug[0,5,3]=0.387888, Ua=0.387888, err=0 
ug[0,5,4]=0.176288, Ua=0.176288, err=0 
ug[0,5,5]=-0.030000, Ua=-0.030000, err=0 
ug[1,0,0]=1.000016, Ua=1.000016, err=0 
ug[1,0,1]=0.960288, Ua=0.960288, err=0 
ug[1,0,2]=0.922960, Ua=0.922960, err=0 
ug[1,0,3]=0.891680, Ua=0.891680, err=0 
ug[1,0,4]=0.871248, Ua=0.871248, err=0 
ug[1,0,5]=0.867616, Ua=0.867616, err=0 
ug[1,1,0]=0.960144, Ua=0.960144, err=0 
ug[1,1,1]=-0.673362, Ua=0.880288, err=-1.55365 
ug[1,1,2]=-0.391448, Ua=0.802672, err=-1.19412 
ug[1,1,3]=-0.426198, Ua=0.730944, err=-1.15714 
ug[1,1,4]=0.008666, Ua=0.669904, err=-0.661238 
ug[1,1,5]=0.625504, Ua=0.625504, err=0 
ug[1,2,0]=0.920592, Ua=0.920592, err=0 
ug[1,2,1]=-0.398433, Ua=0.800160, err=-1.19859 
ug[1,2,2]=-1.346363, Ua=0.681488, err=-2.02785 
ug[1,2,3]=-8.223679, Ua=0.568224, err=-8.7919 
ug[1,2,4]=-0.315073, Ua=0.465168, err=-0.780241 
ug[1,2,5]=0.378272, Ua=0.378272, err=0 
ug[1,3,0]=0.882512, Ua=0.882512, err=0 
ug[1,3,1]=-0.426752, Ua=0.720768, err=-1.14752 
ug[1,3,2]=-8.065959, Ua=0.559984, err=-8.62594 
ug[1,3,3]=0.089203, Ua=0.403808, err=-0.314605 
ug[1,3,4]=0.770416, Ua=0.257040, err=0.513376 
ug[1,3,5]=0.125632, Ua=0.125632, err=0 
ug[1,4,0]=0.847824, Ua=0.847824, err=0 
ug[1,4,1]=0.035156, Ua=0.643744, err=-0.608588 
ug[1,4,2]=-0.252290, Ua=0.439504, err=-0.691794 
ug[1,4,3]=0.931148, Ua=0.238752, err=0.692396 
ug[1,4,4]=0.709835, Ua=0.046288, err=0.663547 
ug[1,4,5]=-0.131936, Ua=-0.131936, err=0 
ug[1,5,0]=0.819216, Ua=0.819216, err=0 
ug[1,5,1]=0.571488, Ua=0.571488, err=0 
ug[1,5,2]=0.322160, Ua=0.322160, err=0 
ug[1,5,3]=0.074880, Ua=0.074880, err=0 
ug[1,5,4]=-0.165552, Ua=-0.165552, err=0 
ug[1,5,5]=-0.393184, Ua=-0.393184, err=0 
ug[2,0,0]=1.000256, Ua=1.000256, err=0 
ug[2,0,1]=0.920560, Ua=0.920560, err=0 
ug[2,0,2]=0.844608, Ua=0.844608, err=0 
ug[2,0,3]=0.777968, Ua=0.777968, err=0 
ug[2,0,4]=0.727360, Ua=0.727360, err=0 
ug[2,0,5]=0.700656, Ua=0.700656, err=0 
ug[2,1,0]=0.921312, Ua=0.921312, err=0 
ug[2,1,1]=0.589043, Ua=0.801488, err=-0.212445 
ug[2,1,2]=-0.162580, Ua=0.685248, err=-0.847828 
ug[2,1,3]=-0.127666, Ua=0.578160, err=-0.705826 
ug[2,1,4]=-0.666326, Ua=0.486944, err=-1.15327 
ug[2,1,5]=0.419472, Ua=0.419472, err=0 
ug[2,2,0]=0.842304, Ua=0.842304, err=0 
ug[2,2,1]=-0.152956, Ua=0.681904, err=-0.83486 
ug[2,2,2]=0.476530, Ua=0.524608, err=-0.0480777 
ug[2,2,3]=4.463798, Ua=0.375984, err=4.08781 
ug[2,2,4]=-0.184698, Ua=0.242752, err=-0.42745 
ug[2,2,5]=0.132784, Ua=0.132784, err=0 
ug[2,3,0]=0.764384, Ua=0.764384, err=0 
ug[2,3,1]=-0.125349, Ua=0.562672, err=-0.688021 
ug[2,3,2]=4.370574, Ua=0.363264, err=4.00731 
ug[2,3,3]=-0.448664, Ua=0.171728, err=-0.620392 
ug[2,3,4]=-0.831530, Ua=-0.005216, err=-0.826314 
ug[2,3,5]=-0.159696, Ua=-0.159696, err=0 
ug[2,4,0]=0.689472, Ua=0.689472, err=0 
ug[2,4,1]=-0.698756, Ua=0.445424, err=-1.14418 
ug[2,4,2]=-0.229656, Ua=0.202560, err=-0.432216 
ug[2,4,3]=-0.932145, Ua=-0.033552, err=-0.898593 
ug[2,4,4]=-0.839404, Ua=-0.256192, err=-0.583212 
ug[2,4,5]=-0.457488, Ua=-0.457488, err=0 
ug[2,5,0]=0.620256, Ua=0.620256, err=0 
ug[2,5,1]=0.332560, Ua=0.332560, err=0 
ug[2,5,2]=0.044608, Ua=0.044608, err=0 
ug[2,5,3]=-0.238032, Ua=-0.238032, err=0 
ug[2,5,4]=-0.508640, Ua=-0.508640, err=0 
ug[2,5,5]=-0.759344, Ua=-0.759344, err=0 
ug[3,0,0]=1.001296, Ua=1.001296, err=0 
ug[3,0,1]=0.881440, Ua=0.881440, err=0 
ug[3,0,2]=0.766288, Ua=0.766288, err=0 
ug[3,0,3]=0.663328, Ua=0.663328, err=0 
ug[3,0,4]=0.581200, Ua=0.581200, err=0 
ug[3,0,5]=0.529696, Ua=0.529696, err=0 
ug[3,1,0]=0.885072, Ua=0.885072, err=0 
ug[3,1,1]=-0.254801, Ua=0.725088, err=-0.979889 
ug[3,1,2]=0.626389, Ua=0.569648, err=0.0567415 
ug[3,1,3]=0.601136, Ua=0.426240, err=0.174896 
ug[3,1,4]=0.929977, Ua=0.303504, err=0.626473 
ug[3,1,5]=0.211232, Ua=0.211232, err=0 
ug[3,2,0]=0.768144, Ua=0.768144, err=0 
ug[3,2,1]=0.620468, Ua=0.567584, err=0.0528837 
ug[3,2,2]=0.926910, Ua=0.371088, err=0.555822 
ug[3,2,3]=4.975641, Ua=0.186144, err=4.7895 
ug[3,2,4]=0.503094, Ua=0.021392, err=0.481702 
ug[3,2,5]=-0.113376, Ua=-0.113376, err=0 
ug[3,3,0]=0.651664, Ua=0.651664, err=0 
ug[3,3,1]=0.598789, Ua=0.409792, err=0.188997 
ug[3,3,2]=4.885922, Ua=0.171184, err=4.71474 
ug[3,3,3]=0.150757, Ua=-0.056672, err=0.207429 
ug[3,3,4]=-0.109234, Ua=-0.265136, err=0.155902 
ug[3,3,5]=-0.444416, Ua=-0.444416, err=0 
ug[3,4,0]=0.537552, Ua=0.537552, err=0 
ug[3,4,1]=0.948202, Ua=0.253344, err=0.694858 
ug[3,4,2]=0.483321, Ua=-0.028720, err=0.512041 
ug[3,4,3]=-0.190894, Ua=-0.301152, err=0.110258 
ug[3,4,4]=0.148729, Ua=-0.555312, err=0.704041 
ug[3,4,5]=-0.781408, Ua=-0.781408, err=0 
ug[3,5,0]=0.428496, Ua=0.428496, err=0 
ug[3,5,1]=0.100640, Ua=0.100640, err=0 
ug[3,5,2]=-0.226512, Ua=-0.226512, err=0 
ug[3,5,3]=-0.545472, Ua=-0.545472, err=0 
ug[3,5,4]=-0.847600, Ua=-0.847600, err=0 
ug[3,5,5]=-1.123104, Ua=-1.123104, err=0 
ug[4,0,0]=1.004096, Ua=1.004096, err=0 
ug[4,0,1]=0.843888, Ua=0.843888, err=0 
ug[4,0,2]=0.688960, Ua=0.688960, err=0 
ug[4,0,3]=0.548720, Ua=0.548720, err=0 
ug[4,0,4]=0.433728, Ua=0.433728, err=0 
ug[4,0,5]=0.355696, Ua=0.355696, err=0 
ug[4,1,0]=0.853344, Ua=0.853344, err=0 
ug[4,1,1]=-0.153596, Ua=0.653008, err=-0.806604 
ug[4,1,2]=-0.671219, Ua=0.457792, err=-1.12901 
ug[4,1,3]=-0.609785, Ua=0.277104, err=-0.886889 
ug[4,1,4]=-0.622370, Ua=0.121504, err=-0.743874 
ug[4,1,5]=0.002704, Ua=0.002704, err=0 
ug[4,2,0]=0.700992, Ua=0.700992, err=0 
ug[4,2,1]=-0.672692, Ua=0.460080, err=-1.13277 
ug[4,2,2]=-1.315432, Ua=0.223808, err=-1.53924 
ug[4,2,3]=-8.097614, Ua=0.001584, err=-8.0992 
ug[4,2,4]=-0.303967, Ua=-0.196032, err=-0.107935 
ug[4,2,5]=-0.357328, Ua=-0.357328, err=0 
ug[4,3,0]=0.548192, Ua=0.548192, err=0 
ug[4,3,1]=-0.609165, Ua=0.265968, err=-0.875133 
ug[4,3,2]=-7.939030, Ua=-0.012416, err=-7.92661 
ug[4,3,3]=0.365060, Ua=-0.277552, err=0.642612 
ug[4,3,4]=0.906816, Ua=-0.518880, err=1.4257 
ug[4,3,5]=-0.724688, Ua=-0.724688, err=0 
ug[4,4,0]=0.396864, Ua=0.396864, err=0 
ug[4,4,1]=-0.615903, Ua=0.072304, err=-0.688207 
ug[4,4,2]=-0.240777, Ua=-0.249536, err=0.00875854 
ug[4,4,3]=1.071128, Ua=-0.559248, err=1.63038 
ug[4,4,4]=0.670488, Ua=-0.846272, err=1.51676 
ug[4,4,5]=-1.098896, Ua=-1.098896, err=0 
ug[4,5,0]=0.249696, Ua=0.249696, err=0 
ug[4,5,1]=-0.118512, Ua=-0.118512, err=0 
ug[4,5,2]=-0.485440, Ua=-0.485440, err=0 
ug[4,5,3]=-0.841680, Ua=-0.841680, err=0 
ug[4,5,4]=-1.176672, Ua=-1.176672, err=0 
ug[4,5,5]=-1.478704, Ua=-1.478704, err=0 
ug[5,0,0]=1.010000, Ua=1.010000, err=0 
ug[5,0,1]=0.809248, Ua=0.809248, err=0 
ug[5,0,2]=0.613968, Ua=0.613968, err=0 
ug[5,0,3]=0.435488, Ua=0.435488, err=0 
ug[5,0,4]=0.286288, Ua=0.286288, err=0 
ug[5,0,5]=0.180000, Ua=0.180000, err=0 
ug[5,1,0]=0.828432, Ua=0.828432, err=0 
ug[5,1,1]=0.587552, Ua=0.587552, err=0 
ug[5,1,2]=0.351984, Ua=0.351984, err=0 
ug[5,1,3]=0.133056, Ua=0.133056, err=0 
ug[5,1,4]=-0.056752, Ua=-0.056752, err=0 
ug[5,1,5]=-0.203808, Ua=-0.203808, err=0 
ug[5,2,0]=0.644112, Ua=0.644112, err=0 
ug[5,2,1]=0.362656, Ua=0.362656, err=0 
ug[5,2,2]=0.086032, Ua=0.086032, err=0 
ug[5,2,3]=-0.174432, Ua=-0.174432, err=0 
ug[5,2,4]=-0.406256, Ua=-0.406256, err=0 
ug[5,2,5]=-0.595808, Ua=-0.595808, err=0 
ug[5,3,0]=0.458192, Ua=0.458192, err=0 
ug[5,3,1]=0.135424, Ua=0.135424, err=0 
ug[5,3,2]=-0.183312, Ua=-0.183312, err=0 
ug[5,3,3]=-0.486688, Ua=-0.486688, err=0 
ug[5,3,4]=-0.762224, Ua=-0.762224, err=0 
ug[5,3,5]=-0.996288, Ua=-0.996288, err=0 
ug[5,4,0]=0.272592, Ua=0.272592, err=0 
ug[5,4,1]=-0.092512, Ua=-0.092512, err=0 
ug[5,4,2]=-0.454704, Ua=-0.454704, err=0 
ug[5,4,3]=-0.802656, Ua=-0.802656, err=0 
ug[5,4,4]=-1.123888, Ua=-1.123888, err=0 
ug[5,4,5]=-1.404768, Ua=-1.404768, err=0 
ug[5,5,0]=0.090000, Ua=0.090000, err=0 
ug[5,5,1]=-0.318752, Ua=-0.318752, err=0 
ug[5,5,2]=-0.726032, Ua=-0.726032, err=0 
ug[5,5,3]=-1.120512, Ua=-1.120512, err=0 
ug[5,5,4]=-1.489712, Ua=-1.489712, err=0 
ug[5,5,5]=-1.820000, Ua=-1.820000, err=0 
xmax=1.000000, ymax=1.000000, zmax=1.000000, nx=6, ny=6, nz=6, npx=6, npy=6, npz=6 
                  maxerr=8.7919, avgerr=0.427343 

calling gnuplot