f = @(x)_b9+b1*(x-x®);
b0 = 420.0;
x0 = 20.0;
x1 = 25.5;
fx0 = 420.0;
fx1 525.5;
b1 = (fx1-fx0)/(x1-x0);
ans = f(21);
ans