function Z=reciprocalmethod(A, B, X, C)
%Z=reciprocalmethod(A, B, X, C)
%function that performs the reciprocal method
n=size(A,1);
Y=inv(eye(n)-A')*X;
Z=C+B'*Y;