请教各路大神！请问大家，我的两个矩阵怎么相乘啊？我的代码怎么都运行不起。能否请知道的大神帮我看看代码哪里有问题？拜托！！

clear

k1 = 0.289;

k2 = 0.793/5000;

k3 = 0.684/100000;

k4 = 0.15;

k0 = 0.727;

k5 = 0.084;

k6 = 0.8;

k7 = 0.531;

A = [E ES ED FNP E1 ES1 FNP1 S Sum_FNP];

B = [1-k1 k1 0 0 0 0 0 0 0;

k0+k4 1-k0-k4-k2 k2 k2+k4 0 0 0 k0 0;

k3 0 1-k3 0 0 0 0 0 0;

0 0 0 1 0 0 0 0 1;

0 0 0 0 1-k5 k5 0 0 0;

0 0 0 0 k6+k7 1-k6-k7 k6 k7 0;

0 0 0 0 0 0 1 0 1;

0 0 0 0 0 0 0 1-k1-k5 0;

0 0 0 0 0 0 0 0 0];

A(1,1) = 2.96e-3;

A(1,2) = 0;

A(1,3) = 0;

A(1,4) = 0;

A(1,5) = 2.96e-3;

A(1,6) = 0;

A(1,7) = 0;

A(1,8) = 0.03e3;

A(1,9) = 0;

C = A * B;

for i=3:200

A = [E ES ED FNP E1 ES1 FNP1 S Sum_FNP];

B = [1-k1 k1 0 0 0 0 0 0 0;

k0+k4 1-k0-k4-k2 k2 k2+k4 0 0 0 k0 0;

k3 0 1-k3 0 0 0 0 0 0;

0 0 0 1 0 0 0 0 1;

0 0 0 0 1-k5 k5 0 0 0;

0 0 0 0 k6+k7 1-k6-k7 k6 k7 0;

0 0 0 0 0 0 1 0 1;

0 0 0 0 0 0 0 1-k1-k5 0;

0 0 0 0 0 0 0 0 0];

C = A * B;

end

PNS_3mm=[150        0.0082

550            0.0244

950          0.0408

1350        0.0559

1750        0.0715

2150        0.0866

2250        0.0903

2950        0.1157

3350        0.1295

3750        0.1431

4150        0.1565

4550        0.1692

4950        0.1818

5350        0.1943

5750        0.2058

6150        0.2174

6550        0.2287

6950        0.2393

7350        0.25

7750        0.2602

];

t=PNS_3mm(:,1)'*12;

n=PNS_3mm(:,2)'/16*1000;

hold on

plot (t,n,'*k')

%plot(FNP(1,1:12001),'DisplayName','FNP(1,1:12001)','YDataSource','FNP(1,1:12001)');figure(gcf)

plot(Sum_FNP(1,:),'r','LineWidth',1.5);

hold on

以上是我的代码，k值是已知的，数据是time--Sum_FNP，现在我代码运行不起，不知道哪里的原因！

请知道的大神帮帮忙！