Make a program for solving Newton's law of cooling T′=−k(T−Ts),T(0)=T0, t∈(0,tend] with the Forward Euler, Backward Euler, and Crank-Nicolson schemes (or a θ scheme). Verify the implementation.