clc;
clear;
Vocs=21.6;
I=0:0.01:2;
k=1.38*10^-23;
q=1.6*10^-19;
Tr=301.18;
Pmax=25;
Vmax=18;
Imax=1.39;
Isc=1.83;
Kt=2.06*10^-2;
K=-0.77;
A=1;
T=298;
G=1000;
%for T=293:5:313
%for T=298
%for G=1000:-200:200
Iph=(Isc+Kt*(T-298))*G/1000;
Vt=A*k*T/q;
Voc=Vocs+K*(T-298);
Io=Iph/(exp(Voc/(50*Vt))-1);
Rs=((50*Vt.*log(((Iph-Imax)/Io)+1))-Vmax)/Imax;
V=50*Vt.*(log((Iph-I+Io)/Io))-I*Rs;
figure(1);
plot(V(1:184),I(1:184));
hold on;
grid;
P=V.*I;
figure(2);
plot(V(1:184),P(1:184));
hold on;
grid;
%end
%end