PROGRAM Pol1;
{Demonstruje se vypocet polynomu a jeho prvni derivace
Koeficienty jsou a[i] = i+1 }
CONST
  np = 100;
TYPE
  PolKoefNP = Array[0..np] OF REAL;
LABEL
  1;
VAR
  a: PolKoefNP;
  n,i:integer;
  x,dpol,pol:real;
PROCEDURE vyppol(x:real;n:integer;a:PolKoefNP;VAR p,dp:real);
BEGIN
  p := a[n];
  dp := 0;
  FOR i:=n-1 DOWNTO 0 DO BEGIN
      dp:= dp*x + p;
      p := p*x + a[i]
  END;
END;
BEGIN
  WRITELN(' ZADEJ RAD POLYNOMU');
  READLN(n);
  FOR i:=0 TO n DO
      a[i] := i+1;
  WHILE true DO BEGIN
      WRITE(' ZADEJ x (x<-1.e+30 konci vypocet) >>');
      READLN(x);
      IF (x < -1.e+30) THEN goto 1;
      vyppol(x,n,a,pol,dpol);
      writeln('X=',x,' POL=',pol,' DPOL=',dpol);
      END;
 1:
END.