include 'horner.f'
c     DOLACZENIE PODPROGRAMU 'horner.f"      
      PROGRAM HOR4
c     "PROGRAM GLOWNY"      
c     OBLICZA WARTOSC WIELOMIANU I JEGO POCHODNEJ 
c     DLA x Z PRZEDZIALU [aa,bb]
c     DEKLARACJE
c     W PROGRAMIE GLOWNYM DEKLARUJEMY MAKSYMALNE WYMIARY
c     SA ONE ZAWSZE STALE
      implicit none
      integer N,j,k,pyt
      real a,b,c,x,P,Q,h,aa,bb
      dimension a(0:100),b(0:100),c(0:100)
   10 write(*,*)'PROGEAM HOR4'
c     PRZYGOTOWANIE DO PISANIA NA DYSK      
      open(unit=10,file='WYNIK')
c     CZYTANIE DANYCH
      write(*,*)'PODAJ N - STOPIEN, 0<N<101'
      read (*,*) N
      write(*,*)'PODAWAJ a(j), j=0,1,...,N'
      do 1,j=0,N,1
      read(*,*) a(j)
    1 continue
      write(*,*)'PODAJ PRZEDZIAL [aa,bb]: aa<x<bb'
      read (*,*)aa,bb
c     OBLICZANIE  SKOKU DLA x      
      h=(bb-aa)/100
c     OBLICZANIE WARTOSCI WIELOMIANU I POCHODNEJ       
      do 2,j=0,100,1
      x=aa+j*h
c     WOLANIE PODPROGRAMU OBLICZAJACEGO WARTOSC WIELOMIANU (P)      
      call horner(N,a,b,x,P)
c     WOLANIE PODPROGRAMU OBLICZAJACEGO WARTOSCI POCHODNEJ (Q)     
      call horner(N-1,b,c,x,Q)
c     PISANIE NA DYSK DO PLIKU "WYNIK"      
      write(10,*) P,Q
    2 continue 
c     ZAMKNIECIE PLIKU "WYNIK"
      endfile(10)
      close(10)
c     ZAKONCZENIE      
      write(*,*)'GOTOWE'
      write(*,*)'<ENTER>'
      read (*,*)
      write(*,*)' 0 - KONIEC'
      write(*,*)' 1 - OD POCZATKU'
      read (*,*) pyt
      if(pyt.eq.0)goto 100
      goto 10
  100 stop 100
      end