这段程序画Kl随P变化图像,P是变量P公式中有Ar,Kl公式中也有Ar,为什么它说识别不出来Ar

E=1.3*10^11
D=1.69
G=1*10^(-10)
K=2.8
Sy=3*10^8
A=1
ac=G^2*(2*E/K/Sy)^(2/(D-1))
al=(2-D)*Ar/D
f1=(D/(3-2*D))*((2-D)/D)^(D/2)
f2=((2-D)/D)^((3-2*D)/2)
f3=(D/(2-D))^((2-D)/2)
P=(0:0.5*10^7:3*10^7)
P=(1-3*pi/16)*pi*G^(D-1)*f1*Ar^(D/2)*(f2*Ar^((3-2*D)/2)-ac^((3-2*D)/2))+K*Sy*f3*Ar^(D/2)*ac^((2-D)/2)
Kl=9*(A-Ar).^3*n*a*v/2./Ar/pi/h.^4
figure(1)
hold on
plot (P,Kl)
xlabel ('P')
ylabel ('Kl')

 

img

第八行 al=(2-D)*Ar/D 用到了 Ar,但是在第八行之前没有定义 Ar