I can plot the intersection line of two surfaces f,g by using: Contourplot3D[f[x,y,z]==0,MeshFunctions->{Function[{x,y,z,f},f[x,y,z]-g[x,y,z]]}] But I need to do some calculations based on the intersection line. Or I need 100 (for example) group data of (x,y,z) on the line. How can I do this?