%%
%Partie 2 : Q4
%On utilisera ici la fonction ode45 qui est une résolution de
%Runge-Kutta-Felhberg d'ordre 4-5
Ftemp= @(t,T) ftemp_chanfreau(t,T,tch,Ta,a,e,hc,rho,a1,a2,a3,a4,Pch)
[tode,Tode]=ode45(@(t,T)ftemp_chanfreau(t,T,tch,Ta,a,e,hc,rho,a1,a2,a3,a4,Pch) ,[0 tfin],Ta);