clear
% SIMPLE CLIMATE MODEL
options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
T0=[15 15];
[t,Temp]=ode45(@scm_fun,[0 100],T0,options);
plot(t,Temp)