import sympy as s
H0,z,g=s.symbols('H0 z g')
def f(H0,g,z):
H1=(6/H0)/(g+(6-g)*(1+z)**3)
return H1import sympy as s
H0,z,g=s.symbols(import sympy as s
H0,z,g=s.symbols('H0 z g')
def f(H0,g,z):
H1=(6/H0)/(g+(6-g)*(1+z)**3)
return H1import sympy as s
H0,z,g=s.symbols(