real y(2),work(15)
	integer iwork(5)
	common omega2,beta
	external func
	open (10,file="vaba_rkf.in")
	read (10,*) nt,dt,omega2,relerr,abserr,y,beta
	
	iflag=0
	t=0.
	tout=0.
	call rkf45(func,2,y,t,tout,relerr,abserr,iflag,work,iwork)
	
	do i=1,nt
	iflag=2
	tout=t+dt
	call rkf45(func,2,y,t,tout,relerr,abserr,iflag,work,iwork)
	t=tout
	print*,iflag
	print *,t,y
	enddo
	stop
	end
	subroutine func(t,y,dy)
	common omega2,beta
	real y(2),dy(2)
	dy(1)=y(2)
	dy(2)=-omega2*y(1)-beta*y(2)
	return
	end