function y = noisysin(x)

y = sin(x) + randn(size(x))*0.001;
end
