clear
format long

% input intial data
S=[.01 .05 .1 .5 1 5 10 50 100];
v0= [6.078e-11 7.595e-9 6.063e-8 5.788e-6 1.737e-5 2.423e-5 2.43e-5 2.431e-5 2.431e-5];

% first create an m-file function to compute the sum of the squares- fSSR()

% the minimization of the function is implemented by 
a= fminsearch(@fSSR, [2e-5, 1], [], S, v0)

% Therefore the best fit model is 
v0p= a(1)*S.^3./(a(2)+S.^3);

% Graphically display the fit along the data using log-log plot due to  the
% wide variation of the magnitudes of the values being displayed.
loglog(S, v0, 'o', S, v0p)