algorithme équation de la chaleur

Bonjour, s’il vous plait j’ai besoin d'une aide urgente, je cherche un algorithme ( ou un programme en langage C ) décrivant la résolution de l’équation de la chaleur d’une fonction u(x,t), x>0 , t>=0 avec comme condition au bord : u(x=0,t)= a+bcos , a et b des coefficients connus.
Merci