erroranalysis.m

fem nonlinear source code using weighted residuals

%use this
function ErrorAnalysis
clc
NN=[10;20;30;40;50;60;70;80;90;100];
lengthNN = length(NN);

for M=1:lengthNN
    N=NN(M);
    [b]=MainFileError1(N);
    [a]=MainFileError2(N);
%     disp(b)
% disp(a)
    error(M) = norm(b-a);
end


disp(error');
% disp(mesh')
plot(NN,error);

end