function of Eq. (2)
syntaxSUBROUTINE DAGRB (NY,NV,NU,NP,NRB,T0,Y0,V0,U0, TF,YF,VF,UF,P,RB)
parametersIMPLICIT NONE