oscillator_experiments.m

function res = oscillator_experiments(step)
%function res = oscillator_experiments(step)
%
% experiments with Steffen's oscillator model
%
% step 1: several detailed simulations and plot

% B.Haasdonk and S. Waldherr 5.8.2010
% Time-stamp: <Last change 2010-09-16 14:52:19 by Steffen Waldherr>



if nargin < 1
  step = 1;
end;


res = [];

switch step
 case 1 
  
  model = oscillator_model;
  model_data = gen_model_data(model);
  model = model.set_mu(model,[0.1]);
  sim_data  = detailed_simulation(model,model_data);
  plot_sim_data(model,model_data,sim_data,[]);
  
    model = model.set_mu(model,[0.2]);
  sim_data  = detailed_simulation(model,model_data);
  plot_sim_data(model,model_data,sim_data,[]);

      model = model.set_mu(model,[0.3]);
  sim_data  = detailed_simulation(model,model_data);
  plot_sim_data(model,model_data,sim_data,[]);
  keyboard

  case 2 
  
%  A = sparse([],[],[],numstates,numstates,size(S,2)*numstates);
%  for j=1:size(S,2)
%    A = A + p{j}*Acoeff{j};
% end;
  
%  [t,P] = ode15s(@(t,x) A*x,[0 1], Pic);
 
 otherwise
  error('step unknown')

end;