遇到一个最小二乘法求曲线拟合的程序，但是是fortran77版编译的，看不懂，求大神用do和enddo给改一下，感谢！！
subroutine hpir1(x,y,a,n,m,dt1,dt2,dt3)
dimension x(n),y(n),a(m),s(20),t(20),b(20)
double precision x,y,a,s,t,b,dt1,dt2,dt3,z,d1,p,c,d2,g,q,dt
 do 5 i=1,m
 5 a(i)=0.0
if (m.gt.n) m=n
if (m.gt.20) m=20
 z=0.0
 do 10 i=1,n
10 z=z+x(i)/n
b(1)=1.0
d1=n
p=0.0
c=0.0
 do 20 i=1,n
 p=p+(x(i)-z)
 c=c+y(i)
20 continue
c=c/d1
p=p/d1
a(1)=c*b(1)
if (m.gt.1) then
t(2)=1.0
t(1)=-p
d2=0.0
c=0.0
g=0.0
 do 30 i=1,n
 q=x(i)-z-p
 d2=d2+q*q
 c=y(i)*q+c
 g=(x(i)-z)*q*q+g
 30 continue
c=c/d2
p=g/d2
q=d2/d1
d1=d2
a(2)=c*t(2)
a(1)=c*t(1)+a(1)
endif
do 100 j=3,m
s(j)=t(j-1)
s(j-1)=-p*t(j-1)+t(j-2)
if (j.ge.4) then
do 40  k=j-2,2,-1
40 s(k)=-p*t(k)+t(k-1)-q*b(k)
endif
s(1)=-p*T(1)-q*b(1)
d2=0.0
c=0.0
g=0.0
do 70 i=1,n
q=s(j)
do 60 k=j-1,1,-1
60 q=q*(x(i)-z)+s(k)
d2=d2+q*q
c=y(i)*q+c
g=(x(i)-z)*q*q+g
70 continue
c=c/d2
p=g/d2
q=d2/d1
d1=d2
a(j)=c*s(j)
t(j)=s(j)
do 80 k=j-1,1,-1
a(k)=c*s(k)+a(k)
b(k)=t(k)
t(k)=s(k)
80 continue
100 continue
dt1=0.0
dt2=0.0
dt3=0.0
do 120 i=1,n
q=a(m)
do 110 k=m-1,1,-1
110 q=q*(x(i)-z)+a(k)
dt=q-y(i)
if (abs(dt).gt.dt3) dt3=abs(dt)
dt1=dt1+dt*dt
dt2=dt2+abs(dt)
120 continue
return
end
把循环用do enddo表示就好，谢谢！