BEM matlab

Dear, write matlab codes for my attached 1st problem( I have 5 more), I need 4 output figures using subplot in matlab for: exact solution, approximate solution, exact-approximate, relative errors (I hope he can provide me the same in figure, i solved the same tests in other method) the code should work fast with nice convergent...........slightly like finite element(triangles, linear basis function living in these triangles, Gauss quadrature for boundary integrals............... Important remark:please follow my method not your own way to do my 1st assignment. Dear, please check carefully my 1st assigment, try to see what T_x==? mean,(normal derevative), generating mesh for square and circle domains. Calculated the integrals and the then assemble the matrix A, right hand side b: then in matlab z=A\b(remember Az=b, mean we trying to find u in Neuman part+interior nodes, since in Dirchlet part the value of u, known. Later you will try to assemble the approximate solution by putting the value of u in Dirchlet part as it is given , and Neumann+interior==z.(I'm using collocation point method) Then, as I was given above, using subplot in matlab(try to see my results in the figure given the pdf). Waiting your response to go ahead.