Möbius Strip

import matplotlib.pyplot as plt
import numpy as np

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

# Möbius strip parametric equations
n = 120
u = np.linspace(0, 2*np.pi, n)
width = 0.4

x1 = (1 + width*np.cos(u/2)) * np.cos(u)
y1 = (1 + width*np.cos(u/2)) * np.sin(u)
z1 = width * np.sin(u/2)

x2 = (1 - width*np.cos(u/2)) * np.cos(u)
y2 = (1 - width*np.cos(u/2)) * np.sin(u)
z2 = -width * np.sin(u/2)

ax.plot(x1, y1, z1, linewidth=3, color='#F5276C')
ax.plot(x2, y2, z2, linewidth=3, color='#4927F5')

# Surface fill
for i in range(0, n, 2):
    ax.plot([x1[i], x2[i]], [y1[i], y2[i]], [z1[i], z2[i]], 
            color='#27F5B0', alpha=0.35, linewidth=0.7)

ax.set_xlabel('X', fontsize=11, color='#1f2937')
ax.set_ylabel('Y', fontsize=11, color='#1f2937')
ax.set_zlabel('Z', fontsize=11, color='#1f2937')
ax.set_title('Möbius Strip', fontsize=14, color='#1f2937', fontweight='bold', pad=15)

ax.xaxis.pane.fill = False
ax.yaxis.pane.fill = False
ax.zaxis.pane.fill = False
ax.xaxis.pane.set_edgecolor('#000000')
ax.yaxis.pane.set_edgecolor('#000000')
ax.zaxis.pane.set_edgecolor('#000000')
ax.tick_params(colors='#000000', labelsize=9)
ax.grid(True, alpha=0.5, linewidth=0.5)
ax.xaxis._axinfo['grid']['color'] = '#000000'
ax.yaxis._axinfo['grid']['color'] = '#000000'
ax.zaxis._axinfo['grid']['color'] = '#000000'

ax.view_init(elev=30, azim=45)
plt.tight_layout()
plt.show()