bordered system solution after govaerts solve (A b; c' d)(x y)=(f g)
0001 function [x,y] = bss(A,b,c,d,f,g) 0002 % bordered system solution after govaerts 0003 % solve (A b; c' d)(x y)=(f g) 0004 vs=A'\c; dels=d-vs'*b; 0005 v=A\b; del=d-c'*v; 0006 y1=(g-vs'*f)/dels; 0007 f1=f-b*y1; g1=g-d*y1; 0008 w=A\f1; y2=(g1-c'*w)/del; 0009 x=w-v*y2; y=y1+y2; 0010 end 0011