function xdot=esim_ode_fun(t,x,u)
% ESIM_ODE_FUN, odefun-tiedostoesimerkki
%    dx/dt(t) = -0.1*x(t) + 3*u(t)

% differentiaaliyhtälö (voi olla myös pystyvektori)
xdot=-0.4*x+3*u;