In MATLAB

A ball at 1200K is allowed to cool down in air at an ambient temperature of 300K. Assuming heat is lost only due to radiation, the differential equation for the temperature of the ball is given by d theta/dt = - 2.2067 times 10^-12 (theta^4 - 81 times 10^8), theta (0) = 1200K Find the temperature at t = 480, using (a) Euler (b) Heun (c) Midpoint (d) Ralston (e) Runge-Kutta 4th order method