schubkurbel.m

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));