clear all
rng(1122334457);
%Replace the above line with the one below for Octave
%randn("seed",1122334457);
global a b c
a=[1 1];
b=[0 0];
c=[1 2];
x=[0:0.1:4*pi];
model=@(t) f(1,f(2,0,t),t);
eta=0.1*randn(size(x));
y=model(x)+eta;
plot(x,y,".");
