function dxy = rugo_egyenlet(t,xy, b)
dxy = zeros(2,1); %kimenet [dx, dy]
x=xy(1);
y=xy(2);
%egyenletek: rugó egyenlet:        dx=y; dy=-x;
%            csillapított rugó egyenlet:  dx=y; dy=-b*y-x;

%dxy(1) = y;
%dxy(2) = -x;

dxy(1)=y; 
dxy(2)=-b*y-x;
end