function x_1=model_simple(x_0,system) eps=0.01; x_1=sin(x_0)./(x_0.^2+eps)+system.sigma*rand(size(x_0)); end