bvp_02_test: FENICS/Python version Solve -u" + u = x, 0 < x < 1 u(0) = 0, u(1) = 0 Exact solution is u(x) = x - sinh(x)/sinh(1) Using 8 elements. Solving linear variational problem. X U(X) 0.000000 0.000000 0.062500 0.009190 0.125000 0.018379 0.187500 0.026734 0.250000 0.035089 0.312500 0.041760 0.375000 0.048431 0.437500 0.052544 0.500000 0.056657 0.562500 0.057298 0.625000 0.057938 0.687500 0.054137 0.750000 0.050336 0.812500 0.041055 0.875000 0.031773 0.937500 0.015886 1.000000 0.000000 Calling FFC just-in-time (JIT) compiler, this may take some time. L2 norm of error is 0.000736338 Calling FFC just-in-time (JIT) compiler, this may take some time. H1 seminorm of error is 0.000160696 Average value of solution (my way) is 0.0373255 Calling FFC just-in-time (JIT) compiler, this may take some time. Average value of solution (FENICS way) is 0.0373255 bvp_02_test: Normal end of execution.