rge in coulombs</span>
nVth = n*Ns*k*(Tcell+273.15)/q;
NumPoints = 1000;
[IVResult] = pvl_singlediode(IL, I0, Rs, Rsh, nVth, NumPoints);

figure
plot(IVResult.V,IVResult.I)
hold <span class="string">on</span>
scatter(IVResult.Vmp,IVResult.Imp,<span class="string">'filled'</span>)
xlabel(<span class="string">'Voltage (V)'</span>)
ylabel(<span class="string">'Current (A)'</span>)
title(<span class="string">'Example IV Curve from Single Diode Model'</span>,<span class="string">'FontSize'</span>,14)
Pmax = IVRe