% derivation of 3 mass system
clear all; close all; clc;

syms F0 m0 x0 c1 d1 m1 x1 c2 d2 m2 x2 Fe s

x1 = (-m2*s^2+d2*s+c2)/(d2*s+c2)*x2;
x0 = (F0+(d1*s+c1)*x1)/(-m0*s^2+d1*s+c1);
F0 = -m1*s^2*x1-m0*s^2*x0-c2*(x2-x1)-d2*(x2-x1);

solve(F0)

