% shell0=[,];
% tdata=[];
% rdata=[];
% 
% [shell,resnorm]=lsqcurvefit(@hoffsolve,shell0,tdata,rdata);
% 
clear all

shell(1)=100; %shear modulus Gs in MPa assuming a 1nm shell and Hoff model
shell(2)=1;   %shear viscosity in Pas assuming a 1nm shell and Hoff model

[a,b]=hoffsolve(shell);

%figure
plot(a,b)

coat=shell*.99;
lb=coat*.8;
ub=coat*1.2;


%[x]=hoffrun(coat,a);

[shellfit,resnorm]=lsqcurvefit(@hoffrun,coat,a,b,lb,ub)