THE CALCULUS OF BIVARIATE FRACTAL INTERPOLATION SURFACES