open (10,file="verlet_sumb.in")
	read (10,*) nt,dt,x,v,rm,rk,beta
	omega2=rk/rm
	call kiirendus(t,x,v,a,omega2,beta)
	do i=1,nt
	x=x+v*dt+a*dt**2/2.
	call kiirendus(t,x,v,a1,omega2,beta)
	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
!	print *,i*dt,x
	enddo
	stop
	end
	
	
	
	subroutine kiirendus(t,x,v,a,omega2,beta)
	a=-omega2*x-beta*v
	return
	end