8. (contd.) Hence solve the problem (D) approximately on the meshes T2, T3, T4, T5 with K given in (4) in the case when the exact solution is u(x) = (2−x1)2 +(2−x2)2. You will need to edit the files elt stiffness.m, source term.m and femg.m. Use the quantity calculated in the last line of femg.m to estimate the rate of convergence of ǫh. (Note, however, that the identity you found in Question 3 does not hold for K 6= I in general.)

附件是前半道 这是后半道