clear();


function dx=der(t,x)
  

 dx(1)=- k*(x)

 
endfunction 

k=1e1;

t=0:0.001:0.1;
x0=[1];

x=ode(x0,0,t,der);


plot2d(t,x);
DAT2(:,1)=t;
DAT2(:,2)=x'+rand(x','normal')/40 ;// creates  the 'experimental data'

plot2d(t,x, style =5);
plot2d(DAT2(:,1),DAT2(:,2), style = -1);

save('DAT2.dat' ,DAT2);

clear()

load('DAT2.dat', 'DAT2') ; //load the 'experimental data'

//now I would like to do the reverse procedure evaluating the k parameter