function solveshit()
k = 0.01;
sol = ode45(@(t,y) dA(t,y,k),[0 20],0);
plot(sol.x,sol.y,'k-');

function dy = dA(t,y,k)
dy =k*exp(y/2)*exp(-t/3);