function plot_temp()

[L,K,T] = av_temp
plot(L,T(:,1),'g',L,T(:,2),'r',L,T(:,3),'b',L,T(:,4),'k')
  title(' The temperature as a function of latitude ')
  xlabel(' Latitude (in degrees) ')
  ylabel(' Temperature ')
  legend('K=.67','K=1.5','K=2.0','K=3.0')

function [L,K,T] = av_temp
% Q&S Table 3.1
L = 65:-10:-55;
K = [.67, 1.5,2.0,3.0];
% T is 13x4
T = [-3.10, 3.52, 6.05, 9.30;...
     -3.22, 3.62, 6.02, 9.30;...
     -3.30, 3.65, 5.92, 9.17;...
     -3.32, 3.52, 5.70, 8.82;...
     -3.17, 3.47, 5.30, 8.10;...
     -3.07, 3.25, 5.02, 7.52;...
     -3.02, 3.15, 4.95, 7.30;...
     -3.02, 3.15, 4.97, 7.35;...
     -3.12, 3.20, 5.07, 7.62;...
     -3.20, 3.27, 5.35, 8.22;...
     -3.35, 3.52, 5.62, 8.80;...
     -3.37, 3.70, 5.95, 9.25;...
     -3.25, 3.70, 6.10, 9.50];