C1 surface interpolation with constraints