clc
clear all
pai0=[11.9900000000000,35.6000000000000,6.75666666700000,26.0100000000000,7.05000000000000,6.04333333300000,11.2350000000000,31.0600000000000,8.50500000000000,24.4100000000000,11.3850000000000,21.7700000000000,6.47333333300000,11.8400000000000,33.8500000000000,5.02000000000000,9.10000000000000,28.9700000000000,9.05000000000000,26.8100000000000,9.34000000000000,14.1200000000000,5.35000000000000,17.4900000000000,7.35000000000000]'
pai1=[1250;2500;833.333333333333;400;200;666.666666666667;1000;2000;250;350;350;250;750;1125;2250;566.666666666667;850;1700;750;1500;300;150;150;200;933.333333333333]'
pai2=[0.00422591402541550;0.00831629085077855;0.00283607074943188;0.0105409578681139;0.00531613445949953;0.00255861859428642;0.00380625907479581;0.00751483752544705;0.00591948718845710;0.00987925113736621;0.00695654193723514;0.00838304550479858;0.00270055666443629;0.00401948344303690;0.00792344989567206;0.00238263128531228;0.00354879534293637;0.00697061153384213;0.00334205278018335;0.00658417615288411;0.00646730257093456;0.00655067556132041;0.00462822217809494;0.00753093361753438;0.00298476709202631]'
kesi0=log(pai0)
kesi1=log(pai1)
kesi2=log(pai2)
y=kesi0;
M=size(pai0,1)
A=zeros(M,3)
A

,1)=1
A

,2)=kesi1
A

,3)=kesi2
x=(inv(A'.*A))*A'*y
C=exp(x(1,1));
alpha=x(2,1);
beta=x(3,1);
pai11=pai1.^alpha;
pai22=pai2.^beta
pai0sec=C.*(pai11).*pai22.*pai33
ضرب ماتریس ها رو نمیدونم چه جوری باید با حلقه for بنویسم میشه کمک کنید ؟؟