Numerical approximations to extremal metrics on toric surfaces