Optimisation d'un probleme

Salut,
je voudrais obtenir la solution de ce problème en utilisant la méthode de primal-dual interior point :
Pièce jointe 129098
Sous les contraintes :

Pi -Pwi -Phi=0
Pwi- Ppi=pvi
Ei -Ei+1+ 0.84.Ppi - 1.2 Phi=0
Ei -∆ EM ≤ 0
Phi-0.9Ei ≤ 0
0 ≤ Ppi ≤ 2
0 ≤ Pi ≤ 6
0 ≤ Phi ≤ 2
∆ EM > 0
avec 1≤i≤48 et x=(P,Pp,E,Pw,Ph,EM) total=241 variables
j'ai essayé avec fmincon et linprog mais j 'obtiens pas les mêmes courbes que celles de la source :Pièce jointe 129155
Un coup de main serait tres gentil de votre part ;)

Amitiés,