点击查看代码importnumpyasnpimportpylabasplta=1-0.2**(1/12);m=1.109*10**5w3=17.86;w4=22.99X=[];Z=[];N=[]forkinnp.arange(0,1.0,0.001):x1=1.22*10**11*(1-1/(m*(1-a-0.42*k)**8*(1-a)**24*(1/2+(1-a-k)**8*(1-a)**4