Below is a small part of what I’m trying to run, it gets stuck on the final line

Below is a small part of what I’m trying to run, it gets stuck on the final line

I’m trying to create a custom symbol in sympy that behaves like r = √(x² + y² + z²) — specifically, ∂r/∂x = x/r:

Consider the following sympy expression

After a page of sympy manipulations I get and expression involving LambertW, but the answer given is for the zeroth branch cut. I can tell that I wanted the -1 branch. So, I’d like to change every occurrence of LambertW(y) into LambertW(y, -1).real, where y is an arbitrary sympy expression (that is different in each instance), in a set of expressions that are large enough so that doing it by hand is error prone.

I am using the following code, which returns the integral itself and not its result:

According to the documentation one can make substitutions as:

How does sympy‘s nsolve do what it does?

I’m trying to display a sympy matrix in scientific notation, and it is only working for values less than 1. I can’t understand why applyfunc is behaving differently. Is there a better way to do this?

I’m trying to display a sympy matrix in scientific notation, and it is only working for values less than 1. I can’t understand why applyfunc is behaving differently. Is there a better way to do this?