1;

function richtungsfeld(dgl)
% dgl ist die erste Ableitung von y nach x und ist i.A. eine Funktion von x und y
% Ausschnitt und Abstand zwischen den Vektoren
y = -5:.5:5; x = -5:.5:5;
for y_n = 1:length(y)
 for x_n = 1:length(x)
   len = sqrt( dgl(y(y_n), x(x_n))^2 + 1 ); % Länge des Vektors für Normierung
   dx(y_n,x_n) = 1 / len;                   % Länge des Vektors entlang der Abszisse
   dy(y_n,x_n) = dgl(y(y_n), x(x_n)) / len; % Länge des Vektors entlang der Ordinate
 end
end 
h=quiver(x, y, dx, dy,0.5,"r","linewidth",1); % Vektoren zeichnen 
set (h, "maxheadsize", 0.1); 
xlabel ("x");
ylabel("y");
print('field.svg', '-dsvg')  % Plot als svg-Datei exportieren
% Ende des Files
endfunction

%source("richtungsfeld.m")
dgl = @(y, x) y-x   % Funktionsdefinition 
richtungsfeld(dgl)
