xbinavg(x, y, xedge, medopt)

Simple function for summarizing a large number of 2d-data.
Calculating mean and standard deviation of y-datas within each bin of x-data.
Example code:
--
x = rand(1,1280)*100;
y = x.^2 + 5000*randn(size(x));
[meanv, stdv, xv, nv] = xbinavg(x, y, 0:10:100, true);
scatter(x,y, [], [.8 .8 .8], '.');
hold on;
errorbar(xv,meanv,stdv, 'k','LineWidth', 5/2)
---