sympy | A computer algebra system written in pure Python

Good morning,

I am creating an animation by connecting 6 nodes coordinates at 10 different time steps.

I start with a test: I first create a static plot for a time equal to 2 (for example) and I get the exact static plot that I am expecting: all and only the consecutive nodes are connected.

Then I create the animation. Unfortunately, the animated plot connects the nodes in the wrong way. You can see that the animation connects the consecutive nodes, but also the second and the second to last nodes.

Any idea why? Thanks

Sympy lambdify can be used to bridge functionality between numpy and sympy. However, I could not find what shortcomings exist when using lambdify. For example, I am interested in using Max with sympy and numpy:

If sigma[i, j] == 0 or sigma[i, j] == 1 it doesnt matter. It only takes the position of non-zero elements. There is a code that provides finding coeffs y1, y2, w1, w2, which are sympy.symbols, of func f by sympy.solve. C - is a matrix.

The process of creating C-matrix:

In Python, I do some SymPy calculations that yield an array full of polynomials such as:

a*(a*(a*(a + b) + b*(a + b)) + b*(a*(a + b) + b*(a + b))) + b*(a*(a*(a + b) + b*(a + b)) + b*(a*(a + b) + b*(a + b)))

Note that this example happens to simplify to (a+b)**4, but this won't always be the case obviously. So how do I convert this expression to the following form:

c_1*a**4 + c_2*a**3*b + ... + c_n*b**4

And once I have such an expression, how would I extract the exponents c_1, ..., c_n? All I have is the .exp command, but it only works on expressions of the form a**n (i.e. no mixture of a and b and a coefficient of 1).

Any help would be majorly appreciated.

The following code does exactly what I want; however, the for loop is far too slow. On my machine, the wall time for the for loop is 1min 5s. I'm looking for an alternative to the for loop that is much faster.

the code below is from the book "Python for Probability, Statistics, and Machine Learning. The clarification needed is for the plotting section. The problem is that "logJ" in the script is not-defined. However, the book provides this as the code to plot the graph. How do you correct (code) the plotting part of the script so it plots the output shown?

I am doing some operations on two matrices in sympy and I want to record how the result was obtained. For example in a isympy session:

I am using lcapy together with sympy and trying to process complex numbers from a circuit.

I have the following sympy expression:

I am getting this sympy related comparison error. Do you know how I can fix it? It doesn't seem that authors in the original repo are facing it. I see that for one of the pkl files, namely 001.pkl, the y2d is imaginary (not sure why), and is equal to 1.05147176800692e-6 + 3.72470439655723*I.

I am running this code: https://github.com/DirtyHarryLYL/DJ-RN/blob/master/script/generate_3D_obj_GT.py