help-octave
[Top][All Lists]

## qz

 From: Stefan van der Walt Subject: qz Date: Thu, 15 Jul 2004 21:02:51 +0200 User-agent: Mutt/1.5.6+20040523i

```I am trying to solve the generalised eigenvalue problem using qz, but
the following behaviour looks suspicious:

octave:5> a = rand(3,3); b = rand(3,3); [AA, BB, Q, Z, V, W] =
qz(a,b); Q*a*Z - AA, Q*b*Z - BB
ans =

-0.506316   0.133307   1.509299
-0.425943  -0.168888  -1.856508
0.203352  -0.028521  -0.327870

ans =

-0.377653  -0.108805   1.424638
0.274852  -0.178300  -0.189498
0.098925   0.051424  -0.448613

If I understand correctly, the answers should be much closer to
zero(3,3)?

Regards
Stefan

-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------

```

reply via email to