29 lines
677 B
Python
29 lines
677 B
Python
import pybertini as pb
|
|
import numpy as np
|
|
|
|
pb.default_precision(500)
|
|
|
|
x = pb.Variable('x')
|
|
y = pb.Variable('y')
|
|
z = pb.Variable('z')
|
|
|
|
vg = pb.VariableGroup([x,y,z])
|
|
|
|
sys = pb.System()
|
|
|
|
sys.add_function(x-y)
|
|
sys.add_function(x**2 + y**2 - 1)
|
|
sys.add_function(5*x**3 + 16*x*y**4 - 17*x*y*z - z**3)
|
|
|
|
sys.add_variable_group(vg)
|
|
|
|
random_complex = lambda : pb.random.complex_in_minus_one_to_one()
|
|
|
|
n_iterations = 10000
|
|
print(f'evaluating system at random complex point {n_iterations} times at precision {pb.default_precision()}')
|
|
for ii in range(n_iterations):
|
|
variable_values = np.array([random_complex(), random_complex(), random_complex()])
|
|
result = sys.eval(variable_values)
|
|
|
|
|