disp('Running MATLAB script file SiCur1fig3.m') ;
%
%    Makes Figure 3 for the SIgnificance of CURvature paper
%
%    Copied from Old\SiConToyEG1.m
%
%    This analyzes the Toy SiCon Data



pstr = 'SiCur1fig3' ;



figh = figure(1) ;
clf ;


%  Generate data
%
left = 0 ;
right = 1 ;
bottom = 0 ;
top = 1 ;

ndat = 200 ;
sig = 2 ;

seed = 32876429 ;
randn('seed',seed) ;

datx = linspace(left,right,ndat)' ;
daty = sin(2 * pi * 4 * datx) ;
daty = daty + 2 + 20 * datx ;
daty = daty + sig * rand(ndat,1) ;

data = [datx,daty] ;





%  do SiZer analysis
%
paramstruct = struct('iout',2,...
                     'imovie',0,...
                     'hhighlight',0,...
                     'ibigdot',1,...
                     'famoltitle',[],...
                     'sizertitle',[],...
                     'curvsizertitle',[],...
                     'xlabelstr',['x'],...
                     'ylabelstr',['y'],...
                     'labelfontsize',12,...
                     'ibdryadj',0,...
                     'iscreenwrite',1) ;

sizerSM(data,paramstruct) ;





  orient portrait ;
  set(gcf,'PaperSize',[5.0, 9.0]) ; 
  set(gcf,'PaperPosition',[0.25, 0.25, 4.5, 8.5]) ; 
print('-dpsc', [pstr '.ps']) ;


