#Importacion de librerias
from [Link] import *
import [Link] as plt
import numpy as np
from math import *
from sympy import *
from [Link] import *
from time import *
from scipy import *
asctime()
clock()
#Introduccion de datos
while True:
try:
β=float(input("introduzca el ángulo en grados β = "))
except ValueError:
print ()
print("Cantidad Incorrecta")
print ()
continue
break
while True:
try:
vi=float(input("introduzca la velocidad inicial en (m/s) vi = "))
except ValueError:
print ()
print("Cantidad Incorrecta")
print ()
continue
break
#Ejecucion de calculos
print ()
θ=((β*pi)/180)
print ("θ=",θ)
g=9.81 # m/s**2
print ()
a=[Link](θ)
b=((g)/((2*vi**2)*[Link](θ)**2))
print ()
ymax=(vi**2)*([Link](θ)*[Link](θ))/(2*g)
xmax=(vi**2)*([Link](2*θ))/(g)
print ("................................................................................")
print ("................................................................................")
tmax=(vi*[Link](θ))/(g)
tv=2*(tmax)
#Salida de datos calculados
print (str("La altura máxima alcanzada por el proyectil es: Ymax")+" = "+str(ymax)+" m")
print (str("El alcance máximo horizontal del proyectil es: Xmax")+" = "+str(xmax)+" m")
print ("................................................................................")
print ("................................................................................")
print ("La altura máxima (m) alcanzada por el proyectil es: Ymax =",format(ymax,".2f"))
print ("El alcance máximo horizontal(m) del proyectil es: Xmax =",format(xmax," .2f"))
print ("................................................................................")
print ("................................................................................")
print ("El tiempo máximo t1max (s) que alcanza el proyectil para el ángulo β es: t1max
=",format(tmax,".2f"))
print ("El tiempo de vuelo t1v(s) que alcanza el proyectil para el angulo β es: t1v
=",format(tv,".2f"))
print ("................................................................................")
print ("................................................................................")
# Definimos la ecuación de la trayectoria
def f(x):
return(a*x-b*x**2)
x=[Link](0,xmax,500)
#Agregamos cuadricula
[Link](True)
[Link](color = '0.5', linestyle = '--', linewidth = 1)
#Proyeccion de la curva
[Link](x, f(x), "m", linewidth = 2, label = (str(β)+"º"))
#nombres de los ejes
[Link]("Xmax")
[Link]("Ymax")
# mostramos en pantalla la gráfica
[Link]()