Résolution des équations différentielles stochastiques

Fermé
IMEN BELMABROUK - 19 avril 2016 à 17:41
Bonjour,
Bonjour,j'ai besoin un programe en matlab qui résoudre le systéme des équations différentiel suivant avec monte carlo simulation:
j'ai utilisé ode45 !
je cherche de dessiner la valeur moyenne et l'ecart type en fonction du temps
function dmc = MC(tc,mc)
nsample=100000;
dmc=zeros(4,1);
a=0.262;
b=0.37;
cp=a+rand(nsample,1)*(b-a);
Mont=[];
for i = 1:nsample
dmc(1)=mc(2);
dmc(3)=mc(4);
dmc(2)=-0.5*mc(1)-1/4*0.2
dmc(4)=P11*(mc(3)-mc(1))+1/0.5*cp(i)*(mc(2).^2);
end