function [timeode45 ode45mg] = rugo_numgyakhoz(t_null, t_max, x_null, y_null, b)

rugo_ode=@(t,xy) [xy(2); -xy(1)];
%rugo_ode=@(t, xy) [xy(2); -b*xy(2)-xy(1)]; %csillapított

[timeode45, ode45mg]=ode45(rugo_ode, [t_null, t_max], [x_null, y_null]);

end