Вы можете сделать что-то вроде следующего, используя два логических значения вместо реляционного 0 <= x <= p
.
import sympy as sp
x, h, p = sp.symbols('x h p')
fx = sp.Piecewise(
((x*h)/p, ((0 <= x) & (x <= p))),
((h*(sp.pi - x))/(sp.pi - p), ((p <= x) & (x <= sp.pi)))
)
, что дает