disp('Esempio di equazione diffusione-trasporto stazionaria omogenea')
disp('-a*ddu+b*du=0')
disp('numero di Peclet generale Peg=|b|/2a  e Peclet locale Peh=Peg*h')
Pe=input('Peg=')
while 0<1
   N=30
h=0.25/N
x=.75:h:1
m=length(x)
u(1)=.5; u(m)=1;
A=2*eye(m-2)-(Pe*h+1)*diag(ones(m-3,1),-1)+(Pe*h-1)*diag(ones(m-3,1),1)
F=zeros(m,1);
F(2)=.5*(Pe*h+1);
F(m-1)=-Pe*h+1;
u(2:m-1)=A\F(2:m-1);
plot(x,u,'r')
hold on

N=input('numero di passi N=')
h=0.25/N
x=.75:h:1
m=length(x)
v(1)=.5; v(m)=1;
A=2*eye(m-2)-(Pe*h+1)*diag(ones(m-3,1),-1)+(Pe*h-1)*diag(ones(m-3,1),1)
F=zeros(m,1);
F(2)=.5*(Pe*h+1);
F(m-1)=-Pe*h+1;
v(2:m-1)=A\F(2:m-1);
plot(x,v,'d',x,v)
title(['numero di Peclet Peh=',num2str(Pe*h)])
pause
hold off
clear v
end
