Universidad Nacional Mayor de San Marcos
Facultad de Ciencias Matemáticas
Castillo Melchor Julios Deciderio
1. Representar gráficamente las siguientes funciones:
a) 𝒇(𝒙) = 𝒔𝒆𝒏(𝒙), 𝒈(𝒙) = 𝒙𝟐 + 𝟑𝒙 en el intervalo [𝟎, 𝟐𝝅].
import numpy as np
import [Link] as plt
pi= [Link]
def f(x):
y = [Link](x)
return y
def h(x):
y = x**2 + 3*x
return y
x = [Link](0,2*pi,100)
y_1=f(x)
y_2=h(x)
[Link](2,2,1)
[Link](x,y_1)
[Link](r"$f(x)=sen(x)$")
[Link](2,2,2)
[Link](x,y_2)
[Link](r"$h(x)=x^2+3x$")
[Link]()
b) Dibujar el polinomio 𝒙𝟐 + 𝟓𝒙 − 𝟑 con 𝟐𝟎𝟎 puntos, en color rojo, con trazo
discontinuo, con título y rejilla, en el intervalo [𝟏𝟎, −𝟏𝟎].
import numpy as np
import [Link] as plt
x = [Link](-10,10,200)
y = x**2 + 5*x -2
[Link](x,y,"r--")
[Link]( "Gráfica del polinomio "r"$y=x^2+5x-2$")
[Link]("Eje X")
[Link]("Eje Y")
[Link](True)
[Link]()
c) Representar en [𝟎, 𝟔] y en la misma gráfica las funciones:
i) 𝒇(𝒙) = 𝟑𝒙𝒆𝟐 en color azul.
ii) 𝒈(𝒙) = 𝒔𝒆𝒏(𝒙 + 𝟑) en rojo y con trazo discontinuo con leyendas
apropiadas.
import numpy as np
import [Link] as plt
e=np.e
def f(x):
y= 3*x*e**x
return y
def h(x):
y= [Link](x+3)
return y
x=[Link](0,6,200)
y_1=f(x)
y_2=h(x)
[Link](x,y_1,"b",label=r"$f(x)=3xe^x$")
[Link](x,y_2,"r--",label=r"$h(x)=sen(x+3)$")
[Link]()
[Link]()
2. Graficar los siguientes campos escalares:
𝟏
a) 𝒇(𝒓) = , en dos dimensiones se puede escribir: 𝒓 = √𝒙𝟐 + 𝒚𝟐
𝒓
import [Link] as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
fig = [Link]()
ax = Axes3D(fig)
x = [Link](-5,5,100)
y = [Link](-5,5,100)
X, Y= [Link](x, y)
def f(x, y):
z = 1/[Link](x**2 + y**2)
return z
ax.plot_surface(X, Y, f(X, Y))
[Link]()
𝟏 𝟏
b) 𝒇(𝒙, 𝒚) = −
√𝒙𝟐+(𝒚−𝟏)𝟐 √𝒙𝟐 +(𝒚+𝟏)𝟐
import [Link] as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
fig = [Link]()
ax = Axes3D(fig)
x = [Link](-5,5,100)
y = [Link](-5,5,100)
X, Y= [Link](x, y)
def f(x, y):
z = (1/[Link](x**2 + (y-1)**2))-(1/[Link](x**2 + (y+1)**2))
return z
ax.plot_surface(X, Y, f(X, Y))
[Link]()
c) 𝒇(𝒙, 𝒚) = 𝒙 − 𝒚 + 𝟐
import [Link] as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
fig = [Link]()
ax = Axes3D(fig)
x = [Link](-5,5,100)
y = [Link](-5,5,100)
X, Y= [Link](x, y)
def f(x, y):
z = x-y+2
return z
ax.plot_surface(X, Y, f(X, Y))
[Link]()
d) 𝒇(𝒙, 𝒚) = √𝟏𝟔 − 𝒙𝟐 − 𝒚𝟐
import [Link] as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
fig = [Link]()
ax = Axes3D(fig)
x = [Link](-5,5,100)
y = [Link](-5,5,100)
X, Y= [Link](x, y)
def f(x, y):
z = [Link](16-x**2-y**2)
return z
ax.plot_surface(X, Y, f(X, Y))
[Link]()
3. Graficar los siguientes campos vectoriales:
⃗ = 𝟑𝒊̂ − 𝟓𝒋̂
a) 𝒗
import [Link] as plt
import numpy as np
x = [0]
y = [0]
u = [3]
v = [-5]
[Link](x, y, u, v, color="g",
units="xy", scale=1)
[Link](-10,10)
[Link](-10,10)
[Link]("Eje X")
[Link]("Eje Y")
[Link]()
[Link]()
⃗ =𝒓
b) 𝒗 ⃗
import [Link] as plt
import numpy as np
x,y = [Link]([Link](-
10,10,10),[Link](-10,10,10))
u = x
v = y
[Link](x, y, u, v, color="g")
[Link](-10,10)
[Link](-10,10)
[Link]("Eje X")
[Link]("Eje Y")
[Link]()
⃗𝒓
⃗ =
c) 𝒗 , en dos dimensiones, usando coordenadas cartesianas donde:
𝒓𝟐
⃗ = 𝒙𝒊̂ + 𝒚𝒋̂, donde 𝒗
𝒓 ⃗ puede escribirse como:
𝒓̂ ⃗
𝒓 𝒙𝒊̂ + 𝒚𝒋̂
⃗ =
𝒗 𝟐 = 𝟑= 𝟑
𝒓 𝒓
(𝒙𝟐 + 𝒚𝟐 )𝟐
import [Link] as plt
import numpy as np
x,y = [Link]([Link](-
10,10,20),[Link](-10,10,20))
u = x/((x**2+y**2)**(3/2))
v = y/((x**2+y**2)**(3/2))
[Link](x, y, u, v, color="r")
[Link](-10,10)
[Link](-10,10)
[Link]("Eje X")
[Link]("Eje Y")
[Link]()
𝟑𝒙𝒚 𝟐𝒚𝟐 −𝒙𝟐
d) ⃗ =
𝒗 𝒊̂ + 𝒋̂
𝒓𝟓 𝒓𝟓
import [Link] as plt
import numpy as np
x,y = [Link]([Link](-
30,30,30),[Link](-30,30,30))
u = (3*x*y)/(x**2+y**2)**(5/2)
v = (2*y**2-x**2)/(x**2+y**2)**(5/2)
[Link](x, y, u, v, color="g")
[Link](-13,13)
[Link](-13,13)
[Link]("Eje X")
[Link]("Eje Y")
[Link]()
⃗ = 𝒙𝒊̂ − 𝒚𝒋̂
e) 𝒗
import [Link] as plt
import numpy as np
x,y = [Link]([Link](-
10,10,10),[Link](-10,10,10))
u = x
v = -y
[Link](x, y, u, v, color="g")
[Link](-10,10)
[Link](-10,10)
[Link]("Eje X")
[Link]("Eje Y")
[Link]()
𝟏
⃗ =
f) 𝒗 (𝒊̂ − 𝒋̂)
√𝟐
import [Link] as plt
import numpy as np
x,y = [Link]([Link](-
10,10,10),[Link](-10,10,10))
u = 1/([Link](2))
v = -1/([Link](2))
[Link](x, y, u, v, color="g")
[Link](-10,10)
[Link](-10,10)
[Link]("Eje X")
[Link]("Eje Y")
[Link]()
𝒙𝒊̂−𝒚𝒋̂
g) ⃗ =
𝒗
√𝟐
import [Link] as plt
import numpy as np
x,y=[Link]([Link](-
10,10,10),[Link](-10,10,10))
u = x/([Link](2))
v = -y/([Link](2))
[Link](x,y,u,v,color="r")
[Link]("Eje X")
[Link]("Eje Y")
[Link](-10,10,20)
[Link](-10,10,20)
[Link]()
h) ⃗𝒗 = 𝟐𝒚𝒊̂
import [Link] as plt
import numpy as np
x,y=[Link]([Link](-
10,10,20),[Link](-10,10,20))
u = 2*y
v = 0
[Link](x,y,u,v,color="r")
[Link]("Eje X")
[Link]("Eje Y")
[Link](-10,10,20)
[Link](-10,10,20)
[Link]()
⃗ = 𝒙𝟐 𝒊̂ + 𝒚𝟐 𝒋̂
i) 𝒗
import [Link] as plt
import numpy as np
x,y=[Link]([Link](-
10,10,20),[Link](-10,10,20))
u = x**2
v = y**2
[Link](x,y,u,v,color="r")
[Link]("Eje X")
[Link]("Eje Y")
[Link](-10,10,20)
[Link](-10,10,20)
[Link]()
𝒙𝒊̂−𝒚𝒋̂
j) ⃗ =
𝒗
√𝒙𝟐 +𝒚𝟐
import [Link] as plt
import numpy as np
x,y=[Link]([Link](-
10,10,20), [Link](-10,10,20))
u = x/([Link](x**2 + y**2))
v = -y/([Link](x**2 + y**2))
[Link](x,y,u,v,color="r")
[Link]("Eje X")
[Link]("Eje Y")
[Link](-10,10,20)
[Link](-10,10,20)
[Link]()
⃗ = 𝒙𝒚𝒊̂ − 𝒙𝒋̂
k) 𝒗
import [Link] as plt
import numpy as np
x,y=[Link]([Link](-
10,10,20),[Link](-10,10,20))
u = x*y
v = -x
[Link](x,y,u,v,color="r")
[Link]("Eje X")
[Link]("Eje Y")
[Link](-10,10,20)
[Link](-10,10,20)
[Link]()
⃗ = 𝒄𝒐𝒔𝒙𝒊̂ + 𝒔𝒊𝒏𝒙𝒋̂
l) 𝒗
import [Link] as plt
import numpy as np
x,y= [Link]([Link](-
10,10,20),[Link](-10,10,20))
u = [Link](x)
v = [Link](y)
[Link](x,y,u,v,color="r")
[Link]("Eje X")
[Link]("Eje Y")
[Link](-10,10,20)
[Link](-10,10,20)
[Link]()