function [yp] = heiluri(t,Y)
% Heiluri-systeemin m"a"arittely
%   
x=Y(1);y=Y(2);
g=9.81;L=g;
yp=[y;
    -(g/L)*sin(x)];
end

