x=u_x(:,2)
%hist(x)
t = linspace (0,1,11)
xr=zeros(10,1)
for k=1:10
  sum=0
  for i=1:300
    if(x(i)>=t(k) & x(i)<t(k+1))
      sum=sum+1
    end
    xr(k)=sum/300
  end
end
tc=t(1:10)+0.05
%plot(tc,xr)
pb = xr

GHelp(2)
min = fminsearch(@(alpha) G, [0], pb, t)

