function aufgabe_1 (method)
%% Parameter + ABs
p.m1=1;
p.m2=2;
p.c1=500;
p.c2=700;
p.method=method;
x1=0;
x2=0;
x1p=0.2;
x2p=0.3;

%% Lösen der Dgl
[t,X] = ode45(@dgl,[0 6],[x1 x2 x1p x2p],[],p);
%% Plot
fig=figure(1);
subplot(3,1,1);plot(t,h_von_t(t,p.method));
title('Bewegung des Aufhängepunktes'); 
grid on; xlabel('Zeit in [s]'); ylabel('Weg in [m]');
subplot(3,1,2);plot(t,X(:,1:2)); 
legend('x1','x2'); title('Bewegungen der Massen'); 
grid on; xlabel('Zeit in [s]'); ylabel('Weg in [m]');
subplot(3,1,3);plot(t,X(:,3:4));
legend('v1','v2'); title('Geschwindigkeit der Massen'); 
grid on; xlabel('Zeit in [s]'); ylabel('Geschwindigkeit in [m/s]');
filename = p.method;
print(fig, '-dpdf', filename);
end

%% Dgl
function zp=dgl(t,z,p)
h=h_von_t(t,p.method);
zp=[z(3);
    z(4);
    1/p.m1*(-p.c1*(z(1)-h)-p.c2*(z(1)-z(2)));
    -p.c2/p.m2*(z(2)-z(1))];
end

%% Bewegung des Aufhängepunktes
function h=h_von_t(t,method)
switch lower(method)
   case '1b_i'
      h=zeros(1,length(t));
   case '1b_ii'
      h=0.5*sin(8*t);
   case '1b_iii'
      ts=[0    1 2.5 4   6]; 
      hs=[0 -0.1 0.3 0.5 0];
      h=interp1(ts,hs,t,'spline');
   otherwise
      disp('Unknown h(t).')
end
end