clear
close
%givens
E0 = 3; %input voltage
ls = 7*10^-14; %value of ls
R = 200; %value of resistor
C = .026; %value of KT/q
V_vec = [];
Vdiode(1) = 1; %initial guess for Vdiode
Vdiode(2)= 2;
i=1;
%iterate until difference between last 2 Vdiodes < 10^-6
while abs(Vdiode(i + 1) - Vdiode(i)) > 10^-6
% derivative of f
fprime = -1-ls*R/C * exp(Vdiode(i)/C);
% f at Vdiode
f = E0 - Vdiode(i) - ls * R * (exp(Vdiode(i)/C) - 1);
% new Vdiode
Vdiode(i) = Vdiode(i) - f/fprime;
% store in vector
V_vec(i) = Vdiode(i);
% increment i
i = i + 1;
end
Vdiode
length(V_vec)