surf_area function

Surface area calculation