disp（‘牛頓下山法‘）?
p0=0.3;?
for?i=1:n0?????
????k=1;
????p=p0-k*（600*（p0.^4）-550*（p0.^3）+200*（p0.^2）-20*p0-1）/（2400*（p0.^3）-1650*p0.^2+400*p0-20）;?
????while（abs（600*（p1.^4）-550*（p1.^3）+200*（p1.^2）-20*p1-1）>=abs（600*（p0.^4）-550*（p0.^3）+200*（p0.^2）-20*p0-1）&&k>tol）
????????k=k/2;
????????p1=p0-k*（600*（p0.^4）-550*（p0.^3）+200*（p0.^2）-