MATLAB's default startup directory is: "/Users/jburkardt/Documents/MATLAB". Moving to preferred directory: "/Users/jburkardt/public_html/m_src". >> cd fem1d_bvp_linear >> fem1d_bvp_linear_test 18-Feb-2012 10:56:11 FEM1D_BVP_LINEAR_TEST MATLAB version Test the FEM1D_BVP_LINEAR library. FEM1D_BVP_LINEAR_TEST01 A1(X) = 1.0 C1(X) = 0.0 F1(X) = X * ( X + 3 ) * exp ( X ) U1(X) = X * ( 1 - X ) * exp ( X ) Number of nodes = 11 I X U Uexact Error 1 0.000000 0.000000 0.000000 2.220446e-16 2 0.100000 0.099466 0.099465 1.334229e-07 3 0.200000 0.195425 0.195424 2.475629e-07 4 0.300000 0.283471 0.283470 3.394330e-07 5 0.400000 0.358038 0.358038 4.056126e-07 6 0.500000 0.412181 0.412180 4.421874e-07 7 0.600000 0.437309 0.437309 4.446805e-07 8 0.700000 0.422888 0.422888 4.079761e-07 9 0.800000 0.356087 0.356087 3.262308e-07 10 0.900000 0.221364 0.221364 1.927749e-07 11 1.000000 0.000000 0.000000 0.000000e+00 L2 norm of error = 0.00400665 Seminorm of error = 0.138667 FEM1D_BVP_LINEAR_TEST02 A1(X) = 1.0 C2(X) = 2.0 F2(X) = X * ( 5 - X ) * exp ( X ) U1(X) = X * ( 1 - X ) * exp ( X ) Number of nodes = 11 I X U Uexact Error 1 0.000000 0.000000 0.000000 1.782298e-15 2 0.100000 0.099598 0.099465 1.321791e-04 3 0.200000 0.195686 0.195424 2.610606e-04 4 0.300000 0.283852 0.283470 3.818454e-04 5 0.400000 0.358526 0.358038 4.876318e-04 6 0.500000 0.412749 0.412180 5.689040e-04 7 0.600000 0.437921 0.437309 6.129042e-04 8 0.700000 0.423491 0.422888 6.028696e-04 9 0.800000 0.356604 0.356087 5.171057e-04 10 0.900000 0.221692 0.221364 3.278658e-04 11 1.000000 0.000000 0.000000 0.000000e+00 L2 norm of error = 0.00369835 Seminorm of error = 0.138675 FEM1D_BVP_LINEAR_TEST03 A1(X) = 1.0 C3(X) = 2.0 * X F3(X) = - X * ( 2 * X * X - 3 * X - 3 ) * exp ( X ) U1(X) = X * ( 1 - X ) * exp ( X ) Number of nodes = 11 I X U Uexact Error 1 0.000000 0.000000 0.000000 1.621196e-15 2 0.100000 0.099549 0.099465 8.350349e-05 3 0.200000 0.195591 0.195424 1.664831e-04 4 0.300000 0.283718 0.283470 2.473411e-04 5 0.400000 0.358361 0.358038 3.227375e-04 6 0.500000 0.412567 0.412180 3.868178e-04 7 0.600000 0.437739 0.437309 4.302058e-04 8 0.700000 0.423327 0.422888 4.386892e-04 9 0.800000 0.356478 0.356087 3.914985e-04 10 0.900000 0.221623 0.221364 2.590522e-04 11 1.000000 0.000000 0.000000 0.000000e+00 L2 norm of error = 0.00377892 Seminorm of error = 0.138671 FEM1D_BVP_LINEAR_TEST04 A2(X) = 1.0 + X * X C1(X) = 0.0 F4(X) = ( X + 3 X^2 + 5 X^3 + X^4 ) * exp ( X ) U1(X) = X * ( 1 - X ) * exp ( X ) Number of nodes = 11 I X U Uexact Error 1 0.000000 -0.000000 0.000000 6.639207e-16 2 0.100000 0.099820 0.099465 3.548374e-04 3 0.200000 0.196115 0.195424 6.903995e-04 4 0.300000 0.284455 0.283470 9.850737e-04 5 0.400000 0.359254 0.358038 1.215952e-03 6 0.500000 0.413540 0.412180 1.359969e-03 7 0.600000 0.438703 0.437309 1.394547e-03 8 0.700000 0.424186 0.422888 1.297708e-03 9 0.800000 0.357134 0.356087 1.047774e-03 10 0.900000 0.221987 0.221364 6.228182e-04 11 1.000000 0.000000 0.000000 0.000000e+00 L2 norm of error = 0.00338872 Seminorm of error = 0.138705 FEM1D_BVP_LINEAR_TEST05 A3(X) = 1.0 + X * X for X <= 1/3 = 7/9 + X for 1/3 < X C1(X) = 0.0 F5(X) = ( X + 3 X^2 + 5 X^3 + X^4 ) * exp ( X ) for X <= 1/3 = ( - 1 + 10/3 X + 43/9 X^2 + X^3 ) .* exp ( X ) for 1/3 <= X U1(X) = X * ( 1 - X ) * exp ( X ) Number of nodes = 11 I X U Uexact Error 1 0.000000 0.000000 0.000000 1.150796e-15 2 0.100000 0.099981 0.099465 5.151509e-04 3 0.200000 0.196432 0.195424 1.007893e-03 4 0.300000 0.284924 0.283470 1.453835e-03 5 0.400000 0.359566 0.358038 1.528433e-03 6 0.500000 0.413603 0.412180 1.422913e-03 7 0.600000 0.438574 0.437309 1.265587e-03 8 0.700000 0.423939 0.422888 1.051364e-03 9 0.800000 0.356861 0.356087 7.740815e-04 10 0.900000 0.221791 0.221364 4.264543e-04 11 1.000000 0.000000 0.000000 0.000000e+00 L2 norm of error = 0.00349352 Seminorm of error = 0.138709 FEM1D_BVP_LINEAR_TEST Normal end of execution. 18-Feb-2012 10:56:11 >>