clear all;
close all;
%
% N = # of intervals,
% ax,bx = the endpoints of the interval [a,b]
%
N = 10;
ax = 0;
bx = 4;
%
% set grid points:  
gridchoice = 'chebyshev';          % see xgrid.m for other choices
x = xgrid(ax,bx,N,gridchoice);
%
plot(x,zeros(length(x)),'bo');
