def fun_lin(x, a, b):
    y = a*x + b
    return y

p0 = 2.5
p1 = 3.0