задача об изгибе консоли (задача Коши) для fortran77

real M
M(x)=1/sqrt((1+(0.02*(13+6)*x)**2)**3)
f(x,z)=0.02*(13+6)*M(x)*sqrt((1+z**2)**3)
open(1,file=dmtr4.dat')
print*,'vvesti dl,N,x0,y0,z0'
read*,dl,N,x0,yo,zo
h=dl/N
print*,'dl=',dl,'N=',N,'x0=',x0,'y0=',y0,'z0=',z0,'h=',h
write(1,*)'dl=',dl,'N=',N,'x0=',x0,'y0=',y0,'z0=',z0,'h=',h
print100
write(1,100)
100format(7x,'x',10x,'y',10x,'z'/)
do 2 i=0,N
print 3,x0,y0,z0
3format(3E12.4)
write(1,3)x0,y0,z0
z=z0+h*(z0+z)/2
x0=x0+h
z0=z
2y0=y
end

ЛАБА3!!!
parameter(n1=8,n2=6)
dimension u(0:n1,0:n2)
open(1,file='dmtr3.dat')
print*,'vvesti dl1,dl2,c='
read*,dl1,dl2,c
print*,'dl1=',dl1,'dl2=',dl2,'c=',c
write(1,*)'dl1=',dl1,'dl2=',dl2,'c=',c
h1=dl1/n1
h2=dl2/n2
do 1 i=0,n1
x=h1*j
do 1 j=0,n2
y=h2*i
u(i,j)=0
if(i.eq.0.or.i.eq.n1)u(i,j)=4*13/dl2*y*(dl2-y)
1if(j.eq.0.or.j.eq.n2)u(i,j)=4*6/dl1*x*(dl1-x)
k=1
2z=0
do 3 i=1,n1-1
do 3 j=1,n2-1
u(i,j)=((i-1,j)+u(j+1,j))/h1**2+(u(i,j-1)+u(i,j+1))/h2**2-c)/(2/h1**2+2/h2**2)
3z=z+abs(u(i,j)-r)
k=k+1
if(z.gt.0.001.and.k.lt.100) goto 2
print 4,k
write(1,4)k
4format(2x,'chislo iteracii k=',i3//5x,'reshenie'//)
do 5 j=0,n2
write(1,6)(u(i,j),i=0,n1)
5print 6,(u(i,j),i=0,n1)
6format(11f7.2)
end