Method of constructing bivariate osculatory rational interpolation function