Applying Mercator projection to ECEF data

Let [x,y,z] be a data point lying on the surface of a wgs84Ellipsoid. I would like to find the lattitude and longitude of this data point in Mercator Projection. If I were to do this using the standard Equirectangular projection, I would simply run

[lat, lon, ~] = ecef2geodetic(wgs84Ellipsoid, x,y,z);

However, I can't find any existing code to do the same thing to get the coordinates in Mercator Projection. Thank you very much for your help!