(* Richtungsfeld fuer ein Differentialgleichung  *)

(* Definition von y':  y' = y+1-x  *)
yp[x_,y_]:=y+1-x;

xlow=-4; xup=4; xdel=0.5;
ylow=-4; yup=4; ydel=0.5;  

pair[x_,y_]:=(der=yp[x,y]; eps=0.1/Sqrt[der^2+1];
              {N[{x-eps,y-eps*der}],
               N[{x+eps,y+eps*der}]});

pts=Table[Point[N[{x,y}]],
               {x,xlow,xup,xdel},{y,ylow,yup,ydel}];

vectors=Table[Line[pair[x,y]],
               {x,xlow,xup,xdel},{y,ylow,yup,ydel}];

Show[Graphics[{vectors, pts}],
              Axes -> Automatic,
              AspectRatio -> Automatic]

(* Schneller geht es mit der eingebauten 
  Routine PlotVectorField: *)

(* Version f. Mathematice 5.x :*)
<<Graphics`PlotField`
fvec[x_,y_] := ({1,yp[x,y]}/Sqrt[1+yp[x,y]^2]);
PlotVectorField[fvec[x,y],{x,xlow,xup},{y,ylow,yup}]

(* Version f. Mathematice 6.x :*)
Needs["VectorFieldPlots`"];
fvec[x_,y_] := ({1,yp[x,y]}/Sqrt[1+yp[x,y]^2]);
VectorFieldPlot[fvec[x, y], {x, xlow, xup}, {y, ylow, yup}]
