dydt=@(y,t) (y+3*t)/t^2;              % funcao anonima
tspan = [1 4];                        % vetor intervalo de integracao
y0 = -2;                              % condicao inicial 
[tout,yout] = ode45(dydt,tspan,y0);   % resolve o problema
plot(tout,yout)