def suspensionDE(Y, t, m, c, k):
	"""
		Colin Simpson's code...
		But the beautiful docstring is all mine

		Y: array_like
			contains initial conditions
		
		t: array_like
			time values

		m: float
			mass parameter

		c: float
			damping coeffcient parameter
		
		k: float
			spring constant parameter
	"""

	return [Y[1], -c*Y[1]-k*Y[0]/m]