Solving Equations with Python

[newbieAuggie2019]

More on the problem...

I have been working also on other versions of the program. One of them, is the following one (going back to numpy, even if the precision and the ability to handle so many numbers or just long numbers, are not very good). The positive point of this version, is that at least prints (plots) the values (x, y), but not all of them are correct, as we can see clearly (check by running the program) that the values for 'y' seem to be correct, going through a range from '0' to '18' as line 7 indicates:

y = np.arange(0,18,.01)

but for the values of 'x', they all seem to be changed to '1'. I still have to find out why, and to correct it:

import numpy as np
import matplotlib.pyplot as plt

x = np.arange(4,12,.01)
bigx = len(x)
print("\nThe matrix 'x' has", bigx, "elements.")
y = np.arange(0,18,.01)
ynumbers = [y]
bigy = len(y)
print("The matrix 'y' has", bigy, "elements.")
mother_of_matrices = bigx*bigy
print("Therefore, the interaction of matrices 'x' and 'y' produces a huge matrix of:", mother_of_matrices, "elements.")
k = 10/(6+x)

print('\n\nThese are the values for "x":\n\n', x)
print('\n\nThese are the values for "y":\n\n', y)
print('\n\nThese are the values for "k":\n\n', k)

for y in ynumbers:
        validynumbers = list(filter(lambda k: (1-k) != 0 and k != 0, k))
        for k in k:
            x = ((y**(1/k)-(y-6)**(1/k))/6)**(k/(1-k))
            plt.scatter(x,y)
            plt.show()

I tried to avoid problems related to dividing by '0' (or values of 'k' that seem to produce powers and fractions difficult to be handled by numpy), by adding line 20:

validynumbers = list(filter(lambda k: (1-k) != 0 and k != 0, k))

which filters the values of 'k' equal to '0' or that makes (1-k) = 0.
Nevertheless, it still gives us a warning message (it lets the program run, but warns us about something, in this case, still problems with the value '0'):

Error:
Warning (from warnings module):
  File "C:/Users/User1/AppData/Local/Programs/Python/Python37/equations_16.py", line 22
    x = ((y**(1/k)-(y-6)**(1/k))/6)**(k/(1-k))
RuntimeWarning: divide by zero encountered in double_scalars

All the best,

---

I have just noticed, trying to close the matplot window, before closing the Python 3.7.4 Shell window, that another matplot window pops up, showing values for 'y' from around '6' to around '18' and for values for 'x' from around '5' to around '40', which I think they are still incorrect, as the values for 'x', supposedly go from '4' to '12'.