Tengo la siguiente inquietud con sistema de ecuaciones no lineales en python:
from scipy.optimize import fsolve as fs
def nonlinear(z):
'''Solve a sistems 2x2 not linear
'''
x, y = z[0], z[1]
f1 = (x**2)+(y**2)-2*(4.41*x+2.68*y)+25.59
f2 = (x**2)+(y**2)-2*(3.23*x+2.1*y)+14.49
return [f1, f2]
¿Cómo puedo encontrar las múltiples soluciones del sistema anterior con el código que realicé?
Por que si ejecuto la función anterior con unos valores iniciales me arrojaría sólo una solución:
x, y = fs(nonlinear, [1,1])
print(x,y)
3.38998152293 2.67210655679
Todas las soluciones que tiene el sistema son las siguientes:
f1
yf2
. Qué tipo de sistema es ese y cómo quieres resolverlo. Si a lo que te refieres es a obtener las dos raíces de cada incógnita, deberías resolver de forma explícita para las dos raíces.