function dy = state_equations(y, f, M)
% allapot-egyenlet, metastabil periodikus forgohullamok
% 2013. 12. 10.

    dy = M * f(y) - y;

end
