bordered elim +1 after govaerts'90
0001 function [x,y]=belpo(A,b,c,d,f,g) 0002 % bordered elim +1 after govaerts'90 0003 [x0,y0]=bel(A,b,c,d,f,g); % block-elim 0004 fs=f-(A*x0+b*y0); gs=g-(c*x0+d*y0); % resi 0005 [x1,y1]=bel(A,b,c,d,fs,gs); % corrector 0006 x=x0+x1;y=y0+y1;