function b=seqfit(sequenzliste,verschiebungen,wichtung);
%
% function b=seqfit(sequenzliste,verschiebungen,wichtung)
%
% Funktion zur Berechnung der Abstandsmatrix zwischen
% den Verschiebungen identifizierter Reste in den 
% NMR-Spektren "verschiebungen" und den aufgrund
% einer Kristallstruktur vorhergesagten Verschiebungen.
%
% In "sequenzliste" und "verschiebungen" wird jede Zeile
% als ein Rest im NMR-Spektrum interpretiert. Stellen,
% an denen "0" steht, werden ignoriert.
%
% "wichtung ist ein optionaler Parameter der in einem
% Zeilenvektor die Wichung zwischen den Dimensionen
% festlegt.
%
% Eine Spalte entspricht jeweils einem Rest in der Sequenz.
%
%
% Es folgt die extra ausführliche Anleitung auf ausdrücklichen 
% Wunsch eines einzelnen Herrn im Wasserwerk:
%
% Einladen von "Sequenzliste" mit
% <sequenzliste_variablein_name> = loadascii ('<sequenzliste_filename>');
% 
verschanz=size(sequenzliste,2);
sequenzliste=sequenzliste(:,4:verschanz);
verschanz=verschanz-3;
seqnmranz=size(verschiebungen,1);
if nargin==2
 wichtung=ones(1,verschanz);
end
seqlaenge=size(sequenzliste,1);
ignorieren=sum(verschiebungen(:)==0);
for j=1:seqnmranz
 h=(ones(seqlaenge,1)*(wichtung.*verschiebungen(j,:))-(ones(seqlaenge,1)*wichtung).*sequenzliste).^2;
 h(:,verschiebungen(j,:)==0)=0;
 b0(:,j)=sum(h,2);
end
h=prod(size(b0));
for j=1:seqlaenge-seqnmranz+1
 b(j)=sum(b0(j:seqlaenge+1:h));
end
if prod(size(verschiebungen))-ignorieren~=0
 b=sqrt(b./(prod(size(verschiebungen))-ignorieren));
end