(*ALGORITHMUS 8.5*)
(*$B-*)
PROGRAM fktplt;
VAR i, xg, yg, xt, yt, oldx, oldy,  farbe, keine, fall: integer;
    xmin, xmax, ymin, ymax, x, h : real ;

(*$I plotproz *)

  FUNCTION f(x:real):real;
  CONST eps = 0.1;
  VAR a,b,y,ya,yn : real;
  BEGIN
   CASE fall OF
     1: BEGIN
          a:=-5; b:=5;
          y:=(a+b)/2;
          WHILE b-a > eps DO
          BEGIN
            IF y+exp(y)<x THEN a:=y ELSE b:=y;
            y:=(a+b)/2;
          END;
          f:=y;
        END;
     2: BEGIN
          yn:=0; (* Startwert  *)
          REPEAT
            ya:=yn;
            yn := ya-(ya+exp(ya)-x)/(1+exp(ya));
          UNTIL abs(ya-yn)<eps;
          f:=yn;
        END;
  END;
END;
BEGIN
   writeln('Fall eingeben: (1) = Bisektion ');
   writeln('               (2) = Newton ');
   readln(fall);
   writeln('xmin, xmax, ymin, ymax eingeben');
   readln(xmin, xmax, ymin, ymax);
   plotbegin( xg, yg, xt, yt, farbe, keine);
   achsen(xmin, xmax, ymin, ymax, xg, yg, xt, yt, farbe, keine);
   pl(xmin, f(xmin), keine);
   h := (xmax-xmin)/xg;
   for i := 0 to xg do
   begin
      x :=xmin+ i*h;
      pl(x, f(x), farbe)
   end;
END.