open (10,file="verlet.in")
	read (10,*) nt,dt,x,v,rm,rk
	omega2=rk/rm
	call kiirendus(t,x,a,omega2)
	do i=1,nt
	x=x+v*dt+a*dt**2/2.
	call kiirendus(t,x,a1,omega2)
	v=v+(a1+a)*dt/2.
	a=a1
	print *,i*dt,x,v,rm*v**2/2.,rk*x**2/2.,(rm*v**2+rk*x**2)*0.5
	enddo
	stop
	end
	
	
	
	subroutine kiirendus(t,x,a,omega2)
	a=-omega2*x
	return
	end