# Simplify to real and positive numbers

eng8 = [solve(eng==(k-1)*md*Vm1^2/(2*Patm*Vs*((Pc-(P2*s6))-(Pc-(P2*s6))^(1/k))),Vm1) for s6 in srange(r9,r10,s5,include_endpoint=1)]
v4=[v[0].rhs().n() for v in eng8]
v5=[v for v in v4 if v >= 0]


I can't seem to get v5 to output both real and positive numbers, only an empty list. All the inputs are positive, and I know the square root of the squared gives both positive and negative output. Any help?

edit retag close merge delete

Have you verified that the results in eng8 are indeed positive real numbers? The statement "the square root of the squared gives both positive and negative output" does not make mathematical sense. What you probably mean is that the equation $x^2 = a$ has two solutions, one positive and one negative when $a > 0$. These are the principle square root of a (usually represented by sqrt(a)) and the negative of the principle square root of a.