module calculations
public::harmonic
contains
subroutine harmonic(t,p,v,ac)
real,intent(in)::t
real,intent(out)::p,v,ac
real,parameter::n=3.14159265,eps=0,a=2.5
p= a*sin(n*t+eps)
v= n*a*cos(n*t+eps)
ac= -a*(n**2)*sin(n*t+eps)
end subroutine harmonic
end module calculations