Make a program for solving Newton's law of cooling T' = -k(T-T_s),\quad T(0)=T_0,\ t\in (0,t_{\mbox{end}}] with the Forward Euler, Backward Euler, and Crank-Nicolson schemes (or a \theta scheme). Verify the implementation.