clear all, close all, clc;

p=[10 16 25 40 60]';
q=[94 118 147 180 230]';
A=[ones(size(q)) log(p)];
c=A\log(q);
c1=exp(c(1))
c2=c(2)

kortP=p(1):p(end);
kortQ=c1*kortP.^c2;
plot(kortP,kortQ,p,q,'o')