pde44h_eq.f90 running solve uxxxx(x,y,z,t) + uyyyy(x,y,z,t) + uzzzz(x,y,z,t) + utttt(x,y,z,t) + 2 uxx(x,y,z,t) + 2 uyy(x,y,z,t) + 2 uzz(x,y,z,t) + 2 u(x,y,z,t) = boundary conditions computed using u(x,y,z,t) analytic solution is u(x,y,z,t) = sin(x+y+z+t) xmin= 0.00000000000000 , xmax = 1.00000000000000 ymin= 0.00000000000000 , ymax = 1.00000000000000 zmin= 0.00000000000000 , zmax = 1.00000000000000 tmin= 0.00000000000000 , tmax = 1.00000000000000 nx= 6 , ny = 6 , nz = 6 , nt= 6 x grid and analytic solution, at ymin,zmin,tmin i= 1 , Ua( 0.00000000000000 )= 0.00000000000000 i= 2 , Ua( 0.200000000000000 )= 0.198669330795061 i= 3 , Ua( 0.400000000000000 )= 0.389418342308651 i= 4 , Ua( 0.600000000000000 )= 0.564642473395035 i= 5 , Ua( 0.800000000000000 )= 0.717356090899523 i= 6 , Ua( 1.00000000000000 )= 0.841470984807897 y grid and analytic solution, at xmin,zmin,tmin ii= 1 , Ua( 0.00000000000000 )= 0.00000000000000 ii= 2 , Ua( 0.200000000000000 )= 0.198669330795061 ii= 3 , Ua( 0.400000000000000 )= 0.389418342308651 ii= 4 , Ua( 0.600000000000000 )= 0.564642473395035 ii= 5 , Ua( 0.800000000000000 )= 0.717356090899523 ii= 6 , Ua( 1.00000000000000 )= 0.841470984807897 z grid and analytic solution, at xmin,ymin,tmin iii= 1 , Ua( 0.00000000000000 )= 0.00000000000000 iii= 2 , Ua( 0.200000000000000 )= 0.198669330795061 iii= 3 , Ua( 0.400000000000000 )= 0.389418342308651 iii= 4 , Ua( 0.600000000000000 )= 0.564642473395035 iii= 5 , Ua( 0.800000000000000 )= 0.717356090899523 iii= 6 , Ua( 1.00000000000000 )= 0.841470984807897 t grid and analytic solution, at xmin,ymin,zmin iiii= 1 , Ua( 0.00000000000000 )= 0.00000000000000 iiii= 2 , Ua( 0.200000000000000 )= 0.198669330795061 iiii= 3 , Ua( 0.400000000000000 )= 0.389418342308651 iiii= 4 , Ua( 0.600000000000000 )= 0.564642473395035 iiii= 5 , Ua( 0.800000000000000 )= 0.717356090899523 iiii= 6 , Ua( 1.00000000000000 )= 0.841470984807897 compute A matrix and Y RHS check rderiv1(4,nx,1,hx,cxxxx) cxxxx( 1 )= 1875.00000000000 cxxxx( 2 )= -8750.00000000000 cxxxx( 3 )= 16250.0000000000 cxxxx( 4 )= -15000.0000000000 cxxxx( 5 )= 6875.00000000000 cxxxx( 6 )= -1250.00000000000 check rderiv1(4,nx,nx,hx,cxxxx) cxxxx( 1 )= -1250.00000000000 cxxxx( 2 )= 6875.00000000000 cxxxx( 3 )= -15000.0000000000 cxxxx( 4 )= 16250.0000000000 cxxxx( 5 )= -8750.00000000000 cxxxx( 6 )= 1875.00000000000 A and R computed check check_soln against solution check_soln maxerr= 0.128796313247735 check computed solution check_soln maxerr= 1.556468959163837E-003 write file pde44h_eq_f90.dat file written U computed Galerkin, Ua analytic, error U(2,2,2,2)= 0.7172, Ua= 0.7174, err= -0.1351092E-03 U(2,2,2,3)= 0.8412, Ua= 0.8415, err= -0.2236294E-03 U(2,2,2,4)= 0.9318, Ua= 0.9320, err= -0.2220546E-03 U(2,2,2,5)= 0.9853, Ua= 0.9854, err= -0.1318366E-03 U(2,2,3,2)= 0.8412, Ua= 0.8415, err= -0.2258835E-03 U(2,2,3,3)= 0.9317, Ua= 0.9320, err= -0.3499555E-03 U(2,2,3,4)= 0.9851, Ua= 0.9854, err= -0.3460670E-03 U(2,2,3,5)= 0.9994, Ua= 0.9996, err= -0.2183779E-03 U(2,2,4,2)= 0.9318, Ua= 0.9320, err= -0.2240139E-03 U(2,2,4,3)= 0.9851, Ua= 0.9854, err= -0.3456320E-03 U(2,2,4,4)= 0.9992, Ua= 0.9996, err= -0.3400854E-03 U(2,2,4,5)= 0.9736, Ua= 0.9738, err= -0.2134588E-03 U(2,2,5,2)= 0.9853, Ua= 0.9854, err= -0.1312666E-03 U(2,2,5,3)= 0.9994, Ua= 0.9996, err= -0.2153047E-03 U(2,2,5,4)= 0.9736, Ua= 0.9738, err= -0.2106951E-03 U(2,2,5,5)= 0.9092, Ua= 0.9093, err= -0.1230096E-03 U(2,3,2,2)= 0.8412, Ua= 0.8415, err= -0.2258836E-03 U(2,3,2,3)= 0.9317, Ua= 0.9320, err= -0.3499555E-03 U(2,3,2,4)= 0.9851, Ua= 0.9854, err= -0.3460669E-03 U(2,3,2,5)= 0.9994, Ua= 0.9996, err= -0.2183777E-03 U(2,3,3,2)= 0.9317, Ua= 0.9320, err= -0.3538019E-03 U(2,3,3,3)= 0.9849, Ua= 0.9854, err= -0.5277484E-03 U(2,3,3,4)= 0.9991, Ua= 0.9996, err= -0.5191528E-03 U(2,3,3,5)= 0.9735, Ua= 0.9738, err= -0.3368304E-03 U(2,3,4,2)= 0.9851, Ua= 0.9854, err= -0.3492281E-03 U(2,3,4,3)= 0.9991, Ua= 0.9996, err= -0.5182070E-03 U(2,3,4,4)= 0.9733, Ua= 0.9738, err= -0.5069806E-03 U(2,3,4,5)= 0.9090, Ua= 0.9093, err= -0.3270356E-03 U(2,3,5,2)= 0.9994, Ua= 0.9996, err= -0.2170721E-03 U(2,3,5,3)= 0.9735, Ua= 0.9738, err= -0.3312333E-03 U(2,3,5,4)= 0.9090, Ua= 0.9093, err= -0.3221623E-03 U(2,3,5,5)= 0.8083, Ua= 0.8085, err= -0.1996449E-03 U(2,4,2,2)= 0.9318, Ua= 0.9320, err= -0.2240142E-03 U(2,4,2,3)= 0.9851, Ua= 0.9854, err= -0.3456322E-03 U(2,4,2,4)= 0.9992, Ua= 0.9996, err= -0.3400853E-03 U(2,4,2,5)= 0.9736, Ua= 0.9738, err= -0.2134585E-03 U(2,4,3,2)= 0.9851, Ua= 0.9854, err= -0.3492283E-03 U(2,4,3,3)= 0.9991, Ua= 0.9996, err= -0.5182072E-03 U(2,4,3,4)= 0.9733, Ua= 0.9738, err= -0.5069805E-03 U(2,4,3,5)= 0.9090, Ua= 0.9093, err= -0.3270353E-03 U(2,4,4,2)= 0.9992, Ua= 0.9996, err= -0.3428021E-03 U(2,4,4,3)= 0.9733, Ua= 0.9738, err= -0.5057800E-03 U(2,4,4,4)= 0.9088, Ua= 0.9093, err= -0.4919810E-03 U(2,4,4,5)= 0.8082, Ua= 0.8085, err= -0.3154995E-03 U(2,4,5,2)= 0.9736, Ua= 0.9738, err= -0.2117953E-03 U(2,4,5,3)= 0.9090, Ua= 0.9093, err= -0.3212247E-03 U(2,4,5,4)= 0.8082, Ua= 0.8085, err= -0.3106060E-03 U(2,4,5,5)= 0.6753, Ua= 0.6755, err= -0.1915141E-03 U(2,5,2,2)= 0.9853, Ua= 0.9854, err= -0.1312668E-03 U(2,5,2,3)= 0.9994, Ua= 0.9996, err= -0.2153048E-03 U(2,5,2,4)= 0.9736, Ua= 0.9738, err= -0.2106950E-03 U(2,5,2,5)= 0.9092, Ua= 0.9093, err= -0.1230094E-03 U(2,5,3,2)= 0.9994, Ua= 0.9996, err= -0.2170724E-03 U(2,5,3,3)= 0.9735, Ua= 0.9738, err= -0.3312335E-03 U(2,5,3,4)= 0.9090, Ua= 0.9093, err= -0.3221621E-03 U(2,5,3,5)= 0.8083, Ua= 0.8085, err= -0.1996446E-03 U(2,5,4,2)= 0.9736, Ua= 0.9738, err= -0.2117955E-03 U(2,5,4,3)= 0.9090, Ua= 0.9093, err= -0.3212248E-03 U(2,5,4,4)= 0.8082, Ua= 0.8085, err= -0.3106060E-03 U(2,5,4,5)= 0.6753, Ua= 0.6755, err= -0.1915139E-03 U(2,5,5,2)= 0.9092, Ua= 0.9093, err= -0.1217466E-03 U(2,5,5,3)= 0.8083, Ua= 0.8085, err= -0.1961340E-03 U(2,5,5,4)= 0.6753, Ua= 0.6755, err= -0.1886838E-03 U(2,5,5,5)= 0.5154, Ua= 0.5155, err= -0.1087933E-03 U(3,2,2,2)= 0.8412, Ua= 0.8415, err= -0.2258834E-03 U(3,2,2,3)= 0.9317, Ua= 0.9320, err= -0.3499555E-03 U(3,2,2,4)= 0.9851, Ua= 0.9854, err= -0.3460670E-03 U(3,2,2,5)= 0.9994, Ua= 0.9996, err= -0.2183779E-03 U(3,2,3,2)= 0.9317, Ua= 0.9320, err= -0.3538016E-03 U(3,2,3,3)= 0.9849, Ua= 0.9854, err= -0.5277483E-03 U(3,2,3,4)= 0.9991, Ua= 0.9996, err= -0.5191529E-03 U(3,2,3,5)= 0.9735, Ua= 0.9738, err= -0.3368307E-03 U(3,2,4,2)= 0.9851, Ua= 0.9854, err= -0.3492281E-03 U(3,2,4,3)= 0.9991, Ua= 0.9996, err= -0.5182071E-03 U(3,2,4,4)= 0.9733, Ua= 0.9738, err= -0.5069806E-03 U(3,2,4,5)= 0.9090, Ua= 0.9093, err= -0.3270357E-03 U(3,2,5,2)= 0.9994, Ua= 0.9996, err= -0.2170721E-03 U(3,2,5,3)= 0.9735, Ua= 0.9738, err= -0.3312334E-03 U(3,2,5,4)= 0.9090, Ua= 0.9093, err= -0.3221622E-03 U(3,2,5,5)= 0.8083, Ua= 0.8085, err= -0.1996448E-03 U(3,3,2,2)= 0.9317, Ua= 0.9320, err= -0.3538019E-03 U(3,3,2,3)= 0.9849, Ua= 0.9854, err= -0.5277482E-03 U(3,3,2,4)= 0.9991, Ua= 0.9996, err= -0.5191528E-03 U(3,3,2,5)= 0.9735, Ua= 0.9738, err= -0.3368306E-03 U(3,3,3,2)= 0.9849, Ua= 0.9854, err= -0.5337199E-03 U(3,3,3,3)= 0.9988, Ua= 0.9996, err= -0.7780414E-03 U(3,3,3,4)= 0.9731, Ua= 0.9738, err= -0.7610383E-03 U(3,3,3,5)= 0.9088, Ua= 0.9093, err= -0.4994000E-03 U(3,3,4,2)= 0.9990, Ua= 0.9996, err= -0.5237721E-03 U(3,3,4,3)= 0.9731, Ua= 0.9738, err= -0.7592327E-03 U(3,3,4,4)= 0.9086, Ua= 0.9093, err= -0.7383141E-03 U(3,3,4,5)= 0.8080, Ua= 0.8085, err= -0.4814375E-03 U(3,3,5,2)= 0.9735, Ua= 0.9738, err= -0.3342832E-03 U(3,3,5,3)= 0.9088, Ua= 0.9093, err= -0.4901216E-03 U(3,3,5,4)= 0.8080, Ua= 0.8085, err= -0.4735799E-03 U(3,3,5,5)= 0.6752, Ua= 0.6755, err= -0.3012212E-03 U(3,4,2,2)= 0.9851, Ua= 0.9854, err= -0.3492283E-03 U(3,4,2,3)= 0.9991, Ua= 0.9996, err= -0.5182071E-03 U(3,4,2,4)= 0.9733, Ua= 0.9738, err= -0.5069805E-03 U(3,4,2,5)= 0.9090, Ua= 0.9093, err= -0.3270353E-03 U(3,4,3,2)= 0.9990, Ua= 0.9996, err= -0.5237722E-03 U(3,4,3,3)= 0.9731, Ua= 0.9738, err= -0.7592328E-03 U(3,4,3,4)= 0.9086, Ua= 0.9093, err= -0.7383140E-03 U(3,4,3,5)= 0.8080, Ua= 0.8085, err= -0.4814372E-03 U(3,4,4,2)= 0.9733, Ua= 0.9738, err= -0.5109078E-03 U(3,4,4,3)= 0.9086, Ua= 0.9093, err= -0.7361416E-03 U(3,4,4,4)= 0.8078, Ua= 0.8085, err= -0.7114330E-03 U(3,4,4,5)= 0.6750, Ua= 0.6755, err= -0.4608183E-03 U(3,4,5,2)= 0.9090, Ua= 0.9093, err= -0.3239676E-03 U(3,4,5,3)= 0.8080, Ua= 0.8085, err= -0.4719044E-03 U(3,4,5,4)= 0.6750, Ua= 0.6755, err= -0.4529894E-03 U(3,4,5,5)= 0.5152, Ua= 0.5155, err= -0.2862094E-03 U(3,5,2,2)= 0.9994, Ua= 0.9996, err= -0.2170725E-03 U(3,5,2,3)= 0.9735, Ua= 0.9738, err= -0.3312335E-03 U(3,5,2,4)= 0.9090, Ua= 0.9093, err= -0.3221621E-03 U(3,5,2,5)= 0.8083, Ua= 0.8085, err= -0.1996446E-03 U(3,5,3,2)= 0.9735, Ua= 0.9738, err= -0.3342835E-03 U(3,5,3,3)= 0.9088, Ua= 0.9093, err= -0.4901217E-03 U(3,5,3,4)= 0.8080, Ua= 0.8085, err= -0.4735797E-03 U(3,5,3,5)= 0.6752, Ua= 0.6755, err= -0.3012210E-03 U(3,5,4,2)= 0.9090, Ua= 0.9093, err= -0.3239678E-03 U(3,5,4,3)= 0.8080, Ua= 0.8085, err= -0.4719044E-03 U(3,5,4,4)= 0.6750, Ua= 0.6755, err= -0.4529893E-03 U(3,5,4,5)= 0.5152, Ua= 0.5155, err= -0.2862093E-03 U(3,5,5,2)= 0.8083, Ua= 0.8085, err= -0.1973165E-03 U(3,5,5,3)= 0.6752, Ua= 0.6755, err= -0.2951012E-03 U(3,5,5,4)= 0.5152, Ua= 0.5155, err= -0.2813683E-03 U(3,5,5,5)= 0.3348, Ua= 0.3350, err= -0.1709907E-03 U(4,2,2,2)= 0.9318, Ua= 0.9320, err= -0.2240138E-03 U(4,2,2,3)= 0.9851, Ua= 0.9854, err= -0.3456320E-03 U(4,2,2,4)= 0.9992, Ua= 0.9996, err= -0.3400855E-03 U(4,2,2,5)= 0.9736, Ua= 0.9738, err= -0.2134590E-03 U(4,2,3,2)= 0.9851, Ua= 0.9854, err= -0.3492279E-03 U(4,2,3,3)= 0.9991, Ua= 0.9996, err= -0.5182070E-03 U(4,2,3,4)= 0.9733, Ua= 0.9738, err= -0.5069808E-03 U(4,2,3,5)= 0.9090, Ua= 0.9093, err= -0.3270358E-03 U(4,2,4,2)= 0.9992, Ua= 0.9996, err= -0.3428020E-03 U(4,2,4,3)= 0.9733, Ua= 0.9738, err= -0.5057800E-03 U(4,2,4,4)= 0.9088, Ua= 0.9093, err= -0.4919811E-03 U(4,2,4,5)= 0.8082, Ua= 0.8085, err= -0.3154996E-03 U(4,2,5,2)= 0.9736, Ua= 0.9738, err= -0.2117955E-03 U(4,2,5,3)= 0.9090, Ua= 0.9093, err= -0.3212248E-03 U(4,2,5,4)= 0.8082, Ua= 0.8085, err= -0.3106060E-03 U(4,2,5,5)= 0.6753, Ua= 0.6755, err= -0.1915139E-03 U(4,3,2,2)= 0.9851, Ua= 0.9854, err= -0.3492279E-03 U(4,3,2,3)= 0.9991, Ua= 0.9996, err= -0.5182070E-03 U(4,3,2,4)= 0.9733, Ua= 0.9738, err= -0.5069807E-03 U(4,3,2,5)= 0.9090, Ua= 0.9093, err= -0.3270357E-03 U(4,3,3,2)= 0.9990, Ua= 0.9996, err= -0.5237719E-03 U(4,3,3,3)= 0.9731, Ua= 0.9738, err= -0.7592327E-03 U(4,3,3,4)= 0.9086, Ua= 0.9093, err= -0.7383142E-03 U(4,3,3,5)= 0.8080, Ua= 0.8085, err= -0.4814376E-03 U(4,3,4,2)= 0.9733, Ua= 0.9738, err= -0.5109077E-03 U(4,3,4,3)= 0.9086, Ua= 0.9093, err= -0.7361416E-03 U(4,3,4,4)= 0.8078, Ua= 0.8085, err= -0.7114331E-03 U(4,3,4,5)= 0.6750, Ua= 0.6755, err= -0.4608184E-03 U(4,3,5,2)= 0.9090, Ua= 0.9093, err= -0.3239677E-03 U(4,3,5,3)= 0.8080, Ua= 0.8085, err= -0.4719044E-03 U(4,3,5,4)= 0.6750, Ua= 0.6755, err= -0.4529894E-03 U(4,3,5,5)= 0.5152, Ua= 0.5155, err= -0.2862094E-03 U(4,4,2,2)= 0.9992, Ua= 0.9996, err= -0.3428021E-03 U(4,4,2,3)= 0.9733, Ua= 0.9738, err= -0.5057800E-03 U(4,4,2,4)= 0.9088, Ua= 0.9093, err= -0.4919810E-03 U(4,4,2,5)= 0.8082, Ua= 0.8085, err= -0.3154995E-03 U(4,4,3,2)= 0.9733, Ua= 0.9738, err= -0.5109078E-03 U(4,4,3,3)= 0.9086, Ua= 0.9093, err= -0.7361416E-03 U(4,4,3,4)= 0.8078, Ua= 0.8085, err= -0.7114329E-03 U(4,4,3,5)= 0.6750, Ua= 0.6755, err= -0.4608183E-03 U(4,4,4,2)= 0.9088, Ua= 0.9093, err= -0.4951949E-03 U(4,4,4,3)= 0.8078, Ua= 0.8085, err= -0.7089062E-03 U(4,4,4,4)= 0.6748, Ua= 0.6755, err= -0.6805608E-03 U(4,4,4,5)= 0.5151, Ua= 0.5155, err= -0.4376693E-03 U(4,4,5,2)= 0.8082, Ua= 0.8085, err= -0.3119289E-03 U(4,4,5,3)= 0.6750, Ua= 0.6755, err= -0.4510844E-03 U(4,4,5,4)= 0.5151, Ua= 0.5155, err= -0.4299126E-03 U(4,4,5,5)= 0.3347, Ua= 0.3350, err= -0.2697047E-03 U(4,5,2,2)= 0.9736, Ua= 0.9738, err= -0.2117957E-03 U(4,5,2,3)= 0.9090, Ua= 0.9093, err= -0.3212249E-03 U(4,5,2,4)= 0.8082, Ua= 0.8085, err= -0.3106059E-03 U(4,5,2,5)= 0.6753, Ua= 0.6755, err= -0.1915137E-03 U(4,5,3,2)= 0.9090, Ua= 0.9093, err= -0.3239679E-03 U(4,5,3,3)= 0.8080, Ua= 0.8085, err= -0.4719044E-03 U(4,5,3,4)= 0.6750, Ua= 0.6755, err= -0.4529892E-03 U(4,5,3,5)= 0.5152, Ua= 0.5155, err= -0.2862092E-03 U(4,5,4,2)= 0.8082, Ua= 0.8085, err= -0.3119290E-03 U(4,5,4,3)= 0.6750, Ua= 0.6755, err= -0.4510845E-03 U(4,5,4,4)= 0.5151, Ua= 0.5155, err= -0.4299126E-03 U(4,5,4,5)= 0.3347, Ua= 0.3350, err= -0.2697045E-03 U(4,5,5,2)= 0.6753, Ua= 0.6755, err= -0.1888630E-03 U(4,5,5,3)= 0.5152, Ua= 0.5155, err= -0.2799791E-03 U(4,5,5,4)= 0.3347, Ua= 0.3350, err= -0.2649269E-03 U(4,5,5,5)= 0.1410, Ua= 0.1411, err= -0.1601051E-03 U(5,2,2,2)= 0.9853, Ua= 0.9854, err= -0.1312663E-03 U(5,2,2,3)= 0.9994, Ua= 0.9996, err= -0.2153046E-03 U(5,2,2,4)= 0.9736, Ua= 0.9738, err= -0.2106952E-03 U(5,2,2,5)= 0.9092, Ua= 0.9093, err= -0.1230098E-03 U(5,2,3,2)= 0.9994, Ua= 0.9996, err= -0.2170719E-03 U(5,2,3,3)= 0.9735, Ua= 0.9738, err= -0.3312333E-03 U(5,2,3,4)= 0.9090, Ua= 0.9093, err= -0.3221624E-03 U(5,2,3,5)= 0.8083, Ua= 0.8085, err= -0.1996450E-03 U(5,2,4,2)= 0.9736, Ua= 0.9738, err= -0.2117954E-03 U(5,2,4,3)= 0.9090, Ua= 0.9093, err= -0.3212248E-03 U(5,2,4,4)= 0.8082, Ua= 0.8085, err= -0.3106061E-03 U(5,2,4,5)= 0.6753, Ua= 0.6755, err= -0.1915141E-03 U(5,2,5,2)= 0.9092, Ua= 0.9093, err= -0.1217467E-03 U(5,2,5,3)= 0.8083, Ua= 0.8085, err= -0.1961341E-03 U(5,2,5,4)= 0.6753, Ua= 0.6755, err= -0.1886838E-03 U(5,2,5,5)= 0.5154, Ua= 0.5155, err= -0.1087932E-03 U(5,3,2,2)= 0.9994, Ua= 0.9996, err= -0.2170720E-03 U(5,3,2,3)= 0.9735, Ua= 0.9738, err= -0.3312333E-03 U(5,3,2,4)= 0.9090, Ua= 0.9093, err= -0.3221624E-03 U(5,3,2,5)= 0.8083, Ua= 0.8085, err= -0.1996450E-03 U(5,3,3,2)= 0.9735, Ua= 0.9738, err= -0.3342831E-03 U(5,3,3,3)= 0.9088, Ua= 0.9093, err= -0.4901216E-03 U(5,3,3,4)= 0.8080, Ua= 0.8085, err= -0.4735800E-03 U(5,3,3,5)= 0.6752, Ua= 0.6755, err= -0.3012215E-03 U(5,3,4,2)= 0.9090, Ua= 0.9093, err= -0.3239677E-03 U(5,3,4,3)= 0.8080, Ua= 0.8085, err= -0.4719044E-03 U(5,3,4,4)= 0.6750, Ua= 0.6755, err= -0.4529894E-03 U(5,3,4,5)= 0.5152, Ua= 0.5155, err= -0.2862095E-03 U(5,3,5,2)= 0.8083, Ua= 0.8085, err= -0.1973166E-03 U(5,3,5,3)= 0.6752, Ua= 0.6755, err= -0.2951014E-03 U(5,3,5,4)= 0.5152, Ua= 0.5155, err= -0.2813683E-03 U(5,3,5,5)= 0.3348, Ua= 0.3350, err= -0.1709906E-03 U(5,4,2,2)= 0.9736, Ua= 0.9738, err= -0.2117954E-03 U(5,4,2,3)= 0.9090, Ua= 0.9093, err= -0.3212248E-03 U(5,4,2,4)= 0.8082, Ua= 0.8085, err= -0.3106060E-03 U(5,4,2,5)= 0.6753, Ua= 0.6755, err= -0.1915140E-03 U(5,4,3,2)= 0.9090, Ua= 0.9093, err= -0.3239677E-03 U(5,4,3,3)= 0.8080, Ua= 0.8085, err= -0.4719043E-03 U(5,4,3,4)= 0.6750, Ua= 0.6755, err= -0.4529893E-03 U(5,4,3,5)= 0.5152, Ua= 0.5155, err= -0.2862094E-03 U(5,4,4,2)= 0.8082, Ua= 0.8085, err= -0.3119290E-03 U(5,4,4,3)= 0.6750, Ua= 0.6755, err= -0.4510845E-03 U(5,4,4,4)= 0.5151, Ua= 0.5155, err= -0.4299126E-03 U(5,4,4,5)= 0.3347, Ua= 0.3350, err= -0.2697046E-03 U(5,4,5,2)= 0.6753, Ua= 0.6755, err= -0.1888632E-03 U(5,4,5,3)= 0.5152, Ua= 0.5155, err= -0.2799791E-03 U(5,4,5,4)= 0.3347, Ua= 0.3350, err= -0.2649269E-03 U(5,4,5,5)= 0.1410, Ua= 0.1411, err= -0.1601050E-03 U(5,5,2,2)= 0.9092, Ua= 0.9093, err= -0.1217467E-03 U(5,5,2,3)= 0.8083, Ua= 0.8085, err= -0.1961341E-03 U(5,5,2,4)= 0.6753, Ua= 0.6755, err= -0.1886837E-03 U(5,5,2,5)= 0.5154, Ua= 0.5155, err= -0.1087931E-03 U(5,5,3,2)= 0.8083, Ua= 0.8085, err= -0.1973167E-03 U(5,5,3,3)= 0.6752, Ua= 0.6755, err= -0.2951013E-03 U(5,5,3,4)= 0.5152, Ua= 0.5155, err= -0.2813682E-03 U(5,5,3,5)= 0.3348, Ua= 0.3350, err= -0.1709905E-03 U(5,5,4,2)= 0.6753, Ua= 0.6755, err= -0.1888632E-03 U(5,5,4,3)= 0.5152, Ua= 0.5155, err= -0.2799791E-03 U(5,5,4,4)= 0.3347, Ua= 0.3350, err= -0.2649269E-03 U(5,5,4,5)= 0.1410, Ua= 0.1411, err= -0.1601050E-03 U(5,5,5,2)= 0.5154, Ua= 0.5155, err= -0.1069057E-03 U(5,5,5,3)= 0.3348, Ua= 0.3350, err= -0.1672339E-03 U(5,5,5,4)= 0.1410, Ua= 0.1411, err= -0.1573623E-03 U(5,5,5,5)=-0.0585, Ua=-0.0584, err= -0.8991495E-04 maxerr= 7.780413724703950E-004 , avgerr= 3.457225480969355E-004 end pde44h_eq.f90