>> %where are we?
>> pwd
ans =
/Applications/MATLAB6p5
>> %change to directory with data
>> cd ~/Documents/Courses/2003_04/CIS401-Matlab/Lecture03/Tutorial2/
>> %see what's there
>> ls
ans =
Omaha.html      corn.txt        omaha.txt       omaha2.xls
Tutorial2.doc   omaha.mat       omaha.xls       tut1diary.txt

>> %load omaha.txt
>> load omaha.txt
>> %see what we have
>> whos omaha
  Name        Size                   Bytes  Class

  omaha      53x20                    8480  double array

Grand total is 1060 elements using 8480 bytes

>> omaha(1:5,1:5)
ans =
   1.0e+03 *
    1.9830    1.9840    1.9850    1.9860    1.9870
       NaN    0.0031    0.0026    0.0024    0.0014
       NaN    0.0032    0.0026    0.0024    0.0015
    0.0024    0.0031    0.0026    0.0024    0.0015
    0.0025    0.0031    0.0026    0.0023    0.0014
>> %extract years
>> year=omaha(1,:);
>> omaha=omaha(2:53,:);
>> whos
  Name        Size                   Bytes  Class

  ans         5x5                      200  double array
  omaha      52x20                    8320  double array
  year        1x20                     160  double array

Grand total is 1085 elements using 8680 bytes

>> %test (:)
>> arr2d=[1 2 3; 4 5 6]
arr2d =
     1     2     3
     4     5     6
>> arr2d(:)
ans =
     1
     4
     2
     5
     3
     6
>> %notice how column2 is below 1 and column 3 is after 2
>> %flatten omaha
>> omahav=omaha(:);
>> %look at first and last 5 values of omahav
>> omahav([1:5,end-4:end])
ans =
       NaN
       NaN
    2.4400
    2.5000
    2.5600
       NaN
       NaN
       NaN
       NaN
       NaN
>> 
>> %look at first and last 5 values of omaha
>> omaha(1:5,1)
ans =
       NaN
       NaN
    2.4400
    2.5000
    2.5600
>> omaha(end-4:end,end)
ans =
   NaN
   NaN
   NaN
   NaN
   NaN
>> %looks good
>> %get sizes
>> [m,n]=size(omaha)
m =
    52
n =
    20
>> %create YR and weeks
>> YR=ones(m,1)*year;
>> weeks=(0:m-1)';%note transpose
>> WK=weeks*ones(1,n);
>> %merge year and week
>> YR=YR+WK/m;
>> %look at first and last 5 values of YR
>> YR(1:5,1)
ans =
   1.0e+03 *
    1.9830
    1.9830
    1.9830
    1.9831
    1.9831
>> YR(end-4:end,end)
ans =
   1.0e+03 *
    2.0029
    2.0029
    2.0029
    2.0030
    2.0030
>> %flatten YR
>> time=YR(:);
>> %see what we have
>> whos
  Name         Size                   Bytes  Class

  WK          52x20                    8320  double array
  YR          52x20                    8320  double array
  ans          5x1                       40  double array
  arr2d        2x3                       48  double array
  m            1x1                        8  double array
  n            1x1                        8  double array
  omaha       52x20                    8320  double array
  omahav    1040x1                     8320  double array
  time      1040x1                     8320  double array
  weeks       52x1                      416  double array
  year         1x20                     160  double array

Grand total is 5285 elements using 42280 bytes

>> %plot omahav(time):
>> plot(time,omahav)
>> %save the data
>> save omaha time omahav
>> %clear & load
>> clear
>> load omaha
>> whos
  Name         Size                   Bytes  Class

  omahav    1040x1                     8320  double array
  time      1040x1                     8320  double array

Grand total is 2080 elements using 16640 bytes

>> diary off
>> edit tut1diary.txt
>>