P（x）为 A*t(s（x）)在0到x范围对x的变上限积分,s(x)为P（x）/B在0到x范围对x的变上限积分，t（s（x））=-3.11302293666198e-5*(0.00331486138220175*(2*s（x） - 0.487)*(s（x）**2 - 0.487*s（x） + 17.3687)*(14.1924*s（x）**3 + 471.841*s（x）**2 + 1048.2*s（x） + 26.3996) - (0.0575748328890475*s（x）**2 - 0.0280389436169662*s（x） + 1)**2*(42.5772*s（x）**2 + 943.682*s（x） + 1048.2))*(14.1924*s（x）**3 + 471.841*s（x）**2 + 1048.2*s（x） + 26.3996)/(0.0575748328890475*s（x）**2 - 0.0280389436169662*s（x） + 1)**4，A = 3.14*0.18，B = (2e5*3.14*0.01*0.01*3)，牛顿迭代法求解P（x）并绘图。