Matematiikka1

import numpy as np
import math

def f(x,y):
    apu = math.sqrt((6-x)/(7-x))
    apu += math.sqrt(2*x*y/(x+y))
    apu += 2 * math.sqrt((8-2*y)/(8-y))
    apu -= 5 * math.sqrt(2/3)
    return apu

range1 = np.arange(0.00, 6.00, 0.01)
d = -100.0
for a in range1:
    range2 = np.arange(0.00, 4.00, 0.01)
    for b in range2:
        c = f(a,b)
        if c > d:
            d = c
            aa = a
            bb = b

print("Max=", d, " x=", aa, " y=", bb)