clc
clear all
v=[3.73 4.17 4.97 6.06 6.71 7.17 7.51 7.98 8.67 9.39 9.89]
dp=[30 35.5 50.5 75.0 92.0 105.0 115 130 153.5 180 199.5]
n=length(v);
plot(dp,v)
y=log(v);
x=log(dp);
A=[n sum(x)
sum(x) sum(x.^2)];
B=[sum(y) sum(x.*y)]';
sol=inv(A)*B;
b=sol(2);
a=exp(sol(1));
dp;
v1=a*dp.^b
hold on
plot(dp,v1,'o')