%EXS 4.14 
p= [10 16 25 40 60]'
q = [94 118 147 180 230]'
A = [ones(size(p)) log(p)]
b = log(q)
x= A\b
C=exp(x(1))
k=x(2)

%residual
residual=q-C*p.^k
r = norm(q-C*p.^k)/norm(q)

%Plot
P =10:60;
funk_q=C*P.^k;
plot(P, funk_q, p,q,'*')
