Hello I have one function that creates polynomials (i.e legendre.polynomials) I want to use this one to create polynomials for variable x and variable y. legendre.polynomials(2) [[1]] 1 [[2]] x [[3]] -0.5 + 1.5*x^2 the ideal would be to receive the same output but for another variable (eg. y) Then having two equation with different x and y I can create the multinomial I want to. I checked the multipolynom package but it can not do what I am looking for. Best Regards Alex