Chapter 9, Problem 32P

Use a MATLAB built-in funtion to numerically solve:

$\frac{dy}{dx}=-x^2+\frac{x^3{e}^{-y}}{4}$ for $1\le x\le 5$ with y(1)=1

Plot the solution.