function [] = prog1main() 
n=3;  !N dim.

for i= 1:n
  a=0;
   for j = 1:n 
       x=1/(i+j-1);
       y=str2num(sprintf('%.4f\n',x));
       A(i,j) =y;
      a=a+x;
   end 
   B(i,1)=str2num(sprintf('%.4f\n',a));
end
! Hilbert Matrix zuweisen

disp(A)
disp(B)
C=A^(-1)*B;
disp(C)
! benutzt die Funktion aus matlap 
for j=1:n
    g=abs(A(j,j));
     for i=j:n
         if abs(A(i,j))>g
             g=A(i,j);
             x=i;
         end
     end
     
     if A(j,j)==g
     else    
        t=A(j,j:n);
        A(j,j:n)=A(x,j:n);
        A(x,j:n)=t;
        t=B(j,1);
        B(j,1)=B(x,1);
        B(x,1)=t;
     end
    ! Spalten-Pivot-Suche
     for i=j+1:n
      L=A(i,j)/A(j,j);
      A(i,j:n)=A(i,j:n)-A(j,j:n)*L;
      B(i,1)=B(i,1)-B(j,1)*L;
     end 
    ! Gauss Elimination
 end 
  X(1:n,1)=0;
   for i=1:n
      BB=A(n+1-i,1:n)*X(1:n,1);
      X(n+1-i,1)=(B(n+1-i,1)-BB)/A(n+1-i,n+1-i);
   end
  
  disp(X)
  disp(A)
  disp(B)
  ! die Ergebniss mit obige zu vergleichen
return