% M-File grxxm025.m
% Slope field of  y' = f(x,y)

figure('Position',[50,100,800,800]);

% Funktionsdefinition f(x,y)
f=inline('x.^2+y.^2','x','y');
%f=inline('y','x','y');

% Rechteckgebiet, Rechteckgitter
a=2;
b=3;
hx=1/8;
hy=1/8;
[x,y]=meshgrid(-a:hx:a,-b:hy:b);
[m,n]=size(x);
dxx=ones(m,n);
dyy=f(x,y);

% Normierung der Vektoren (dxx,dyy)
for i=1:m
    for j=1:n
        u=dyy(i,j);
        dxx(i,j)=1/sqrt(1+u^2);
        dyy(i,j)=u/sqrt(1+u^2);
    end;
end;

% Plot 
hold on all;
quiver(x,y,dxx,dyy,'b');
axis([-a a -b b]);
plot([-a;a],[0;0],'k',[0,0],[-b,b],'k');

% Titel
ff=formula(f);
T=['Richtungsfeld von dy/dx =  ' ff];
title(T,'FontSize',20);
