x=10:2:34;y=[9.6 11.5 12.7 13.8 14.4 15.7 16.1 17.2 17.8 19.0 19.8 20.3 21.8];myfunc=inline('a(1)*(1-exp(-a(2).*x)).^a(3)','a','x');a=nlinfit(x,y,myfunc,[25 1 0.05])x1=10:2:34;y1=a(1).*(1-exp(-a(2).*x)).^a(3);plot(x,y,'o',x1,y1)