erroranalysis.asv

fem nonlinear source code using weighted residuals

%use this
function ErrorAnalysis

NN=[10;20;30;40;50;60;70;80;90;100]

for M=length(NN)
    N=NN(M);
   [b]=MainFileError1(N); 
   [a]=MainFileError2(N);
   error(N) = abs(b-a); 
   mesh(N) = N;
end


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

end