function x=test_cg_2dbvp(n)
[A,b]=get2Dbdvp(n);
tic;
[x,flag,relres,iter]=pcg(A,b,10^(-2),2*10^3);
toc
flag
iter

