PV characteristic IV curve plotting

Here's a code to plot the characteristic IV curve of PV.

The function is used to calculate Isc and Voc in addition to plotting the characteristic IV curve.
The functions takes the cell's parameters as inputs. The performance of the solar cell and its characteristic curves are determined by the cell's parameters.
These Parameters are:
- short circuit current density(J_sc), reverse saturation current density (J_o), Area (A), and ideality factor(If), Temparature in kelvin (T), Irradiance (G),series resistance(R_s) and shunt resistance(R_sh).
An example;
J_sc=34.3*10^-3;
J_o=10*10^-10;
A=pi*(3*2.5)^2;
G=[1000,800,600];
R_sh=100;
R_s=0;
T=298;
If=1.5;
Plot_IV(J_sc,J_o,A,T,If,G,R_sh,R_s)