function Abl = diffqt ( y , x )
#
# Funktion zur Berechnung der Ableitung von y(x)
# y(x) basiert auf nicht-equidistanten x-Werten
# x ist Vektor der x-Werten
#
maxAbl=0;
#Dummy=size(y);
#imax=int32(Dummy(2));
imax=length(y);
#
A=((y(3)-y(2))*(x(2)-x(1))-(y(2)-y(1))*(x(3)-x(2)))/ ...
  ((x(3)^2-x(2)^2)*(x(2)-x(1))-(x(2)^2-x(1)^2)*(x(3)-x(2)));
B=((y(3)-y(2))-A*(x(3)^2-x(2)^2))/(x(3)-x(2));
C=y(2)-A*x(2)^2-B*x(2);
Abl(1)=2*A*x(1)+B;
for i=2:imax-1
  A=((y(i+1)-y(i))*(x(i)-x(i-1))-(y(i)-y(i-1))*(x(i+1)-x(i)))/ ...
    ((x(i+1)^2-x(i)^2)*(x(i)-x(i-1))-(x(i)^2-x(i-1)^2)*(x(i+1)-x(i)));
  B=((y(i+1)-y(i))-A*(x(i+1)^2-x(i)^2))/(x(i+1)-x(i));
  C=y(i)-A*x(i)^2-B*x(i);
  Abl(i)=2*A*x(i)+B;
endfor
Abl(imax)=2*A*x(imax)+B;
#
endfunction
#