program newton_symbol
implicit none

integer :: n ,k
integer :: silnia, licznik, mianownik


write(*,*) 'Wprowadz liczby z ktorych mam policzyc symbol Newtona N(n,k)'
read(*,*) n, k

if (n.gt.k) then

   licznik = silnia(n)
   mianownik = silnia(k)*silnia(n-k)

   write(*,*)'Symbol Newtona jest rowny', licznik/mianownik

else

   write(*,*)'n musi byc wieksze od k!'
end if


end

integer function silnia(a)
integer :: a
integer :: i

if (a==0) then
   silnia = 1.0
else 
   silnia = 1.0
   do i = 1, a
      silnia = silnia* i
   end do
end if

end function