function [t,y] = psbvp(f,N)
% Solve linear BVP y'' = f(t), y(-1)=0, y(1)=0
[D,t] = cheb(N);        
D2 = D^2;
D2 = D2(2:N,2:N);           
y = D2\f(t(2:N));
y = [0;y;0];              
