read "rga_variant.mpl"; t1:=time(): E1:=[u[x]^2-4*u,u[x,y]*v[y]- u + 1,-u[x]+v[x,x]]: MainProc(E1,R): t1:=time()-t1;