function [list]=liste(mat)
% Variable "integral": Integral des Poynoms aus Aufgabenteil b)
% Variable "list": Liste mit den Koeffizienten der orthogonalen Polynome
for i=1:min(size(mat)-[0 1]);
if abs(mat(i,i))<1E-13
warning('Pivotelement (%i,%i) entspricht Null!',i,i);
zeila=abs(mat(i+1:end,i)); zeila=i+find(zeila==max(zeila));
mat([i,zeila],:)=mat([zeila,i],:); %interchange rows
if abs(mat(i,i))<1E-13
warning(['Kein Pivot-Element gefunden!']);
end
end
if abs(mat(i,i))>=1E-13
zeila=[i+1:size(mat,1)];
mat(zeila,:)= mat(zeila,:)-(mat(zeila,i)./mat(i,i))*mat(i,:);
end
end
mat(abs(mat)<1E-13*max(max(abs(mat))))=0;
x=zeros(1,size(mat,2)-1);
for i=(size(mat,1):-1:1)
if i==size(mat,1);
x(i)=mat(i,end);
else x(i)=mat(i,end)-mat(i,[i+1:end-1])*x(i+1:end)';
end
if 1E-13>abs(x(i));
x(i)= 0; % frei wählbar
else x(i)=x(i)/mat(i,i);
end
end
x(isnan(x))=0; % 0 * x(i) = 0: frei wählbar
if any(x==Inf | x==-Inf); warning('Gleichungssystem nicht lösbar!');
list=mat; x=[];
else list=mat(:,1:size(mat,2)-1);
end;
% Ausgabe:
list
x