% L=size(p);
%
% u_a=exp(p(2,:)).*cos(p(1,:))-(exp(p(1,:)).*cos(p(2,:)));
% U=u-(u_a)';
%
% %X=(sum(U).sum(U));
% er = U.^2;
% z = sum(er);
% %iaya = sum(z)
%
% format long
% E=sqrt((1/L(:,2)*sum(u-u_a.').^2))
N=[1/177 1/665 1/2577 1/10145];
E=[0.019373732087035 0.006016619545069 0.002605405357387   0.001261971686290];

plot(N,E,'ro')
grid on
ylabel('RMS Error, e')
xlabel('1/L')
title('Graph showing dependence of RMS error with respect to 1/L using the Finite Element Method')