Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

There are several little problems in your code.

First you do not need to use SR in build_variables (avoid SR as much as possible when using polynomials). You can do

def build_variables(n):
    names = []

    for i in range(n):
        for j in range(n):
            names += ['x_' + str(i) + '_' + str(j)]


Then in the lambda_siep function, you should only use the result my_variables of build_variables once, inside the PolynomialRing constructor.

Then you should not use my_variables anymore, but instead your X in the line

J[i,j] = derivative(Y[0][i], X[j][j])

The two lines


are useless and can be removed. In particular, inject_variables is only to be used in interactive mode.