/////////////////
// Equation de la chaleur 1D
// avec conditions de Dirichlet
// schémas différences finies
/////////////////////////////////////
//
deff('y=u0(x)','y=sin(%pi*x)');
N=10;
dx=1/(N+1);
dt=0.9*(1/2*dx^2);
Tfin=1;
nt=Tfin/dt;
//
uapp=u0(0:dx:1);
uapp=uapp(2:$-1)';
//
A=toeplitz([1-2*dt/(dx^2),dt/(dx^2),zeros(1,N-2)]);
//
for i=1:nt
drawlater
uapp=A*uapp;
clf()
plot2d(0:dx:1,[0;uapp;0],rect=[0,0,1,1])
drawnow
end
