Bessel Function Wireframe

import matplotlib.pyplot as plt
import numpy as np
from scipy.special import jv

# Bessel function J0
x = np.linspace(-10, 10, 70)
y = np.linspace(-10, 10, 70)
X, Y = np.meshgrid(x, y)
R = np.sqrt(X**2 + Y**2)
Z = jv(0, R)

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='#5314E6', linewidth=0.4, 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('J₀(r)', fontsize=11, color='#374151', labelpad=10)
ax.set_title('Bessel Function 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()