The basic idea behind the method above:
You'll see that we're given 4 data points, and our function f(x) is composed of the addition of 4 parts. You can think of each part separately. With the first part we concentrate on the first point (0,0.0499), we're looking for a polynomial that vanishes (is 0) at the other 3 points, that's why the numerator has
)
etc, so when I plug in

this part is zero. Now, we also require that at x = 0, f(x) = 0.0499, this is why there are fractions in the denominator, so when I plug in x = 0, the numerator and denominators all cancel out, leaving the coefficient 0.0499.
The other parts are constructed similarly, and it is easy to see why this works.