Isabel Lucı́a Constantino Preciado

Arturo Jáuregui Dı́az

José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

2.1 Segundos
a) Escriba un guión que calcule el numero de segundos, s, dados el número de horas, h, de
acuerdo a la fórmula:
s = 3600h

b) Use el guión para encontrar el número de segundos en 1.5, 12 y 24 horas.

Algoritmo 1: Segundos
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19

h=f l o a t ( i n p u t ( ) )
s =3600∗h
p r i n t ( ” Para ” , h , ”h , son ” , s , ” s ” )

Algoritmo 2: Segundos
>> Para 1 . 5 h , son 5 4 0 0 . 0 s
>> Para 1 2 . 0 h , son 4 3 2 0 0 . 0 s
>> Para 2 4 . 0 h , son 8 6 4 0 0 . 0 s

2.2 Masa esférica

a) Escriba un guión que calcule la masa de una esfera dados un radio r y una densidad de
masa ρ, de acuerdo a la fórmula:
4π
m = ( )ρr3
3
b) Use el guión para encontrar la masa de una esfera de acero con radio r = 1 mm, r = 1 m
y r = 10 m.
Algoritmo 3: Masa esférica
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19

1
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

from math import p i

p=7850
r=f l o a t ( i n p u t ( ) )
m=round ( ( ( 4 ∗ p i / 3 ) ∗ ( p∗ r ∗ ∗ 3 ) ) , 5 )
p r i n t ( ”La masa para r=” , r , ”m e s : ” , m, ” kg ” )

Algoritmo 4: Masa esférica

>> La masa para r= 0 . 0 0 1 m e s : 3 e −05 kg
>> La masa para r= 1 m e s : 3 2 8 8 2 . 0 0 3 1 1 kg
>> La masa para r= 10 m e s : 3 2 8 8 2 0 0 3 . 1 0 7 5 7 kg

2.3 Ángulo
a) Escriba una función que para un punto (x, y) devuelva el angulo θ del eje x usando la
formula:
y
θ = arctan( )
x
b) Encuentre los ángulos θ para los puntos (1, 1), (−1, 1), (−1, −1), (1, −1)
c) Cómo se podrı́a cambiar la función para regresar valores de θ en el rango [0, 2π]?
Algoritmo 5: Angulo
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19

from math import atan2

def theta (x , y ) :
#Esta funcion devuelve el angulo entre el eje -x y un punto dado
t h e t a=atan2 ( y , x )
print ( theta )
return ;

Algoritmo 6: Angulo
>> t h e t a ( 1 , 1 )
0.7853981633974483
>> t h e t a ( −1 ,1)
2.356194490192345
>> t h e t a ( −1 , −1)
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

−2.356194490192345
>> t h e t a (1 , −1)
−0.7853981633974483

Algoritmo 7: Angulo
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19

from math import atan2 , p i

def theta (x , y ) :
#Esta funcion devuelve el angulo entre el eje -x y un punto dado
t h e t a=atan2 ( y , x )
i f t h e t a <0:
p r i n t ( 2 ∗ p i+t h e t a )
else :
print ( theta )
return ;

2.4 Vector unitario

a) Escriba una función que devuelva el vector unitario bidimensional, (ux , uy ), correspondi-
ente al angulo θ en el eje x. Puede usar la fórmula

(ux , uy ) = (cos θ, sin θ)

Donde θ esta dado en radianes.

b) Encuentre los vectores unitarios para θ = 0, θ = π6 , θ = π3 , θ = π2 , θ = 3π
2
.
c) Reescriba la función para tomar en su lugar el argumento θ en grados.
Algoritmo 8: Vector unitario
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19

from math import cos , s i n

Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

def vector (x ) :
#Esta funcion devuelve un vector unitario bidimensional
u x=round ( c o s ( x ) , 7 )
u y=round ( s i n ( x ) , 7 )
vec =[ u x , u y ]
p r i n t ( vec )
return ;

Algoritmo 9: Vector unitario

from math import p i
>> v e c t o r ( 0 )
[1.0 , 0.0]
>> v e c t o r ( p i / 6 )
[0.8660254 , 0.5]
>> v e c t o r ( p i / 3 )
[0.8660254 , 0.5]
>> v e c t o r ( p i / 2 )
[0.0 , 1.0]
>> v e c t o r ( 3 ∗ p i / 2 )
[ − 0 . 0 , −1.0]

