from physconst import *
from isotope import ion as I


from vegie import *
P, v, u = variables('$P--, v+, u+')
i,j,k,l,m,n = indices('i,j,k,l,m,n')