% Plot the monomial basis polynomials 1, x, x^2, ...
% for [a,b] = [0,1] with x_j = j/5. 

 figure(1), clf
 xx = linspace(0,1,300);

 N = 5;
 x = linspace(0,1,N+1);
 plot(x,zeros(size(x)), 'k.','markersize',20)
 axis([0 1 -0.1 1.1])
 set(gca,'fontsize',16)
 title('monomial basis functions')
 hold on
 
 for j=0:N
    plot(xx, xx.^j, 'k-','linewidth',2), hold on
    input('press <return>')
 end


% axis equal
 axis([0 1 -0.1 1.1])
 set(gca,'fontsize',16)
 print -deps2 mono_basis