[Octave-bug-tracker] [bug #41192] Basic matrix inversion fails

[anonymous]

URL:
  <http://savannah.gnu.org/bugs/?41192>

                 Summary: Basic matrix inversion fails
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Fri 10 Jan 2014 07:26:51 PM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Crash
                  Status: None
             Assigned to: None
         Originator Name: Saul Mtakula
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 3.6.4
        Operating System: Microsoft Windows

    _______________________________________________________

Details:

The following code fails in the last statement; the matrix inversion crashes
Octave. Changing the code to inv() also results in a crash. Matlab or Scilab
execute code without a crash.


clear
clc

A = [ 0 0 1.1320 0 -1 
      0 -0.0538 -0.1712 0 0.0705 
      0 0 0 1 0 
      0 0.0485 0 -0.8556 -1.013 
      0 -0.2909 0 1.0532 -0.6859 ];
B = [ 0 0 0 
      -0.120 1 0 
      0 0 0 
      4.4190 0 -1.665 
      1.5750 0 -0.0732 ] ;
      
C = [ eye(3,3) zeros(3,2) ] ;

D = zeros(3,3) ;

Gam=B;

      Aw=-0.001*eye(3,3);
      Bw=eye(3,3);
      Cw=eye(3,3);
      Dw=zeros(3,3);
      
      Aa=[A B*Cw; zeros(size(Aw,1), size(A,2)) Aw];
      
      
      Ga=[B*Dw;Bw];
      
     
      
      j=i;
      Z=(0.001*j*eye(size(Aa))-Aa);
      Z1=Z\Ga;
     




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?41192>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/