function [ f ] = mittsystem_f( t,y )
% Funktionsfil motsvarande BVP: y'(t)=f(t,y(t))

    f1= 3*y(1)+y(2)*y(2);
    f2= y(1)-4*y(2);
    f=[f1;f2];
end

