myfun = inline('[1./x x.^2]');
x=linspace(0.1,20,100)';
c0=[1;-1];
y=myfun(x)*c0+randn(size(x))/10;
[c,R2,rout] = normalfit(x,y,myfun);
c
R2
figure(1)
plot(x,1./y,'ro',x,1./(myfun(x)*c),'b-')
figure(2)
plot(x,rout,'bo-')