... while(p>z[j-1 ]-0 .05) { dx=x[j-1]-x[j]; w[j-1]=(b[j-1]+(z[j-1]-a[j-1])/m[j-1])*(z[j-1]-a[j-1])/m[j-1]; s[j-1]=2*b[j-1]+2*(z[j-1]-a[j-1])/m[j-1]+2*(z[j-1]-a[j-1])*sqrt(pow(m[j-1],2)+1)/m[j-1]; r[j-1]=w[j-1]/s[j-1]; ... r[j-1]=w[j-1]/s[j-1]; c[j-1]=pow(r[j-1],y)/n[j-1]; v[j-1]=q[j-1]/w[j-1]; p=z[j]+(pow(v[j],2)-pow(v[j-1],2))/19.6+dx/2*(pow(v[j-1],2)/pow(c[j-1],2)/r[j-1]+pow(v[j],2)/pow(c[j],2)/r[j]); z[j-1]=z[j-1]+0.01; ... jp+1. T-ơng tự nh vậy lần lợt tiến hành đánh số các mặt cắt cho những dòngnhập lu tiếp theo. 12 vj-12/2g J J.lxj-1 hj-1 - hj vj2/2g hj-1 Zj-1 hj = Zj Jd lj-1...