function S=Schockley_funktie(I_s,V,V_T);

n=1;

S=I_s.*(exp(V./(n.*V_T))-1);
