close all
clear
disp('Metodo di E.E. per l''equazione del calore con coeff=1')
dati_calore
%N='numero di punti spaziali interni in [0,1]=';
h=1/(N+1);
x=h:h:1-h;
A=1/h^2*(diag(-2*ones(N,1))+diag(ones(N-1,1),1)+diag(ones(N-1,1),-1));
%delta='passo temporale=';
%tf='punto finale=';
M=tf/delta;
Y=zeros(N,M+1);
%g=funzione iniziale
Y(:,1)=g(pi*x');
for i=1:M
   Y(:,i+1)=Y(:,i)+delta*A*Y(:,i);
end
Y=[zeros(1,M+1);Y;zeros(1,M+1)];
subplot(2,1,1)
contour(Y)
pause
mesh(Y)
title(['passo temporale=',num2str(delta)])
subplot(2,1,2)
clear i
z=-1+exp(i*(0:0.01:2*pi));
plot(z)
hold on
plot(real(delta*eig(A)),imag(delta*eig(A)),'o')
axis('equal')

   
   
