propforce.f90

来自「动力学计算程序」· F90 代码 · 共 13 行

real function PropForce(t)
!-------------------------------------------------------------
!  Purpose: calculate proportional force at t
!  Input:   t - time
!  Output:  PropForce - value of force at t.
!-------------------------------------------------------------
	use module_parameter, only: a=>prop_para
	implicit none
	real t
	PropForce = a(1) + a(2)*t + a(3)*sin(a(4)*t+a(5))

end function PropForce