[laosam280]

x=[0.7 0.1312 0.02685 0.00943 0.0033];
y=[0.3636 0.1682 0.4741 1.2704 2.07];
plot(x,y,'g-'); hold on
x0=[1 0]';
x=lsqnonlin('myfun',x0)
y1=y(1);
y=a*x.^b;
plot(x,y,'bo-');
%define function myfun,saved at work
function f=myfun(x,y)
x=[0.7 0.1312 0.02685 0.00943 0.0033];
y=[0.3636 0.1682 0.4741 1.2704 2.07];
y1=y(1);
z=y1-x(1)*x.^x(2);
f=z-y;