function [t,y]=finale03
  
tspan = [0:1:79]*365;
x0 = [0;0;0];

% solve the ode
[t,y] = ode45(@foo,tspan,x0);

function dx = foo(t,x)
  A = [-0.0361 0.0124 0.000035;
        0.0111 -0.0286 0;
        0.0039 0 -0.000035];
  
  f = [49.3;0;0];
  
  dx = A*x + f;