Algoritmia - algritmo del eje

import math

def calcular_esfuerzo_cortante(HP, D, N):

    # Convertir HP a libras-plg/min

    T = 63000 * HP / N

    # Calcular el esfuerzo cortante permitido

    S = (16 * T / (3.1426 * D)) ** (1/3)

    return S

# Lectura de los datos

HP = float(input("Ingrese la cantidad de caballos de fuerza (HP): "))

D = float(input("Ingrese el diámetro del eje (plg): "))

N = float(input("Ingrese la velocidad rotacional (rev/min): "))

# Calcular el esfuerzo cortante permitido

esfuerzo_cortante_permitido = calcular_esfuerzo_cortante(HP, D, N)

# Imprimir el resultado

print("El esfuerzo cortante permitido por el material del eje es:", esfuerzo_cortante_permitido, "lb/plg^2")