def Um(m): U = sympy_utils.eye(N) U[1,0] = m*s**2 return U
def Uf(F): U = sympy_utils.eye(N) U[1,2] = -F return U
def Us(k): U = sympy_utils.eye(N) U[0,1] = 1/k return U