[y708321017]

x=[0 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90];
y1=[0.3 0.983 0.984 0.983 0.986 0.985 0.984 0.982 0.985 0.986 0.983 0.982 0.984 0.983 0.985 0.985 0.984 0.982 0.988 0.981 0.984 0.983 0.981 0.983 0.986 0.985 098.7 0.984 0.985 0.986 0.982];
y=872.6*y1* 0.00652/(8.31*278);
myfunc=inline('-(beta(1)-beta(3)*exp(beta(2)*x))/beta(2)','beta','x');
beta=nlinfit(x,y,myfunc,[01 -0.1 0.1]);
a=beta(1),b=beta(2),c=beta(3)


Warning: The Jacobian at the solution is ill-conditioned, and some model parameters may not be estimated well (they are not identifiable). Use caution in making predictions.
> In nlinfit at 223

a = 0.0148
b = -1.4162
c = 0.0137