 n=input('n=');

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

u=col_bih1d(n,x,'fbih','p','q',0,-pi,0,-pi);

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