matrix close to singular

[amit Kumar nigam]

hi all,
i am working on a problem which involves the generalized eigenvalue
problem:

        given any two square matrices A and B, find matrices X and
        D that satisfy the following equation:

                A*X = B*X*D

        where D is a diagonal matrix of generalized eigenvalues and
        X is a matrix whose columns are the generalized
        eigenvectors.

on the octave prompt if i type the following

        [X,D] = geig(A,B)

where geig.m is a user-defined function (below) that uses the qz
function

        function [V,D] = geig(A,B)
          [AA, BB, Q, Z, V, W, lambda] = qz (A, B);
          D = diag(lambda);
        endfunction

this should be equivalent to the matlab `eig' function.

for certain matrices i get the following warning:

 warning: inverse: matrix singular to machine precision, rcond =
 5.33665e-45

in matlab i get a similar message:

 matrix is close to singular or badly scaled. results may be
 inaccurate. rcond = 3.45674e-56

is there a way of going round the problem of getting inaccurate
results. if there is a way of rescaling the matrix appropriately
please let me know?

any help will be greatly appreciated.


---
Amit K. Nigam
Department of Mechanical Engineering, Concordia university,
Montreal 
Quebec Canada
Ring:1-514-8484199 (O), 1-514-9346381(R)
Great minds have purpose....Little minds have wishes..



Get 250 color business cards for FREE!
http://businesscards.lycos.com/vp/fastpath/



-------------------------------------------------------------
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
-------------------------------------------------------------