function [x,h,y] = ini(g)
x = rand(1,g);
x = sort(x);
noise = rand(1,g)/5-0.1;
h = 0.5+0.4*sin(2*pi*x)+noise;
x11 = [0:0.02:0.98];
y = 0.5+0.4*sin(2*pi*x11);
plot(x,h,'ro',x11,y,'b');