function Richtungsfeld
% dgl ist die erste Ableitung von y nach x, entspricht also
% dem Ausdruck y'(x), bzw y'(y,x)

y = -4:1:4;
x = -4:1:4;
 
 
for yn = 1:length(y)
  for xn = 1:length(x)
    len = sqrt( dgl(y(yn), x(xn))^2 + 1 );  % Länge des Vektors für Normierung
    dx(yn,xn) = 1 / len;                    % Länge des Vektors entlang der Abszisse
    dy(yn,xn) = dgl(y(yn), x(xn)) / len;    % Länge des Vektors entlang der Ordinate
  end
end
 
quiver(x, y, dx, dy, '-r');  % Vektoren zeichnen
grid on

%print('field.png', '-dpng')
end


function y=dgl(x,t)
	
y = (t+x)./(t-x);

end


