function coeff=mypolyfit(x,z,n)
% selbstgestricktes polyfit, ohne polyval
m=length(x);
xx=reshape(x,m,1);    % Spaltenform sichern
zz=reshape(z,m,1);    % Spaltenform sichern
valmat=zeros(m,n+1); % wird Matrix der Polynomwerte 
valmat(:,n+1)=ones(m,1); % hinten Einsen 
for i=n:-1:1         % rekursive Potenzen
    valmat(:,i)=xx.*valmat(:,i+1);
end
coeff=(valmat\zz)';