フォートランに詳しい方、プログラムの説明をお願いしますm(__)m
open(11,file='dmvd.d')
dx=(xf-xi)/dfloat(nn)
do 10 i=0,nn
xx=xi+dx*dfloat(i)
ptx=fac*xx*xx*dexp(-xx*xx)
c calculate Ptilde(x)
write(11,21)xx,ptx
21 format(d15.7,2x,d15.7)
10continue
close(11)