program harmonic_motion
use calculations
real::t,position,velocity,acceleration
print*,"Please enter the time:"
read*,t
call harmonic(t,position,velocity,acceleration)
print*,"position=",position
print*,"velocity=",velocity
print*,"acceleration=",acceleration
end program harmonic_motion