A method is presented for solving a succession of complex matrix equations in which the phase of the real and imaginary components changes. The method is more efficient than the technique obtained by using complex Gaussian elimination on each of the matrix equations separately. In addition, some interesting theoretical relationships are presented for the solution of complex matrix equations in general, using only real-valued arithmetic operations.Key Words. Linear systems, unconstrained minimization, mathematical programming, finite-difference-finite-element methods, pseudoinverse solutions.t T h i s w o r k w a s p e r f o r m e d while the a u t h o r w a s with the Electrical E n g i n e e r i n g D e p a r t m e n t ,McGitl University, M o n t r e a l , C a n a d a . F i n a n c i a l s u p p o r t w a s p r o v i d e d b y the N a t i o n a l