options=optimset('fminunc');
options = optimset(options,'Display','off');
[Theta,FVAL,EXITFLAG,OUTPUT]=fminunc('apost',[40/3.6,20,4],options, y,t,m,sigmav);
disp(['Solutions par optimisation globale de l''a posteriori :' ])
disp(['Nombre d''itérations effectuées : ' num2str(OUTPUT.iterations) ])
disp(['Vitesse estimée : ' num2str(Theta(1)*3.6) ' (Exacte :' num2str(vv) ')' ])
disp(['Ecart type estimé : ' num2str(sqrt(Theta(2))) ' (Exact :' num2str(sigma) ')'  ])
disp(['Position initiale estimée : ' num2str((Theta(3))) ' (Exact :' num2str(xo) ')'  ])