import numpy as np

def f(x):
    return x**2 + 1

dx = 0.1
x_array = np.arange(-2, 2+dx, dx)
y_array = f(x_array)

min_y = min(y_array)
min_x = x_array[y_array.argmin()]

print('min_y =', min_y)
print('min_x =', min_x)