postdf1d.m

１Ｄ力学模型

%*****************************************************************************
%  PostDF1D.m
%
%Post-traitement pour le probleme de differences
%finies 1D
%**************************************************************************


function []=PostDF1D(Uresult,Nx,f,E,nx,L,A,P)

%1)Affichage de la solution calculee
%=======================================================================

figure
plot(Nx,Uresult,'ko');

%2) Affichage de la solution exacte
%=========================================================================

for i=1:nx
    Uex(i)=(-f*Nx(i)^2)/(2*E*A)+((P+f*L)/(E*A))*Nx(i);
    
end

hold on
plot(Nx,Uex,'k-');

legend('u^h','u^(ex)')

xlabel('x','Fontsize',12);
ylabel('u^h(x)','Fontsize',12);


%3=Calcul de l'erreur en norme L^(inf)
%=========================================================================

Err=Uex-Uresult';
err=max(abs(Err))