With all numpy floating-point errors ignored, compute and return np.ones(1) / 0 (which is array([inf])).
np.ones(1) / 0
array([inf])
Sign in to write and run your own code.
solve(...)