10 common/pyrvpm/rm(0:3),a(2),b(2),resm(2),resw(2),mflag
12 DOUBLE PRECISION x,
y,
pyrvg4, rm, a, b, resm, resw, rvs,
pyrvs
14 SAVE /pyrvpm/,/pyg2dx/
16 rvs=
pyrvs(
x,
y,resm(1),resw(1),resm(2),resw(2))
20 pyrvg4 = rvs*(rm(1)*rm(3)*a(1)*a(2)*(
x+
y-rm(1)**2-rm(3)**2)
21 & + rm(1)*rm(0)*b(1)*a(2)*(
y-rm(2)**2-rm(3)**2)
22 & + rm(3)*rm(0)*a(1)*b(2)*(
x-rm(1)**2-rm(2)**2)
23 & + b(1)*b(2)*(
x*
y-(rm(1)*rm(3))**2-(rm(0)*rm(2))**2))