 n=input('n=');

for i=1:n+1
	x(i)=-cos((i-1)*pi/n);
end

u=col_bvp(n,x,'f','p','q',1,0,0,1,0,0);

for i=1:n+1
u(i)=u(i)-sin(pi*x(i));
end
err=max(abs(u))
