function x = schubkurbel(t, om, l1, l2) % Ortskurve x(t) der Schubkurbelschwingung % om Kreisfrequenz % l1 Laenge der kuerzeren Schubstange % l2 Laenge der laengeren Schubstange la = l1/l2; % Schubstangenverhaeltnis x = l2*(la*cos(om*t) + sqrt(1 - la^2*sin(om*t).^2));