% Schubkurbelschwingung, Ortskurve
l_2 = 1.0; % beliebiger Wert fuer l_2
w = 1.0; % beliebiger Wert fuer omega
tmax = 4*pi;
t = [0:tmax/400:tmax];
la = 0.2;
x1 = l_2*(la*cos(w*t) + sqrt(1 - la^2*sin(w*t).^2));
la = 0.3;
x2 = l_2*(la*cos(w*t) + sqrt(1 - la^2*sin(w*t).^2));
la = 0.4;
x3 = l_2*(la*cos(w*t) + sqrt(1 - la^2*sin(w*t).^2));
plot(t, x1, t, x2, t, x3);
axis([0 tmax 0.5 1.5]);
ylabel('x(t)');
xlabel('t');
legend('\lambda=0.2', '\lambda=0.3','\lambda=0.4', 0)