E=1.310^11
D=1.69
G=110^(-10)
K=2.8
Sy=310^8
n=0.072
a=0.421
v=310^(-3)
h=110^(-5)
A=1
ac=G^2(2E/K/Sy)^(2/(D-1))
f1=(D/(3-2D))((2-D)/D)^(D/2)
f2=((2-D)/D)^((3-2D)/2)
f3=(D/(2-D))^((2-D)/2)
P=(0:0.510^7:310^7)
P=(1-3pi/16)pi.G.^(D-1)f1Ar.^(D/2).(f2.Ar.^((3-2D)/2)-ac.^((3-2D)/2))+KSyf3.Ar.^(D/2)ac.^((2-D)/2)
al=(2-D)Ar/D
Kl=9(A-Ar).^3nav/2./Ar/pi/h.^4
figure(1)
hold on
plot (P,Kl)
xlabel ('P')
ylabel ('Kl')