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)
Ar=solve('(1-3pi/16)piG^(D-1)f1Ar^(D/2)(f2Ar^((3-2D)/2)-ac^((3-2D)/2))+KSyf3Ar^(D/2)*ac^((2-D)/2)-p','Ar')
建议用vpasolve试一试