Hi,

In fact, I would like to plot C in function of Vlow, the problem is that I am not getting the results displayed on the plot, Here is the code


Vlow=350;
fs=2*100000;
n=0.56;
Lk=41.98*0.000001;
P=10000;
Vhigh=900;
C=0;
s=poly(0,"s");
while Vlow < 450   
       Polynome=((P*Lk*n*fs)./(Vhigh*Vlow))-s+s^2;
solution=roots(Polynome);
out =  solution(solution<0.35 & solution>0.129);
I1=(1./(2*Lk*fs))*(2*(Vlow./n)*out+Vhigh-(Vlow./n));
I2=(1/(2*Lk*fs))*(2*Vhigh*out-Vhigh+(Vlow./n));

 if I1 > I2 then
    C= I1;
 else
     C=I2;
 end
plot2d(Vlow,C,4);
Vlow=Vlow+20;
end