program hubble
!
vx1=0
vx2=1
vx3=2
vx4=0
vx5=0
vy1=0
vy2=0
vy3=0
vy4=1
vy5=2
do nt=1,3
  x1=0
  x2=vx2*nt
  x3=vx3*nt
  x4=0
  x5=0
  y1=0
  y2=0
  y3=0
  y4=vy4*nt
  y5=vy5*nt
!
  d20=x2
  d30=x3
  d40=y4
  d50=y5
!
  d23=x3-x2
  d45=y5-y4
!
  d24=sqrt(x2**2+y4**2)
  d34=sqrt(x3**2+y4**2)
  d25=sqrt(x2**2+y5**2)
  d35=sqrt(x3**2+y5**2)
!
write(6,100) d20,vx2
write(6,100) d30,vx3
write(6,100) d23,vx3-vx2
!
write(6,100) d40,vy4
write(6,100) d50,vy5
write(6,100) d45,vy5-vy4
!
d24=sqrt(x2**2+y4**2)
d34=sqrt(x3**2+y4**2)
d25=sqrt(x2**2+y5**2)
d35=sqrt(x3**2+y5**2)
!
100 format(1p3e12.5)
!
end do
!
end