7, fix bug of q/(q+1)

the old code was (q+1)/q

Besides, code powerfamilyNETtestfunction.f90 prints variables in Fortran
to see what happened.