function x=uniform(a,b,N);
x=linspace(a,b,N+2)';