-
Notifications
You must be signed in to change notification settings - Fork 1
/
bour2.py
33 lines (24 loc) · 757 Bytes
/
bour2.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
# Define the range of parameters u and v
u = np.linspace(-2 * np.pi, 2 * np.pi, 200)
v = np.linspace(-2, 2, 200)
u, v = np.meshgrid(u, v)
# Define the constant 'a'
a = 1 # You can adjust this constant to change the scale of the surface
# Parametric equations for the Bour surface
x = u - (u**3) / 3 + a * u * v**2
y = -v + (v**3) / 3 - a * v * u**2
z = u**2 - v**2
# Create a 3D plot
fig = plt.figure(figsize=(10, 7))
ax = fig.add_subplot(111, projection='3d')
# Plot the surface
ax.plot_surface(x, y, z, cmap='viridis')
# Set labels and title
ax.set_title('Bour Surface')
ax.set_xlabel('X axis')
ax.set_ylabel('Y axis')
ax.set_zlabel('Z axis')
plt.show()