ex3b.m

% Schubkurbelschwingung, Phasenraum
l_2 = 1.0; % beliebiger Wert fuer l_2
w = 1.0;   % beliebiger Wert fuer omega

tmax = 2*pi;
t = [0:tmax/400:tmax];

la = 0.2;
x1 = l_2*(la*cos(w*t) + sqrt(1 - la^2*sin(w*t).^2));
v1 = -l_2*w*la*sin(w*t) .* (1 + (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));
v2 = -l_2*w*la*sin(w*t) .* (1 + (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));
v3 = -l_2*w*la*sin(w*t) .* (1 + (la*cos(w*t))./sqrt(1 - la^2*sin(w*t).^2));

plot(x1, v1, x2, v2, x3, v3)
%axis([0 tmax 0.5 1.5]);
ylabel('v');
xlabel('x');
legend('\lambda=0.2', '\lambda=0.3','\lambda=0.4', 0)