syms x(t) v(t) I(t) F(t) s
v(t) = diff(x, t);
eq = diff(v) == (562/101)*x+(6744/101)*I+(1/50500)*F;
lap = laplace(eq);
S = subs(lap, [x(0) v(0)], [0 0]);
syms lx lv lI lF
Sx = subs(S, [laplace(x(t), t, s) laplace(I(t), t, s) laplace(F(t), t, s)], [lx lI lF]);
Sv = subs(S, [laplace(x(t), t, s) laplace(I(t), t, s) laplace(F(t), t, s)], [lv/s lI lF]);
g11 = simplify(lx/solve(subs(Sx,lF,0),lI))
g12 = simplify(lx/solve(subs(Sx,lI,0),lF))
g21 = simplify(lv/solve(subs(Sv,lF,0),lI))
g22 = simplify(lv/solve(subs(Sv,lI,0),lF))
Best Answer