subroutine sub1(f,f1,x)
common B1, B2, B3, C1, C2, c0, r21

f =phi(x)+C1

f1 =dphi(x)

END

subroutine sub2(f,f1,x)
common B1, B2, B3, C1, C2, c0, r21

f =phi(x)+C2

f1 =dphi(x)

END