Scherk Surface Wireframe

import matplotlib.pyplot as plt
import numpy as np

# Scherk surface
x = np.linspace(-2, 2, 50)
y = np.linspace(-2, 2, 50)
X, Y = np.meshgrid(x, y)

# Avoid division issues
cos_y = np.cos(Y)
cos_x = np.cos(X)
cos_y[np.abs(cos_y) < 0.01] = 0.01
cos_x[np.abs(cos_x) < 0.01] = 0.01

Z = np.log(np.abs(cos_y / cos_x))
Z = np.clip(Z, -3, 3)

fig = plt.figure(figsize=(10, 8), facecolor='#ffffff')
ax = fig.add_subplot(111, projection='3d', facecolor='#ffffff')

ax.plot_wireframe(X, Y, Z, rstride=2, cstride=2, color='#9C2007', linewidth=0.5, alpha=0.9)

ax.set_xlabel('X', fontsize=11, color='#374151', labelpad=10)
ax.set_ylabel('Y', fontsize=11, color='#374151', labelpad=10)
ax.set_zlabel('Z', fontsize=11, color='#374151', labelpad=10)
ax.set_title('Scherk Surface Wireframe', fontsize=14, color='#1f2937', fontweight='bold', pad=20)

ax.tick_params(colors='#6b7280', labelsize=8)
ax.xaxis.pane.fill = False
ax.yaxis.pane.fill = False
ax.zaxis.pane.fill = False
ax.xaxis.pane.set_edgecolor('#e5e7eb')
ax.yaxis.pane.set_edgecolor('#e5e7eb')
ax.zaxis.pane.set_edgecolor('#e5e7eb')
ax.grid(True, alpha=0.3, color='#9ca3af')

plt.tight_layout()
plt.show()