function xd = ssodefunc(t,x,a,b,c,u)
    xd = a*x + b*u;
end