fortran ループでの二項係数

Question

loopでの二項係数

fortranでloopで正の整数ｎを読み込みr=1,2,...nに対するnCr=n!/r!*(n-r)!　二項係数を計算して表示するプログラム
を作っていますが、階乗を利用するのかと思って階乗のプログラムをうちこみましたが、反応していないと思います。どこを直したらよいでしょうか。
n= 4
nCr= 1
nCr= 0
nCr= 0
nCr= 0

試したこと

program ensyu5
implicit none
integer::n,m,r,s,Y,i

if (n>0) then
end if

write(,) 'n='
read(,) n
Y=1
do r=1,n

do m=1,r
n=n*m
end do

Y=n/(r*(n-r))

write(,)'nCr=', Y
end do
stop
end program ensyu5

Comments

[coco_bauer]

どこに　nの階乗を計算するコードがありますか？ nKaijo=1 // do x=2,n // nKaijo=nKaijo*x // end do (//は改行) のような感じになるかと思うのですが。

[pomme_gramimg]

do r=1,n// do m=1,r// n=n*m //end do　このあたりだと思っていますが後にrがあるのでrを使っています。