Using Levenberg-Marquardt algorithm in the optimization tool box

t = 25; % temperature (C)
T = t + 273; % temperature (K)
i_ph = 3.113; % phase current [A]
e_g = 1.12; % band gap [eV]
c_01 = 170.80; % co-efficient of saturation current [AK ^ (- 3)]
m_1 = 1.00; % diode factor
e = 1.6 * 10 ^ - (19); % electronvolt [j]
k = 1.38 * 10 ^ - (23); % Boltzmann constant [JK ^ (- 1)]
U = (0: 0.001: 0.6); % voltage [V]
Ut = (k * T) / e;
Is = c_01 * T ^ (3) * exp (- ((e_g * e) / (k * T)));
I = i_ph-Is. * (Exp (U ./ (m_1 * Ut)) - 1);
plot (U, I, 'linewidth' , 2.5);
axis ([0,0,6,0,3.5]);
xlabel ( 'voltage U' );
ylabel ( 'current A' );
grid on ;

% how can i use the optimization tool (Levenberg-Marquardt algorithm) to find c_01, using a guess value of c_01 = 870.80
% I have plotted both real and guessed value of c_01