pde44_sparse.adb running Given uxxxx(x,y,z,t) + 2 uyyyy(x,y,z,t) + 3 uzzzz(x,y,z,t) + 4 utttt(x,y,z,t) + 5 uxyt(x,y,z,t) + 6 uyyzz(x,y,z,t) + 7 uztt(x,y,z,t) + 8 uxxx(x,y,z,t) + 9 uyyt(x,y,z,t) + 10 uzt(x,y,z,t) + 11 ut(x,y,z,t) + 12 u(x,y,z,t) = F(x,y,z,t) F(x,y,z,t) = 706 + 120 t^2 + 140 y^2z + 768 x + 180 z^2*t + 200 y^2z*t + 44 t^3 + 110 y^2z^2t + 66 xy + 12t^4 + 24 z^4 + 36 y^4 + 48 x^4 + 60 y^2z^2t^2 + 72 xyt + 84 xz + 96 t xmin<=x<=xmax ymin<=y<=ymax zmin<=z<=zmax Boundaries Analytic solution u(x,y,z,t) = t^4 + 2 z^4 + 3 y^4 + 4 x^4 + 5 y^2 z^2 t^2 + 6 x y t + 7 x z + 8 t + 9 xmin= 0.00000000000000E+00, xmax= 1.00000000000000E+00 ymin= 0.00000000000000E+00, ymax= 1.00000000000000E+00 zmin= 0.00000000000000E+00, zmax= 1.00000000000000E+00 tmin= 0.00000000000000E+00, tmax= 1.00000000000000E+00 nx= 6, ny= 6, nz= 6, nt= 6 x grid, hx= 2.00000000000000E-01 xg( 1)= 0.00000000000000E+00 xg( 2)= 2.00000000000000E-01 xg( 3)= 4.00000000000000E-01 xg( 4)= 6.00000000000000E-01 xg( 5)= 8.00000000000000E-01 xg( 6)= 1.00000000000000E+00 y grid, hy= 2.00000000000000E-01 yg( 1)= 0.00000000000000E+00 yg( 2)= 2.00000000000000E-01 yg( 3)= 4.00000000000000E-01 yg( 4)= 6.00000000000000E-01 yg( 5)= 8.00000000000000E-01 yg( 6)= 1.00000000000000E+00 z grid, hz= 2.00000000000000E-01 zg( 1)= 0.00000000000000E+00 zg( 2)= 2.00000000000000E-01 zg( 3)= 4.00000000000000E-01 zg( 4)= 6.00000000000000E-01 zg( 5)= 8.00000000000000E-01 zg( 6)= 1.00000000000000E+00 t grid, ht= 2.00000000000000E-01 tg( 1)= 0.00000000000000E+00 tg( 2)= 2.00000000000000E-01 tg( 3)= 4.00000000000000E-01 tg( 4)= 6.00000000000000E-01 tg( 5)= 8.00000000000000E-01 tg( 6)= 1.00000000000000E+00 put solution in Ua vector time for previous section = 0.000587000 seconds time for previous section = 0.000025000 seconds compute A matrix and Y RHS check rderiv(4,nx,1,hx,cxxxx cxxxx( 1)= 1.87500000000000E+03 cxxxx( 2)=-8.75000000000000E+03 cxxxx( 3)= 1.62500000000000E+04 cxxxx( 4)=-1.50000000000000E+04 cxxxx( 5)= 6.87500000000000E+03 cxxxx( 6)=-1.25000000000000E+03 check rderiv(4,nx,nx,hx,cxxxx cxxxx( 1)=-1.25000000000000E+03 cxxxx( 2)= 6.87500000000000E+03 cxxxx( 3)=-1.50000000000000E+04 cxxxx( 4)= 1.62500000000000E+04 cxxxx( 5)=-8.75000000000000E+03 cxxxx( 6)= 1.87500000000000E+03 A computed stiffness matrix time for previous section = 0.092202000 seconds solve 256 equations in 256 unknowns equations solved time for previous section = 0.285003000 seconds check pDE against known solution check_soln maxerr= 2.53547516315678E-10 check pDE against computed solution check_soln maxerr= 6.25600087289513E-07 U computed, Ua analytic, error ug( 2, 2, 2, 2)= 1.094432E+01, Ua= 1.094432E+01, err= 2.09485762070472E-11 ug( 2, 2, 2, 3)= 1.261728E+01, Ua= 1.261728E+01, err= 9.02389274415327E-12 ug( 2, 2, 2, 4)= 1.437088E+01, Ua= 1.437088E+01, err= 1.41273659437502E-11 ug( 2, 2, 2, 5)= 1.630112E+01, Ua= 1.630112E+01, err= 1.65982783073559E-11 ug( 2, 2, 3, 2)= 1.127328E+01, Ua= 1.127328E+01, err= 6.42632613789829E-11 ug( 2, 2, 3, 3)= 1.294912E+01, Ua= 1.294912E+01, err= 6.29682972430601E-11 ug( 2, 2, 3, 4)= 1.470752E+01, Ua= 1.470752E+01, err= 7.27204962913675E-11 ug( 2, 2, 3, 5)= 1.664448E+01, Ua= 1.664448E+01, err= 7.11146697085496E-11 ug( 2, 2, 4, 2)= 1.176288E+01, Ua= 1.176288E+01, err= 7.70068453448403E-11 ug( 2, 2, 4, 3)= 1.344352E+01, Ua= 1.344352E+01, err= 9.00008956250531E-11 ug( 2, 2, 4, 4)= 1.520992E+01, Ua= 1.520992E+01, err= 1.04394715094713E-10 ug( 2, 2, 4, 5)= 1.715808E+01, Ua= 1.715808E+01, err= 9.35926891543204E-11 ug( 2, 2, 5, 2)= 1.260512E+01, Ua= 1.260512E+01, err= 5.46691580893821E-11 ug( 2, 2, 5, 3)= 1.429248E+01, Ua= 1.429248E+01, err= 7.05586700178173E-11 ug( 2, 2, 5, 4)= 1.607008E+01, Ua= 1.607008E+01, err= 7.43618500109733E-11 ug( 2, 2, 5, 5)= 1.803392E+01, Ua= 1.803392E+01, err= 7.37294669761468E-11 ug( 2, 3, 2, 2)= 1.106528E+01, Ua= 1.106528E+01, err= 4.56541471294258E-11 ug( 2, 3, 2, 3)= 1.278912E+01, Ua= 1.278912E+01, err= 3.62749830173925E-11 ug( 2, 3, 2, 4)= 1.459552E+01, Ua= 1.459552E+01, err= 4.65139038396956E-11 ug( 2, 3, 2, 5)= 1.658048E+01, Ua= 1.658048E+01, err= 4.86224394080637E-11 ug( 2, 3, 3, 2)= 1.139712E+01, Ua= 1.139712E+01, err= 1.22863497153958E-10 ug( 2, 3, 3, 3)= 1.313248E+01, Ua= 1.313248E+01, err= 1.32459376800398E-10 ug( 2, 3, 3, 4)= 1.495808E+01, Ua= 1.495808E+01, err= 1.54136259311599E-10 ug( 2, 3, 3, 5)= 1.696992E+01, Ua= 1.696992E+01, err= 1.60650159841680E-10 ug( 2, 3, 4, 2)= 1.189152E+01, Ua= 1.189152E+01, err= 1.48400403077176E-10 ug( 2, 3, 4, 3)= 1.364608E+01, Ua= 1.364608E+01, err= 1.70713221336882E-10 ug( 2, 3, 4, 4)= 1.550368E+01, Ua= 1.550368E+01, err= 1.86247461897437E-10 ug( 2, 3, 4, 5)= 1.756032E+01, Ua= 1.756032E+01, err= 2.00529370886215E-10 ug( 2, 3, 5, 2)= 1.274048E+01, Ua= 1.274048E+01, err= 1.02083674846654E-10 ug( 2, 3, 5, 3)= 1.452192E+01, Ua= 1.452192E+01, err= 1.20886411991705E-10 ug( 2, 3, 5, 4)= 1.642432E+01, Ua= 1.642432E+01, err= 1.42104994438341E-10 ug( 2, 3, 5, 5)= 1.854368E+01, Ua= 1.854368E+01, err= 1.35543132273597E-10 ug( 2, 4, 2, 2)= 1.142688E+01, Ua= 1.142688E+01, err= 4.52082815627364E-11 ug( 2, 4, 2, 3)= 1.320352E+01, Ua= 1.320352E+01, err= 5.17221820928171E-11 ug( 2, 4, 2, 4)= 1.506592E+01, Ua= 1.506592E+01, err= 6.44568842744775E-11 ug( 2, 4, 2, 5)= 1.711008E+01, Ua= 1.711008E+01, err= 7.50048911868362E-11 ug( 2, 4, 3, 2)= 1.176352E+01, Ua= 1.176352E+01, err= 1.30929933561674E-10 ug( 2, 4, 3, 3)= 1.356608E+01, Ua= 1.356608E+01, err= 1.54987134237672E-10 ug( 2, 4, 3, 4)= 1.547168E+01, Ua= 1.547168E+01, err= 1.84130044544872E-10 ug( 2, 4, 3, 5)= 1.757632E+01, Ua= 1.757632E+01, err= 1.90475191175210E-10 ug( 2, 4, 4, 2)= 1.226592E+01, Ua= 1.226592E+01, err= 1.62239999212943E-10 ug( 2, 4, 4, 3)= 1.411168E+01, Ua= 1.411168E+01, err= 1.92953208966173E-10 ug( 2, 4, 4, 4)= 1.608928E+01, Ua= 1.608928E+01, err= 2.17529105839276E-10 ug( 2, 4, 4, 5)= 1.829472E+01, Ua= 1.829472E+01, err= 2.09091410852125E-10 ug( 2, 4, 5, 2)= 1.312608E+01, Ua= 1.312608E+01, err= 1.07769793089574E-10 ug( 2, 4, 5, 3)= 1.503232E+01, Ua= 1.503232E+01, err= 1.39579015012714E-10 ug( 2, 4, 5, 4)= 1.711072E+01, Ua= 1.711072E+01, err= 1.51626267097527E-10 ug( 2, 4, 5, 5)= 1.945728E+01, Ua= 1.945728E+01, err= 1.41472611403515E-10 ug( 2, 5, 2, 2)= 1.231712E+01, Ua= 1.231712E+01, err= 3.10791392621468E-11 ug( 2, 5, 2, 3)= 1.414848E+01, Ua= 1.414848E+01, err= 4.12274658856404E-11 ug( 2, 5, 2, 4)= 1.607008E+01, Ua= 1.607008E+01, err= 5.54081225345726E-11 ug( 2, 5, 2, 5)= 1.817792E+01, Ua= 1.817792E+01, err= 6.20907769643964E-11 ug( 2, 5, 3, 2)= 1.266048E+01, Ua= 1.266048E+01, err= 8.29789570389039E-11 ug( 2, 5, 3, 3)= 1.453792E+01, Ua= 1.453792E+01, err= 1.08713038571295E-10 ug( 2, 5, 3, 4)= 1.653632E+01, Ua= 1.653632E+01, err= 1.32313715539567E-10 ug( 2, 5, 3, 5)= 1.875168E+01, Ua= 1.875168E+01, err= 1.38495437340680E-10 ug( 2, 5, 4, 2)= 1.317408E+01, Ua= 1.317408E+01, err= 1.03630881653771E-10 ug( 2, 5, 4, 3)= 1.512832E+01, Ua= 1.512832E+01, err= 1.34445343746847E-10 ug( 2, 5, 4, 4)= 1.725472E+01, Ua= 1.725472E+01, err= 1.57029944602982E-10 ug( 2, 5, 4, 5)= 1.964928E+01, Ua= 1.964928E+01, err= 1.46894052477364E-10 ug( 2, 5, 5, 2)= 1.404992E+01, Ua= 1.404992E+01, err= 8.01296806685059E-11 ug( 2, 5, 5, 3)= 1.611168E+01, Ua= 1.611168E+01, err= 9.09707864593656E-11 ug( 2, 5, 5, 4)= 1.841728E+01, Ua= 1.841728E+01, err= 1.04225961194970E-10 ug( 2, 5, 5, 5)= 2.106272E+01, Ua= 2.106272E+01, err= 9.87903092664055E-11 ug( 3, 2, 2, 2)= 1.136832E+01, Ua= 1.136832E+01, err= 4.49063009000383E-12 ug( 3, 2, 2, 3)= 1.308928E+01, Ua= 1.308928E+01, err= 3.32001093283907E-12 ug( 3, 2, 2, 4)= 1.489088E+01, Ua= 1.489088E+01, err= 9.76818625986198E-12 ug( 3, 2, 2, 5)= 1.686912E+01, Ua= 1.686912E+01, err= 1.99165128833556E-11 ug( 3, 2, 3, 2)= 1.197728E+01, Ua= 1.197728E+01, err= 2.65902855289823E-11 ug( 3, 2, 3, 3)= 1.370112E+01, Ua= 1.370112E+01, err= 7.96163135419192E-12 ug( 3, 2, 3, 4)= 1.550752E+01, Ua= 1.550752E+01, err= 1.81454851144736E-11 ug( 3, 2, 3, 5)= 1.749248E+01, Ua= 1.749248E+01, err= 4.28705959620856E-11 ug( 3, 2, 4, 2)= 1.274688E+01, Ua= 1.274688E+01, err= 2.49844589461645E-11 ug( 3, 2, 4, 3)= 1.447552E+01, Ua= 1.447552E+01, err= 2.00088834390044E-11 ug( 3, 2, 4, 4)= 1.628992E+01, Ua= 1.628992E+01, err= 2.39204211993638E-11 ug( 3, 2, 4, 5)= 1.828608E+01, Ua= 1.828608E+01, err= 5.73905367673433E-11 ug( 3, 2, 5, 2)= 1.386912E+01, Ua= 1.386912E+01, err= 2.85620416207166E-11 ug( 3, 2, 5, 3)= 1.560448E+01, Ua= 1.560448E+01, err= 1.76978431909447E-11 ug( 3, 2, 5, 4)= 1.743008E+01, Ua= 1.743008E+01, err= 2.09432471365290E-11 ug( 3, 2, 5, 5)= 1.944192E+01, Ua= 1.944192E+01, err= 4.29523083766981E-11 ug( 3, 3, 2, 2)= 1.153728E+01, Ua= 1.153728E+01, err= 2.35900188272353E-11 ug( 3, 3, 2, 3)= 1.335712E+01, Ua= 1.335712E+01, err= 4.65458782628048E-11 ug( 3, 3, 2, 4)= 1.525952E+01, Ua= 1.525952E+01, err= 8.41318126276747E-11 ug( 3, 3, 2, 5)= 1.734048E+01, Ua= 1.734048E+01, err= 8.39648350847710E-11 ug( 3, 3, 3, 2)= 1.214912E+01, Ua= 1.214912E+01, err= 7.76942954416882E-11 ug( 3, 3, 3, 3)= 1.398048E+01, Ua= 1.398048E+01, err= 1.01820774034422E-10 ug( 3, 3, 3, 4)= 1.590208E+01, Ua= 1.590208E+01, err= 1.31242572365409E-10 ug( 3, 3, 3, 5)= 1.800992E+01, Ua= 1.800992E+01, err= 1.32679645048484E-10 ug( 3, 3, 4, 2)= 1.292352E+01, Ua= 1.292352E+01, err= 9.49320622112282E-11 ug( 3, 3, 4, 3)= 1.477408E+01, Ua= 1.477408E+01, err= 1.21994858659491E-10 ug( 3, 3, 4, 4)= 1.672768E+01, Ua= 1.672768E+01, err= 1.41643141660097E-10 ug( 3, 3, 4, 5)= 1.888032E+01, Ua= 1.888032E+01, err= 1.34892985670376E-10 ug( 3, 3, 5, 2)= 1.405248E+01, Ua= 1.405248E+01, err= 8.10427280839576E-11 ug( 3, 3, 5, 3)= 1.592992E+01, Ua= 1.592992E+01, err= 8.97522056675371E-11 ug( 3, 3, 5, 4)= 1.792832E+01, Ua= 1.792832E+01, err= 8.35314040159574E-11 ug( 3, 3, 5, 5)= 2.014368E+01, Ua= 2.014368E+01, err= 8.34958768791694E-11 ug( 3, 4, 2, 2)= 1.194688E+01, Ua= 1.194688E+01, err= 3.86126686180432E-11 ug( 3, 4, 2, 3)= 1.386752E+01, Ua= 1.386752E+01, err= 9.53441769979690E-11 ug( 3, 4, 2, 4)= 1.587392E+01, Ua= 1.587392E+01, err= 1.48855150428062E-10 ug( 3, 4, 2, 5)= 1.806208E+01, Ua= 1.806208E+01, err= 1.31645805367953E-10 ug( 3, 4, 3, 2)= 1.256352E+01, Ua= 1.256352E+01, err= 1.10516040763287E-10 ug( 3, 4, 3, 3)= 1.451008E+01, Ua= 1.451008E+01, err= 1.87229787229626E-10 ug( 3, 4, 3, 4)= 1.655968E+01, Ua= 1.655968E+01, err= 2.24964935569005E-10 ug( 3, 4, 3, 5)= 1.880832E+01, Ua= 1.880832E+01, err= 1.97875493768152E-10 ug( 3, 4, 4, 2)= 1.334592E+01, Ua= 1.334592E+01, err= 1.47036161024516E-10 ug( 3, 4, 4, 3)= 1.533568E+01, Ua= 1.533568E+01, err= 2.03396410825007E-10 ug( 3, 4, 4, 4)= 1.745728E+01, Ua= 1.745728E+01, err= 2.20488516333717E-10 ug( 3, 4, 4, 5)= 1.980672E+01, Ua= 1.980672E+01, err= 1.76640924109961E-10 ug( 3, 4, 5, 2)= 1.448608E+01, Ua= 1.448608E+01, err= 1.23494103831945E-10 ug( 3, 4, 5, 3)= 1.653632E+01, Ua= 1.653632E+01, err= 1.38040689989793E-10 ug( 3, 4, 5, 4)= 1.875872E+01, Ua= 1.875872E+01, err= 1.27698740470805E-10 ug( 3, 4, 5, 5)= 2.124928E+01, Ua= 2.124928E+01, err= 9.08570996216440E-11 ug( 3, 5, 2, 2)= 1.288512E+01, Ua= 1.288512E+01, err= 3.64348551329385E-11 ug( 3, 5, 2, 3)= 1.490848E+01, Ua= 1.490848E+01, err= 1.02474473351322E-10 ug( 3, 5, 2, 4)= 1.702208E+01, Ua= 1.702208E+01, err= 1.43515421768825E-10 ug( 3, 5, 2, 5)= 1.932192E+01, Ua= 1.932192E+01, err= 1.28601129745221E-10 ug( 3, 5, 3, 2)= 1.350848E+01, Ua= 1.350848E+01, err= 1.04906305864461E-10 ug( 3, 5, 3, 3)= 1.557792E+01, Ua= 1.557792E+01, err= 1.86993531769986E-10 ug( 3, 5, 3, 4)= 1.776832E+01, Ua= 1.776832E+01, err= 2.29132268714238E-10 ug( 3, 5, 3, 5)= 2.017568E+01, Ua= 2.017568E+01, err= 1.86840765081797E-10 ug( 3, 5, 4, 2)= 1.430208E+01, Ua= 1.430208E+01, err= 1.36557432028894E-10 ug( 3, 5, 4, 3)= 1.644832E+01, Ua= 1.644832E+01, err= 2.02920347192048E-10 ug( 3, 5, 4, 4)= 1.876672E+01, Ua= 1.876672E+01, err= 2.19397833234325E-10 ug( 3, 5, 4, 5)= 2.135328E+01, Ua= 2.135328E+01, err= 1.66163971471178E-10 ug( 3, 5, 5, 2)= 1.545792E+01, Ua= 1.545792E+01, err= 1.06869180171998E-10 ug( 3, 5, 5, 3)= 1.771168E+01, Ua= 1.771168E+01, err= 1.31620936372201E-10 ug( 3, 5, 5, 4)= 2.020928E+01, Ua= 2.020928E+01, err= 1.22930998713855E-10 ug( 3, 5, 5, 5)= 2.304672E+01, Ua= 2.304672E+01, err= 8.16378076251567E-11 ug( 4, 2, 2, 2)= 1.211232E+01, Ua= 1.211232E+01, err= 9.66515756317676E-12 ug( 4, 2, 2, 3)= 1.388128E+01, Ua= 1.388128E+01, err= 3.33386651618639E-11 ug( 4, 2, 2, 4)= 1.573088E+01, Ua= 1.573088E+01, err= 5.51398926518232E-11 ug( 4, 2, 2, 5)= 1.775712E+01, Ua= 1.775712E+01, err= 4.16910950207239E-11 ug( 4, 2, 3, 2)= 1.300128E+01, Ua= 1.300128E+01, err= 3.71294106571440E-11 ug( 4, 2, 3, 3)= 1.477312E+01, Ua= 1.477312E+01, err= 1.01280761555245E-10 ug( 4, 2, 3, 4)= 1.662752E+01, Ua= 1.662752E+01, err= 1.14575016141316E-10 ug( 4, 2, 3, 5)= 1.866048E+01, Ua= 1.866048E+01, err= 6.33910701708373E-11 ug( 4, 2, 4, 2)= 1.405088E+01, Ua= 1.405088E+01, err= 6.43218811546831E-11 ug( 4, 2, 4, 3)= 1.582752E+01, Ua= 1.582752E+01, err= 1.16813225758960E-10 ug( 4, 2, 4, 4)= 1.768992E+01, Ua= 1.768992E+01, err= 1.11967324301077E-10 ug( 4, 2, 4, 5)= 1.973408E+01, Ua= 1.973408E+01, err= 4.10835809816490E-11 ug( 4, 2, 5, 2)= 1.545312E+01, Ua= 1.545312E+01, err= 5.18323162168599E-11 ug( 4, 2, 5, 3)= 1.723648E+01, Ua= 1.723648E+01, err= 7.19850845598557E-11 ug( 4, 2, 5, 4)= 1.911008E+01, Ua= 1.911008E+01, err= 5.22213383646886E-11 ug( 4, 2, 5, 5)= 2.116992E+01, Ua= 2.116992E+01, err= 1.66622271535743E-12 ug( 4, 3, 2, 2)= 1.232928E+01, Ua= 1.232928E+01, err= 3.97530897089382E-11 ug( 4, 3, 2, 3)= 1.424512E+01, Ua= 1.424512E+01, err= 4.42899050767664E-11 ug( 4, 3, 2, 4)= 1.624352E+01, Ua= 1.624352E+01, err= 2.81552559044940E-11 ug( 4, 3, 2, 5)= 1.842048E+01, Ua= 1.842048E+01, err= 5.74118530494161E-12 ug( 4, 3, 3, 2)= 1.322112E+01, Ua= 1.322112E+01, err= 3.62021523869771E-12 ug( 4, 3, 3, 3)= 1.514848E+01, Ua= 1.514848E+01, err= 6.00053340349405E-12 ug( 4, 3, 3, 4)= 1.716608E+01, Ua= 1.716608E+01, err= 8.63309423948522E-13 ug( 4, 3, 3, 5)= 1.936992E+01, Ua= 1.936992E+01, err= 1.05124797755707E-11 ug( 4, 3, 4, 2)= 1.427552E+01, Ua= 1.427552E+01, err= 2.69491096105412E-11 ug( 4, 3, 4, 3)= 1.622208E+01, Ua= 1.622208E+01, err= 2.18491891246231E-11 ug( 4, 3, 4, 4)= 1.827168E+01, Ua= 1.827168E+01, err= 9.05941988094128E-13 ug( 4, 3, 4, 5)= 2.052032E+01, Ua= 2.052032E+01, err= 1.30349064875190E-11 ug( 4, 3, 5, 2)= 1.568448E+01, Ua= 1.568448E+01, err= 3.20721227353715E-11 ug( 4, 3, 5, 3)= 1.765792E+01, Ua= 1.765792E+01, err= 1.41966438604868E-11 ug( 4, 3, 5, 4)= 1.975232E+01, Ua= 1.975232E+01, err= 1.96536120711244E-11 ug( 4, 3, 5, 5)= 2.206368E+01, Ua= 2.206368E+01, err= 4.24229540385568E-11 ug( 4, 4, 2, 2)= 1.278688E+01, Ua= 1.278688E+01, err= 6.96775970254748E-11 ug( 4, 4, 2, 3)= 1.485152E+01, Ua= 1.485152E+01, err= 1.37898581442641E-10 ug( 4, 4, 2, 4)= 1.700192E+01, Ua= 1.700192E+01, err= 1.39674938282042E-10 ug( 4, 4, 2, 5)= 1.933408E+01, Ua= 1.933408E+01, err= 7.38857863780140E-11 ug( 4, 4, 3, 2)= 1.368352E+01, Ua= 1.368352E+01, err= 7.20756787586652E-11 ug( 4, 4, 3, 3)= 1.577408E+01, Ua= 1.577408E+01, err= 1.61978874757551E-10 ug( 4, 4, 3, 4)= 1.796768E+01, Ua= 1.796768E+01, err= 1.67638347647880E-10 ug( 4, 4, 3, 5)= 2.036032E+01, Ua= 2.036032E+01, err= 7.90016940754867E-11 ug( 4, 4, 4, 2)= 1.474592E+01, Ua= 1.474592E+01, err= 5.48148193502129E-11 ug( 4, 4, 4, 3)= 1.687968E+01, Ua= 1.687968E+01, err= 1.36996192168226E-10 ug( 4, 4, 4, 4)= 1.914528E+01, Ua= 1.914528E+01, err= 1.55392143597055E-10 ug( 4, 4, 4, 5)= 2.163872E+01, Ua= 2.163872E+01, err= 8.09166067483602E-11 ug( 4, 4, 5, 2)= 1.616608E+01, Ua= 1.616608E+01, err= 3.74917874523817E-11 ug( 4, 4, 5, 3)= 1.836032E+01, Ua= 1.836032E+01, err= 8.77129480159056E-11 ug( 4, 4, 5, 4)= 2.072672E+01, Ua= 2.072672E+01, err= 1.04630970554354E-10 ug( 4, 4, 5, 5)= 2.336128E+01, Ua= 2.336128E+01, err= 6.10818062796170E-11 ug( 4, 5, 2, 2)= 1.377312E+01, Ua= 1.377312E+01, err= 8.02753419293367E-11 ug( 4, 5, 2, 3)= 1.598848E+01, Ua= 1.598848E+01, err= 1.76470393853378E-10 ug( 4, 5, 2, 4)= 1.829408E+01, Ua= 1.829408E+01, err= 1.97456273554053E-10 ug( 4, 5, 2, 5)= 2.078592E+01, Ua= 2.078592E+01, err= 1.27240440406240E-10 ug( 4, 5, 3, 2)= 1.467648E+01, Ua= 1.467648E+01, err= 1.14130926931466E-10 ug( 4, 5, 3, 3)= 1.693792E+01, Ua= 1.693792E+01, err= 2.51652920724155E-10 ug( 4, 5, 3, 4)= 1.932032E+01, Ua= 1.932032E+01, err= 2.82945222807029E-10 ug( 4, 5, 3, 5)= 2.191968E+01, Ua= 2.191968E+01, err= 1.80271797489695E-10 ug( 4, 5, 4, 2)= 1.575008E+01, Ua= 1.575008E+01, err= 1.06819442180495E-10 ug( 4, 5, 4, 3)= 1.808832E+01, Ua= 1.808832E+01, err= 2.37601938124499E-10 ug( 4, 5, 4, 4)= 2.059872E+01, Ua= 2.059872E+01, err= 2.69416489118157E-10 ug( 4, 5, 4, 5)= 2.337728E+01, Ua= 2.337728E+01, err= 1.68160596558664E-10 ug( 4, 5, 5, 2)= 1.718592E+01, Ua= 1.718592E+01, err= 7.70228325563949E-11 ug( 4, 5, 5, 3)= 1.963168E+01, Ua= 1.963168E+01, err= 1.49473322608173E-10 ug( 4, 5, 5, 4)= 2.232128E+01, Ua= 2.232128E+01, err= 1.65488955872206E-10 ug( 4, 5, 5, 5)= 2.535072E+01, Ua= 2.535072E+01, err= 9.79589742655662E-11 ug( 5, 2, 2, 2)= 1.356032E+01, Ua= 1.356032E+01, err= 1.63780100592703E-12 ug( 5, 2, 2, 3)= 1.537728E+01, Ua= 1.537728E+01, err= 5.87103698990177E-11 ug( 5, 2, 2, 4)= 1.727488E+01, Ua= 1.727488E+01, err= 9.45057365697721E-11 ug( 5, 2, 2, 5)= 1.934912E+01, Ua= 1.934912E+01, err= 6.71676048114023E-11 ug( 5, 2, 3, 2)= 1.472928E+01, Ua= 1.472928E+01, err= 7.93960452938336E-11 ug( 5, 2, 3, 3)= 1.654912E+01, Ua= 1.654912E+01, err= 1.65467639590133E-10 ug( 5, 2, 3, 4)= 1.845152E+01, Ua= 1.845152E+01, err= 1.80794046400479E-10 ug( 5, 2, 3, 5)= 2.053248E+01, Ua= 2.053248E+01, err= 1.05721653653745E-10 ug( 5, 2, 4, 2)= 1.605888E+01, Ua= 1.605888E+01, err= 1.16589404797196E-10 ug( 5, 2, 4, 3)= 1.788352E+01, Ua= 1.788352E+01, err= 1.89732674016341E-10 ug( 5, 2, 4, 4)= 1.979392E+01, Ua= 1.979392E+01, err= 1.76399339579802E-10 ug( 5, 2, 4, 5)= 2.188608E+01, Ua= 2.188608E+01, err= 8.61213322878029E-11 ug( 5, 2, 5, 2)= 1.774112E+01, Ua= 1.774112E+01, err= 8.85940210082481E-11 ug( 5, 2, 5, 3)= 1.957248E+01, Ua= 1.957248E+01, err= 1.17598375481975E-10 ug( 5, 2, 5, 4)= 2.149408E+01, Ua= 2.149408E+01, err= 8.60822524373361E-11 ug( 5, 2, 5, 5)= 2.360192E+01, Ua= 2.360192E+01, err= 2.30713226301305E-11 ug( 5, 3, 2, 2)= 1.382528E+01, Ua= 1.382528E+01, err= 3.82076592586600E-11 ug( 5, 3, 2, 3)= 1.583712E+01, Ua= 1.583712E+01, err= 1.54471990754246E-11 ug( 5, 3, 2, 4)= 1.793152E+01, Ua= 1.793152E+01, err= 3.19779758228833E-11 ug( 5, 3, 2, 5)= 2.020448E+01, Ua= 2.020448E+01, err= 5.33866284513351E-11 ug( 5, 3, 3, 2)= 1.499712E+01, Ua= 1.499712E+01, err= 5.13047382355580E-11 ug( 5, 3, 3, 3)= 1.702048E+01, Ua= 1.702048E+01, err= 8.40927327772079E-11 ug( 5, 3, 3, 4)= 1.913408E+01, Ua= 1.913408E+01, err= 1.03135278095579E-10 ug( 5, 3, 3, 5)= 2.143392E+01, Ua= 2.143392E+01, err= 9.89217596725211E-11 ug( 5, 3, 4, 2)= 1.633152E+01, Ua= 1.633152E+01, err= 1.13320908212700E-10 ug( 5, 3, 4, 3)= 1.837408E+01, Ua= 1.837408E+01, err= 1.28309807223559E-10 ug( 5, 3, 4, 4)= 2.051968E+01, Ua= 2.051968E+01, err= 9.75681757608982E-11 ug( 5, 3, 4, 5)= 2.286432E+01, Ua= 2.286432E+01, err= 6.76578792990767E-11 ug( 5, 3, 5, 2)= 1.802048E+01, Ua= 1.802048E+01, err= 1.05206510170319E-10 ug( 5, 3, 5, 3)= 2.008992E+01, Ua= 2.008992E+01, err= 8.12363509794523E-11 ug( 5, 3, 5, 4)= 2.228032E+01, Ua= 2.228032E+01, err= 1.65023550380283E-11 ug( 5, 3, 5, 5)= 2.468768E+01, Ua= 2.468768E+01, err= 2.17994511331199E-11 ug( 5, 4, 2, 2)= 1.433088E+01, Ua= 1.433088E+01, err= 8.58282334093019E-11 ug( 5, 4, 2, 3)= 1.653952E+01, Ua= 1.653952E+01, err= 1.21321619417358E-10 ug( 5, 4, 2, 4)= 1.883392E+01, Ua= 1.883392E+01, err= 6.56719123526273E-11 ug( 5, 4, 2, 5)= 2.131008E+01, Ua= 2.131008E+01, err= 2.18527418383019E-11 ug( 5, 4, 3, 2)= 1.550752E+01, Ua= 1.550752E+01, err= 3.29301030888018E-11 ug( 5, 4, 3, 3)= 1.774208E+01, Ua= 1.774208E+01, err= 9.10276298782264E-11 ug( 5, 4, 3, 4)= 2.007968E+01, Ua= 2.007968E+01, err= 6.23892049134156E-11 ug( 5, 4, 3, 5)= 2.261632E+01, Ua= 2.261632E+01, err= 3.82769371753966E-11 ug( 5, 4, 4, 2)= 1.684992E+01, Ua= 1.684992E+01, err= 2.87059265247080E-11 ug( 5, 4, 4, 3)= 1.912768E+01, Ua= 1.912768E+01, err= 4.69491112653486E-11 ug( 5, 4, 4, 4)= 2.153728E+01, Ua= 2.153728E+01, err= 7.14415193669993E-11 ug( 5, 4, 4, 5)= 2.417472E+01, Ua= 2.417472E+01, err= 2.94875235340442E-13 ug( 5, 4, 5, 2)= 1.855008E+01, Ua= 1.855008E+01, err= 4.77200501336483E-11 ug( 5, 4, 5, 3)= 2.088832E+01, Ua= 2.088832E+01, err= 2.60378385519289E-11 ug( 5, 4, 5, 4)= 2.339872E+01, Ua= 2.339872E+01, err= 8.87787621195457E-11 ug( 5, 4, 5, 5)= 2.617728E+01, Ua= 2.617728E+01, err= 6.21653839516512E-11 ug( 5, 5, 2, 2)= 1.536512E+01, Ua= 1.536512E+01, err= 1.11938902591646E-10 ug( 5, 5, 2, 3)= 1.777248E+01, Ua= 1.777248E+01, err= 1.84211756959485E-10 ug( 5, 5, 2, 4)= 2.027008E+01, Ua= 2.027008E+01, err= 1.55878865371051E-10 ug( 5, 5, 2, 5)= 2.295392E+01, Ua= 2.295392E+01, err= 5.31095167843887E-11 ug( 5, 5, 3, 2)= 1.654848E+01, Ua= 1.654848E+01, err= 1.02552633052255E-10 ug( 5, 5, 3, 3)= 1.900192E+01, Ua= 1.900192E+01, err= 2.32304842029407E-10 ug( 5, 5, 3, 4)= 2.157632E+01, Ua= 2.157632E+01, err= 2.33413288697193E-10 ug( 5, 5, 3, 5)= 2.436768E+01, Ua= 2.436768E+01, err= 1.03050012967287E-10 ug( 5, 5, 4, 2)= 1.790208E+01, Ua= 1.790208E+01, err= 5.71631630919001E-11 ug( 5, 5, 4, 3)= 2.043232E+01, Ua= 2.043232E+01, err= 2.01456629156382E-10 ug( 5, 5, 4, 4)= 2.313472E+01, Ua= 2.313472E+01, err= 2.39197106566280E-10 ug( 5, 5, 4, 5)= 2.610528E+01, Ua= 2.610528E+01, err= 1.31294086713751E-10 ug( 5, 5, 5, 2)= 1.961792E+01, Ua= 1.961792E+01, err= 2.14157580558094E-11 ug( 5, 5, 5, 3)= 2.225568E+01, Ua= 2.225568E+01, err= 1.22213350550737E-10 ug( 5, 5, 5, 4)= 2.513728E+01, Ua= 2.513728E+01, err= 1.70871317095589E-10 ug( 5, 5, 5, 5)= 2.835872E+01, Ua= 2.835872E+01, err= 1.11622711074233E-10 total CPU time = 0.521446000 seconds maxerr= 2.82945222807029E-10, avgerr= 1.90537517176909E-11 end pde44_sparse.adb