test_bihar3d.java 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) := .24
 
in the cube   xmin<x<xmax and ymin<y<ymax and zmin<z<zmax  ,
 
   subject to boundary conditions of the first kind,
   Dirichlet, must be specified at the boundary.
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.0, xmax=2.0, ymin=0.0, ymax=2.0, zmin=0.0, zmax=2.0
nx=9, ny=9, nz=9, hx=0.25, hy=0.25, hz=0.25
maxerr PDE-f(x,y,z) is 1.1102230246251565E-16, at i_max=0, j_max=0, k_max=0
x=0.0, y=0.0, z=0.0, f(x,y,z)=0.24
check cxxxx  maxerr=3.6598057917558435E-11
check cyyyy  maxerr=2.340128091304905E-11
check cyyyy  maxerr=2.7539526215036858E-11
check cxxyy  maxerr=1.7202100854873947E-11
check cxxzz  maxerr=1.25919275006936E-11
check cyyzz  maxerr=1.8053836203790752E-11
check all terms together in this PDE
nx=9, ny=9, nz=9, hx=0.25, hy=0.25, hz=0.25
maxerr=7.639680554838435E-11, rmserr=1.360180734064457E-11, avgerr=8.699806100658133E-12
nx=9, ny=9, nz=9
hx=0.25, hy=0.25, hz=0.25
umin=-8.12, umax=1.48