If I understand correctly, the G function in f3 is the hypergeometric function 2F1. There appears to be a package to calculate this in R -- look for "hypergeo" with function: ...