Equation de la chaleur, matlab

[saiyuki]

Bonjour à tous,
je dois coder sous matlab une résolution de l équation de la chaleur:
(d/dt)u(t,x)+(d²/dx²)u(t,x)=0
avec t>0 et x dans R par la technique d euler explicite.
on a: u(o,x)=a(x)
u(t,o)=b(x)
u(t,L)=c(x)
ou L représente la borne supérieur de l intervalle de discrétisation des x et T celui des t. et delta t la variation des t sur l intervalle [o,T] de meme pour delta x
On définit alors une suite u à deux indices k et n
avec o<=n<=N (N nbr de pts ds l espace) o<=k<=U(U nbr de pts ds le temps)
la methode d euler explicite nous donne la relation de recurrence pour la suite u:

u(k+1,n)=u(k,n)+(delta t/(delta x)²)(u(k,n+1)+u(k,n-1)-2u(k,n))
Voila j espere ne pas avoir été trop confus, si quelqu un pouvait m aider pour programmer cela sous matlab ca m aiderait beaucoup
Merci d avance et bonne journée