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