如何将分析案例解决为数值案例？

Question

我想根据 p_h 和 p_s 求解 y 变量。从分析案例来看，它没有解决，给定参数的输出是空列表。我想把它作为一个数值案例来解决。我怎样才能做到这一点？我如何用数值解根据 p_h 和 p_s 求解 y？

定义函数（）：

α , h0 , π , τ , θ = 0.7 , 1 , 0.05 , 0.02 , 0.02
y , p_h , p_s = sympy.symbols('y p_h p_s')

c_r  = α*((y/(1+π))+τ)
s_r  = ((1-α)/p_s)*((y/(1+π))+τ)
c0_o = (y/(1+π))-((p_h*θ*h0**2)/(1+π))+τ

expression = (((c_r**α)*(s_r**(1-α))) - ((c0_o**α)*(h0**(1-α))))**2
y00 = solve(expression, y,  dict = True)

return y00

函数（）

我知道我不能通过最小化来解决它，因为这种方法没有给我一个符号解决方案。那么我可以用什么方法解决呢？你能举个例子吗？

注意：我无法更改参数的值。