Algoritmo 10: Vector unitario

# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19

from math import cos , s i n , r a d i a n s

def vector (x ) :
#Esta funcion devuelve un vector unitario bidimensional
x=r a d i a n s ( x )
u x=round ( c o s ( x ) , 7 )
u y=round ( s i n ( x ) , 7 )
vec =[ u x , u y ]
p r i n t ( vec )
return ;
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

2.5 Graficando la distribución normal

La distribución normal, usualmente llamada la distribución Gaussiana, está dada por:
1 2 2
P (x; µ, σ) = √ e−(x−µ) /(2σ )
2πσ 2

Donde µ es el promedio y σ es la desviación estándar.

a) Haga la función (x, µ, σ que devuelva el valor de distribución normal P (x, µ, σ) como dado
por la fórmula.
b) Use esta función para graficar la distribución normal para −5 < x < 5 para µ = 0 y σ = 1.
c) Grafique la distribución normal para −5 < x < 5 para µ = 0 y σ = 2 y para σ = 0.5 en la
misma grafica.
d) Grafique la distribución normal para −5 < x < 5 para σ = 1 y µ = 0, 1, 2 en tres gráficas
sobrepuestas.
Algoritmo 11: Distribucion normal
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 01/02/19

from pylab import ∗

from m a t p l o t l i b import p y p l o t a s p l t
import numpy a s np

#(a
x=np . l i n s p a c e ( −5 ,5 , 1000 )
d e f P( x , mu, sigma ) :
pass
P=(1/(2∗ p i ∗ sigma ∗ ∗ 2 ) ) ∗ exp (( −( x−mu) ∗ ∗ 2 ) / ( 2 ∗ sigma ∗ ∗ 2 ) )
return P

#(b
y=P( x , 0 , 1 )
plt . plot (x , y)
p l t . t i t l e ( ” G r f i c a ( b” , f o n t s i z e =12)
p l t . show ( )

#(c
y2=P( x , 0 , 2 )
y3=P( x , 0 , 0 . 5 )
p l t . p l o t ( x , y2 )
p l t . p l o t ( x , y3 )
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

p l t . t i t l e ( ” G r f i c a ( c ” , f o n t s i z e =12)
p l t . show ( )

#(d
y4=P( x , 0 , 1 )
y5=P( x , 1 , 1 )
y6=P( x , 2 , 1 )
p l t . p l o t ( x , y4 )
p l t . t i t l e ( ” G r f i c a ( d1” , f o n t s i z e =12)
p l t . show ( )
p l t . p l o t ( x , y5 )
p l t . t i t l e ( ” G r f i c a ( d2” , f o n t s i z e =12)
p l t . show ( )
p l t . p l o t ( x , y6 )
p l t . t i t l e ( ” G r f i c a ( d3” , f o n t s i z e =12)
p l t . show ( )
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

Figure 1: Azul: sigma = 2 y mu = 0. Naranja: sigma = 0.5

Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

1
2.6 Graficando xn

La función f (x; n) está dada por:

f (x; n) = x−n
a) Haga una función que devuelva el valor de f (x; n).
b) Use esta función para graficar x1 , 1
x2
, y 1
x3
en la misma gráfica para −1 < x < 1.
Algoritmo 12: Distribucion normal
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 30/01/19
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

from m a t p l o t l i b import p y p l o t a s p l t
import numpy a s np
#a)
x = np . l i n s p a c e ( −10 00 , 10 00 , 10 00 )
def function (x , n ) :
pass
y = ( 1 / x ∗∗n )
return y
#b)
y1=f u n c t i o n ( x , 1 )
y2=f u n c t i o n ( x , 2 )
y3=f u n c t i o n ( x , 3 )
p l t . p l o t ( x , y1 )
p l t . p l o t ( x , y2 )
p l t . p l o t ( x , y3 )
p l t . t i t l e ( ” G r f i c a ” , f o n t s i z e =12)
p l t . show ( )

Verde : n=3; Naranja : n=2; Azul : n=1

Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

2.11 Leyendo datos

El archivo trajectory.dat contiene una lista de números:

t0 x0 y0
t1 x1 y1

... ... ...

tn xn yn
Correspondientes al tiempo ti medido en segundos, y a las posiciones xi y yi medidas en
metros para la trayectoria de un proyectil.
a) Lea y acomode el archivo de datos en los arreglos t, x, y y.
b) Grafique las posiciones x y y como una función del tiempo en dos gráficas una sobre la
otra.
c) Grafique la posición (x, y) del objeto en una gráfica con x y y en los dos ejes.
Algoritmo 13: Trayectoria

# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 04/02/19

import pandas a s pd
from m a t p l o t l i b import p y p l o t a s p l t
d a t o s=pd . r e a d c s v ( ” t r a j e c t o r y . c s v ” , h e a d e r =0)
#a
t=d a t o s . l o c [ : , [ ’ t ’ ] ]
x=d a t o s . l o c [ : , [ ’ x ’ ] ]
y=d a t o s . l o c [ : , [ ’ y ’ ] ]
#b
plt . plot ( t , x)
plt . plot ( t , y)
p l t . t i t l e ( ” G r f i c a de p o s i c i n r e s p e c t o a l tiempo ” )
p l t . x l a b e l ( ”Tiempo ( s ) ” )
p l t . y l a b e l ( ” P o s i c i n (m) ” )
p l t . show ( )
#c
plt . plot (x , y)
p l t . t i t l e ( ” G r f i c a de l a p o s i c i n ” )
p l t . x l a b e l ( ” D i s t a c i a (m) ” )
p l t . y l a b e l ( ” A l t u r a (m) ” )
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

p l t . show ( )

Rojo : x r e s p e c t o a t ; Azul : y r e s p e c t o a t
Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

2.12 Integración numerica de un set de datos

El archivo velocity.dat contiene una lista de numeros:

t0 v0
t1 v1

... ...

tn vn
Correspondientes al tiempo ti medido en segundos, y a la velocidad vi medida en metros por
segundo para la trayectoria de un proyectil.
a) Lea y acomode el archivo de datos en los arreglos t y v.
b) Grafique v(t) como una función del tiempo.
Para un set de datos ti yvi , se puede estimar la función correspondiente a la integral de v(t)
con respecto a t en los tiempos ti usando el esquema iterativo:

y(t1 ) ' y(t0 ) + v(t0 )(t1 − t0 )

y(t2 ) ' y(t1 ) + v(t1 )(t2 − t1 )

Isabel Lucı́a Constantino Preciado
Arturo Jáuregui Dı́az
José Abraham Morales Vidales Mecanica Rotacional, Fluidos y Calor
Ivana Estefanı́a Betancourt Navarro Dr. Jorge Luis Menchaca
Jorge Antonio Ruı́z Zavalza Tarea 1 Febrero 1, 2019

... ' ...

y(tn ) ' y(tn−1 ) + v(tn−1 )(tn − tn−1 )

Donde v(ti ) =0 v(i)0 y ti =0 t(i)0 . Se puede asumir que el movimiento empieza en y(t0 ) = 0.0
m en t = t0 .
plt.plot(x,y,”ro”) plt.ylabel(”V”) plt.xlabel(”t”) plt.show
c) Escriba un guión para calcular la integral del tiempo y(ti ) del set de datos usando esta
fórmula. Implemente usando un ciclo ”repetir”.
d) Grafique la posición y(t) y la derivada v(t) como funciones del tiempo en dos gráficas
sobrepuestas.
Algoritmo 14: Integración
# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 04/02/19

import pandas a s pd
from m a t p l o t l i b import p y p l o t a s p l t
data=pd . r e a d c s v ( ”VC. t x t ” , h e a d e r =1, d e l i m w h i t e s p a c e=True )
x=data . i x [ : , 0 ]
y=data . i x [ : , 1 ]
p r i n t ( data )

Algoritmo 15: Integración

# Equipo Spyders
# Isabel Lucia Constantino Preciado
# Arturo Jauregui Diaz
# Jose Abraham Morales
# Ivana Estefania Betancourt Navarro
# Jorge Antonio Ruiz Zavalza
# Creado el 04/02/19

plt . p l o t (x , y , ” ro ” )
plt . y l a b e l ( ”V” )
plt . xlabel (”t”)
plt . show