# Yoda! -- over 53,756 triangles.
from scipy import io
x = io.loadmat('yodapose.mat')
from sage.plot.plot3d.index_face_set import IndexFaceSet
V = x['V']; F3=x['F3']-1; F4=x['F4']-1
Y = IndexFaceSet(F3,V,color='green') + IndexFaceSet(F4,V,color='green')
Y = Y.rotateZ(-1); Y = Y.rotateX(-1)
show(Y, frame=False)

# 53756 vertices!
len(x['V'])

import matplotlib.pyplot as plt
import numpy as np

# Vector field function (example: a simple radial field)
def vector_field(x, y, z):
    u = x
    v = y
    w = z
    return u, v, w

# Create the grid
x, y, z = np.mgrid[-2:2:20j, -2:2:20j, -2:2:20j]

# Compute the vector field
u, v, w = vector_field(x, y, z)


# Plot the vector field
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.quiver(x, y, z, u, v, w, length=0.1)
ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Z')
plt.title('3D Vector Field')
plt.s()


# LaTeX equation for the vector field (example)
latex_equation = r'\mathbf{F}(x, y, z) = \langle x, y, z \rangle'
print(f"LaTeX equation: ${latex_equation}$")