Path: news.mathworks.com!not-for-mail
From: "Bruno Luong" <brunoluong@yahoo.com>
Newsgroups: comp.soft-sys.matlab
Subject: Re: Matrix Equation solving
Date: Mon, 19 Nov 2007 21:27:19 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 10
Message-ID: <fhsv3n$g3p$1@fred.mathworks.com>
References: <fhe92r$m23$1@fred.mathworks.com> <fhna35$9u4$1@fred.mathworks.com> <fhsnu5$hg7$1@fred.mathworks.com>
Reply-To: "Bruno Luong" <brunoluong@yahoo.com>
NNTP-Posting-Host: webapp-03-blr.mathworks.com
Content-Type: text/plain; charset="ISO-8859-1"
Content-Transfer-Encoding: 8bit
X-Trace: fred.mathworks.com 1195507639 16505 172.30.248.38 (19 Nov 2007 21:27:19 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Mon, 19 Nov 2007 21:27:19 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 1184112
Xref: news.mathworks.com comp.soft-sys.matlab:438416


Alternative solution:

D=-eye(3);
Z=zeros(3);
BIGA1=[A Z Z; Z A Z; Z Z A]';
BIGA2=reshape(permute(reshape(BIGA1,3,3,3,3),[2 1 4 3]),9,9);
P=reshape((BIGA1*BIGA2-eye(9))\D(:),3,3)

;-)
Bruno