% n inervals means n-1 unknonws
% b=x*(1-x)*h^2 
function [A, b]=get2ptbdvp(n)
A=diag(2*ones(n-1,1), 0)-diag(ones(n-2,1), 1)-diag(ones(n-2,1), -1);
h=1/n;
x=[0:h:1];
x=x(2:n);
b=h^2*x.*(1-x); 
b=b';