read "rga_variant.mpl";

t2:=time():
MainProc([u[x,x]^2*v+u[x,x]*v+u[x],u[x,y],u[y,y]-1],R):
t2:=time()-t2;