In [ ]:
from umatgen import *

var('parA, parB, parC, parD, parNu')

MR=umat(Implicit=True, project='Test_MooneyRivlin',
        outfile="MooneyRivlin", outroutine="MR")

parC=1/2*parA+parB
parD=(parA*(5*parNu-2)+parB*(11*parNu-5))/(2*(1-2*parNu))

C=MR.C

MR.W=(parA*(I1(C)-3)+parB*(I2(C)-3)+parC*(I3(C)^(-2)-1)+parD*(I3(C)-1)^2)

MR.write()

In [ ]: