% grxxm013.m (Aufgabe 1.3)

% Gedeampfter Massenschwinger mit zwei gekoppelten Massen
% Massen 0 < m1, m2,

T=10;
m1=2;
m2=3;
d=2;
k1=0;
k2=0;
A=[0,0,1,0;0,0,0,1;-d/m1,d/m1,-(k1+k2)/m1,k2/m1;...
    d/m2,-d/m2,k2/m2,-k2/m2];

F=@(t,y) A*y;

x1=1;
v1=0;

x2=2;
v2=0;

y0=[x1;x2;v1;v2];
tspan=[0,T];
opts=odeset('reltol',1.e-6);
ode23(F,tspan,y0,opts);
