open (10,file="reson_verlet.in")
	read (10,*) nt,dt
	read (10,*) x1,v1,rm1,rk1
	read (10,*) x2,v2,rm2,rk2
	read (10,*) rk
	om2_1=rk1/rm1
	om2_2=rk2/rm2
	om2=rk
	call kiirendus(t,x1,x2,v1,v2,a1,a2,om2_1,om2_2,om2)
	do i=1,nt
	x1=x1+v1*dt+a1*dt**2/2.
	x2=x2+v2*dt+a2*dt**2/2.
	call kiirendus(t,x1,x2,v1,v2,a1_1,a2_1,om2_1,om2_2,om2)
	v1=v1+(a1_1+a1)*dt/2.
	v2=v2+(a2_1+a2)*dt/2.
	a1=a1_1
	a2=a2_1
	print *,i*dt,x1,x2
!	print *,i*dt,x
	enddo
	stop
	end
	
	
	
	subroutine kiirendus(t,x1,x2,v1,v2,a1,a2,om2_1,om2_2,om2)
	a1=-om2_1*x1+om2*(x2-x1)
	a2=-om2_2*x2+om2*(x1-x2)
	return
	end