rendered paste body> restart;
> Iph := 7.0*10^(-3);
> I0 := 1*10^(-14);
> Rsh := 1000;
> Rs := 10;
> n := 1;
>
> alpha := Rsh*(Rs*Iph+Rs*I0+V)/n/0.02585/(Rs+Rsh):
> i := -V/Rs+n*0.02585/Rs*( -LambertW(Rs*I0*Rsh/n/0.02585/(Rs+Rsh)*exp(alpha)) + alpha ) :
>
> N:=2:
> for j from -1*N by 1 while j <= 10*N do
> a[j]:=subs (V=j/10/N,i)*1000;
> od:
> l := [[ t/10/N, -a[t]] $t=-1*N..10*N]:
> plot(l, x, y=-Iph*1000*1.1..Iph*1000*0.5, style=point,symbol=circle):
>
> FF:=0:
> for j from -1*N by 1 while j <= 10*N do
> if (a[j] > 0) then Voc:=(a[j]*(j+1)-a[j+1]*j)/(a[j]-a[j+1])/10/N fi;
> if (j/10/N*a[j] > FF) then FF:=j/10/N*a[j] fi
> od;
> Digits := 5:
> Voc:=evalf(Voc);
> Jsc:=a[0];
> FF:=FF/Jsc/Voc;
> eta:=FF*Jsc*Voc;