load scrippsSST.txt
t=scrippsSST(:,3);temp=scrippsSST(:,5);
[a,b,f,t1]=myfft(t,temp);
edit myfft FourierMat
t2=linspace(1994, 1997,100);
x=FourierMat(a,b,f,t2,t1);
plot(t,temp,'b');hold on
plot(t2,x,'r')
clf
plot(f,sqrt(a.*a+b.*b))
a(1)
ans =
   17.9371
mean(temp)
ans =
   17.9371
A=[3 3 -999 4 4]'
A =
     3
     3
  -999
     4
     4
I=find(A~=-999)
I =
     1
     2
     4
     5
mean(A(I))
ans =
    3.5000
cellar=cell(2,3)
cellar = 
     []     []     []
     []     []     []
cellar{1,2}='abc'
cellar = 
     []    'abc'     []
     []       []     []
cellar{2,3}=2*3
cellar = 
     []    'abc'     []
     []       []    [6]
cellar{1,1}=ones(3,4)
cellar = 
    [3x4 double]    'abc'     []
              []       []    [6]
diary off