33 lines
533 B
Python
33 lines
533 B
Python
import pybertini as pb
|
|
import numpy as np
|
|
|
|
x = pb.Variable('x')
|
|
y = pb.Variable('y')
|
|
|
|
vg = pb.VariableGroup([x,y])
|
|
|
|
sys = pb.System()
|
|
|
|
sys.add_function(x-y)
|
|
sys.add_function(x**2 + y**2 - 1)
|
|
|
|
sys.add_variable_group(vg)
|
|
|
|
C = pb.multiprec.Complex
|
|
|
|
variable_values = np.array([C(0), C(0)])
|
|
|
|
result = sys.eval(variable_values)
|
|
print(result)
|
|
|
|
sys.homogenize()
|
|
sys.auto_patch()
|
|
|
|
|
|
solver = pb.nag_algorithm.ZeroDimCauchyAdaptivePrecisionTotalDegree(sys)
|
|
|
|
solver.solve()
|
|
|
|
for soln in solver.solutions():
|
|
|
|
print(sys.dehomogenize_point(soln)) |