importnumpyasnpimportmatplotlib.pyplotasplt初始化参数a=1-0.2*(1/12)m=1.109*10**5w3=17.86w4=22.99X=[]Z=[]计算X和Z的值forkinnp.arange(0,0.875,0.001):x1=1.221011*(1-1/(m*(1-a-0.42*k)8*(1-a)24*(0.5+