%Aufgabe 4d

%Berechnungsparameter
a=0;
b=5;
h=0.6;

%Parameter
cr=90000;
ca=20000;
dr=100;
da=1000;
mr=20;
ma=250;
yboden=0.1;
yBodenabl=0;

%Anfangswerte
y0=0.1;

%Anfangswertproblem


%Berechnung Konstanten


%Funktionen
%y1=y(2);                        =y(1)
%y2=-ca/ma.*(y(1)-y(3);                 =y2'
%y3=y(4);                        =y3'
%y4=-cr/mr.*(y(3)-yBoden)-dr/mr.*(y(4)-yBodenabl)+ca/mr(y(1)-y(3)+da/mr.*(y(2)-y(4);                               =y(4)


F=@(t,y)  [-ca/ma.*(y(1)-y(3)),y(4),-cr/mr.*(y(3)-yBoden)-dr/mr.*(y(4)-yBodenabl)+ca/mr.*(y(1)-y(3))+da/mr.*(y(2)-y(4))];

[t,y]= euler_alex(a,b,y0,h,F)
y1=y(:,1);
y2=y(:,2);
y3=y(:,3);
y4=y(:,4);


%Diagramm
%subplot(1,2,1)
plot(t,y1,'c' ,t,y2,'b', t,y3,'g',t,y4,'r') 
grid on;
axis([-100 100 -100 100]);

xlabel('Zeit t');
ylabel('y-Wert');
title('Dämpfer');