Modifica di Asse Trasversale Cerniera (sezione)

===='''Script Pyhton: Errore Cuspidale con Apertura a 3 mm e Cuspidi Piatte'''====
</blockquote><br />

<syntaxhighlight lang="python">
import numpy as np
import matplotlib.pyplot as plt

# Definition of the original data
original_center = np.array([11.56, 61.21])  # Original center of the circle
r = 26.38  # Radius of the original circle
theta = np.linspace(0, 2 * np.pi, 10)  # 10 points along the circumference

# Generation of original points on the circumference
points = np.array([
    original_center + r * np.array([np.cos(angle), np.sin(angle)])
    for angle in theta
])

# Definition of the radius for the position of dental cusps
cusp_radius = 75.0  # Distance of the cusp from the center of rotation

# Generation of the positions of the original dental cusps
cusp_positions = np.array([
    original_center + cusp_radius * np.array([np.cos(angle), np.sin(angle)])
    for angle in theta
])

# Simulation of the errors Δ
errors = np.linspace(0, 10, 10)  # Errors Δ from 0 to 10 mm, 10 steps
max_differences_apertura_3mm = []

# Calculation of the cusp error for a mandibular opening of 3 mm
for delta in errors:
    new_center = original_center + np.array([delta, 0])  # Shift along the x-axis
    
    # Calculate the new positions of the dental cusps with error Δ
    new_cusp_positions = np.array([
        new_center + cusp_radius * np.array([np.cos(angle), np.sin(angle)])
        for angle in theta
    ])
    
    # Calculate the spatial differences between the original and calculated positions
    differences = np.linalg.norm(cusp_positions - new_cusp_positions, axis=1)
    max_difference = np.max(differences) * (3 / 30)  # Scaling the error for a 3 mm opening out of 18 mm maximum
    max_differences_apertura_3mm.append(max_difference)
    
    print(f"Error Δ = {delta:.1f} mm: Max Difference in Cusp Positions (3 mm opening) = {max_difference:.2f} mm")

plt.plot(errors, max_differences_apertura_3mm, 'bo-', label='Opening 3 mm')
plt.xlabel('Error Δ (mm) in HA localization')
plt.ylabel('Max Difference in Cusp Positions (mm)')
plt.title('Cusp Error with 3 mm Opening and Flat Cusps')
plt.grid(True)
plt.legend()
plt.show()

</syntaxhighlight><blockquote>