% Harj. 4 AV 10-12.10.2001   HA
%
% Teht. 1
%
format rational;format compact
A=[4 0 1;0 4 1;1 1 -1]
b=[16;32;0]
rref([A b])
A\b
ans =

     1     0     0     2
     0     1     0     6
     0     0     1     8


ans =

     2
     6
     8
%
% Teht. 2
%
A=[1 -1 1 0 0 2;1 1 2 4 1 2;0 2 1 4 3 1;1 -3 0 -4 -2 1]
b= [0; 3; 2; 0]
Ab=[A b]

[m,n]=size(A);
% (1,1)
pivot=Ab(1,1)
for i=2:m
   Ab(i,:)=Ab(i,:)-(Ab(i,1)/pivot)*Ab(1,:)
   pause
end
%(2,2)
pivot=Ab(2,2)
for i=3:m
   Ab(i,:)=Ab(i,:)-(Ab(i,2)/pivot)*Ab(2,:)
   display('Huitaise ENTERiä, kun haluat seuraavan nollan')
   pause
end
% rivi 3
% 1. tukisarake on sarake 5.
pivot=Ab(3,5)
for i=4:m
   Ab(i,:)=Ab(i,:)-(Ab(i,5)/pivot)*Ab(3,:)
   display('Huitaise ENTERiä, kun haluat seuraavan nollan')
   pause
end
break
Nyt on kaikki valmista, yksi 0-rivi, johtavat sarakkeet

Ab =
  Columns 1 through 6 
      1           -1            1            0            0            2      
      0            2            1            4            1            0      
      0            0            0            0            2            1      
      0            0            0            0            0          -1/2     
  Column 7 
      0      
      3      
     -1      
     5/2 
     Nyt on kaikki valmista, ei 0-rivejä, 
     johtavat sarakkeet: 1,2,5,6, r=4.
     Ratk:
     x_1
rref(Ab)
ans =
  Columns 1 through 6 
      1            0           3/2           2            0            0      
      0            1           1/2           2            0            0      
      0            0            0            0            1            0      
      0            0            0            0            0            1      
  Column 7 
    21/2     
     1/2     
      2      
      -5
