import numpy as np
def f12(x, y):
return float(8*x*x+8*x+6*y-14)
def f21(x, y):
return float(13*x*x-7*y-14)
def maxfunction(x, y):
return float(max(abs(f12(x,y)), abs(f21(x,y))))
x1=float(0.01)
x2=float(-0.01)
delta = round(float(maxfunction(x1, x2)),3)
for x in range(-100, 101):
for y in range(-200, 202, 2):
if delta > maxfunction(x / float(100), y / float(100)):
delta = round(float(maxfunction(x / float(100), y / float(100))),3)
x1 =float(x / float(100))
x2 =float(y / float(100))
print("delta = ", delta, " x1 = " , x1, "x2 = ", x2)
aW1wb3J0IG51bXB5IGFzIG5wCmRlZiBmMTIoeCwgeSk6CiAgICByZXR1cm4gZmxvYXQoOCp4KngrOCp4KzYqeS0xNCkKZGVmIGYyMSh4LCB5KToKICAgIHJldHVybiBmbG9hdCgxMyp4KngtNyp5LTE0KQpkZWYgbWF4ZnVuY3Rpb24oeCwgeSk6CiAgICByZXR1cm4gZmxvYXQobWF4KGFicyhmMTIoeCx5KSksIGFicyhmMjEoeCx5KSkpKQp4MT1mbG9hdCgwLjAxKQp4Mj1mbG9hdCgtMC4wMSkKZGVsdGEgPSByb3VuZChmbG9hdChtYXhmdW5jdGlvbih4MSwgeDIpKSwzKQoKZm9yIHggaW4gcmFuZ2UoLTEwMCwgMTAxKToKICAgIGZvciB5IGluIHJhbmdlKC0yMDAsIDIwMiwgMik6CiAgICAgICAgaWYgZGVsdGEgPiBtYXhmdW5jdGlvbih4IC8gZmxvYXQoMTAwKSwgeSAvIGZsb2F0KDEwMCkpOgogICAgICAgICAgICBkZWx0YSA9IHJvdW5kKGZsb2F0KG1heGZ1bmN0aW9uKHggLyBmbG9hdCgxMDApLCB5IC8gZmxvYXQoMTAwKSkpLDMpCiAgICAgICAgICAgIHgxID1mbG9hdCh4IC8gZmxvYXQoMTAwKSkKICAgICAgICAgICAgeDIgPWZsb2F0KHkgLyBmbG9hdCgxMDApKQpwcmludCgiZGVsdGEgPSAiLCBkZWx0YSwgIiB4MSA9ICIgLCB4MSwgIngyID0gIiwgeDIp