VIGILADA MINEDUCACIÓN

Departamento de Matemáticas
ESCUELA
COLOMBIANA Cálculo Vectorial 2023 - 2
DE INGENIERÍA Guı́a de trabajo I
JULIO GARAVITO
Elementos básicos en Wolfram Mathematica
UNIVERSIDAD

Consultar en el texto guı́a los elementos teóricos relativos a vectores, geometrı́a del espacio
y funciones vectoriales.

1. Objetivos
Utilizar un asistente computacional para:
Obtener los gráficos en el espacio de objetos simples: puntos, segmentos, vectores.
Gráficar funciones vectoriales y representar los vectores tangente unitario, normal unitario
y binormal en algunos puntos de una curva.
Obtener expresiones simplificadas para la torsión y la curvatura.
Revisar algunos elementos relativos a parametrización de regiones y presentar algunas
funciones elementales del programa Wolfram Mathematica para familiarizarnos con el tra-
bajo que se hará posteriormente en el curso.

2. Representación de objetos básicos en Wolfram Mathematica

El punto de coordenadas (a, b, c) se grafica de color rojo con un tamaño adecuado, ejecutan-
do la instrucción

Graphics3D[{Red,PointSize[0.1],Point[{a,b,c}]}]

El segmento que une los puntos (a, b, c) y ( p, q, r ) se representa graficamente usando

Graphics3D[Line[{{a,b,c},{p,q,r}}]]

Se pueden incluir varios puntos para obtener un poligono. Con la siguiente instrucción se
dibuja el vector que va desde el punto (a, b, c) hasta el punto ( p, q, r ).

Graphics3D[Arrow[{{a,b,c},{p,q,r}}]]

Con la sentencia Graphics3D[Text["X",a,b,c]] se ubica el texto X en la coordenadas (a, b, c).

1
Ejemplo
Las siguientes tres intrucciones permiten obtener tres vectores de color azul, se les asigna
los nombres ejex, ejey y ejez, el punto y coma al final de linea oculta el producto obtenido.

ejex=Graphics3D[{Blue,Arrow[{{-1,0,0},{5,0,0}}]}];
ejey=Graphics3D[{Blue,Arrow[{{0,-1,0},{0,5,0}}]}];
ejez=Graphics3D[{Blue,Arrow[{{0,0,-1},{0,0,5}}]}];

Ahora se obtiene el gráfico del punto con coordenadas (2, 3, 4), y se agrega el texto ”P” cerca
de este. Estos objetos se nombran como punto y P.

Punto=Graphics3D[{Red,PointSize[0.03],Point[{2,3,4}]}];
P=Graphics3D[Text["P",{2,3,4.3}]];

Enseguida se presentan lı́neas de construcción grises y discontinuas. Estos objetos se nom-

bran como linea1 y linea2.

linea1=Graphics3D[{Gray,Dashed,Line[{{2,3,4},{2,3,0}}]}];
linea2=Graphics3D[{Gray,Dashed,Line[{{2,0,4},{2,0,0},{2,3,0},{0,3,0},{0,3,4},
{2,3,4},{2,0,4},{0,0,4},{0,3,4}}]}];

Por último, se muestran los objetos descritos anteriomente en una solo gráfico.

Show[ejex,ejey,ejez,Punto,linea1,linea2,P]

Si desea ocultar la caja que la cubre, se puede agregar Boxed->False.

Show[ejex,ejey,ejez,Punto,linea2,linea1,P,Boxed->False]

3. Gráficos de ecuaciones
La gráfica (en tres dimensiones) de la ecuación f ( x, y, z) = 0 puede hacerse con la instruc-
ción en Wolfram Mathematica usando la función ContoturPlot3D, para el simbolo de igualdad
se debe incluir doble igual, “==”.

ContourPlot3D[f[x,y,z]==0,{x,a,b},{y,c,d},{z,e,f}]

Cuando se trata de ecuaciones paramétricas de la forma x = f (t), y = g(t), z = h(t), con

t ∈ [ a, b] la instrucción más adecuada es ParametricPlot3D.

ParametricPlot3D[{f(t),g(t),h(t)},{t,a,b}]

2
4. Representación de funciones vectoriales en Wolfram Mathematica
Para definir la función vectorial r(t) = hcos(t), sin (t), cos (2t)i, se puede ejecuta la siguien-
te instrucción

r[t_]:={Cos[t],Sin[t],Cos[2*t]}

El vector posición en t = π/2 se obtiene con la sentencia: r[Pi/2]. El gráfico de la función

r(t) se consigue con la siguiente instrucción, se le da color rojo y un grosor determinado.

ParametricPlot3D[r[t],{t,0,2Pi},PlotStyle->{Thickness->0.02,Red}]

El gráfico del vector posición en t = 2π/3, de color negro, se consigue con las siguientes
instrucciones. Thick presenta un grosor determinado. Con Arrowheads[0.1] se modifica
el tamaño de la cabeza del vector.

Graphics3D[{Black,Thick,Arrowheads[0.1],Arrow[{{0,0,0},r[2*Pi/3]}]}]

La siguiente instrucción permite ver el recorrido de la curva (animar), para esto se usa la
π π
función Manipulate; en ella la variable k satisface que 20 ≤ k ≤ 2π con incrementos de 20 .
El objeto que se manipula es la gráfica de la curva, en la cual t se define de 0 a k.

Manipulate[ParametricPlot3D[r[t],{t,0,k},PlotRange->{-1,1},
PlotStyle->{Thickness->0.01,Red}],{k,Pi/20,2*Pi,Pi/20}]

La siguiente instrucción muestra la gráfica y el vector posición en el punto terminal de la

función r(t); los parámetros para k son los mismos descritos antes.

Manipulate[Show[{ParametricPlot3D[r[t],{t,0,k},PlotRange->{-1,1}],
Graphics3D[{Arrow[{{0,0,0},r[k]}]}]}],{k,Pi/20,2*Pi,Pi/20}]

4.1. Vector Normal

Con ayuda de la función Norm se define el vector tangente unitario T(t), luego se presenta
su enunciado y se halla T( π3 ).

T[t_]:=r’[t]/Norm[r’[t]] T[t] T[Pi/3]

Es posible que el programa incluya unas funciones que dependen de la naturaleza del
parámetro. Para obviar estos valores absolutos, se puede usar Simplify y la función Assuming.

T[t_]:=Assuming[t>=0,Simplify[r’[t]/Norm[r’[t]]]]
T[t]

El gráfico del vector tangente unitario en t = 2π/3 se obtiene con la siguiente instrucción.

3
 Graphics3D[{Blue,Thick,Arrow[{r[2*Pi/3],r[2*Pi/3]+T[2*Pi/3]}]}]

Con los siguientes comandos se obtiene la gráfica de r(t) junto con el vector posición y el
vector tangente en t = π/2.

Show[{ParametricPlot3D[r[t],{t,0,2*Pi},PlotStyle->{Thickness->0.01,Red},
PlotRange->{-1,1}],Graphics3D[{Arrow[{r[0],r[Pi/2]}]}],
Graphics3D[{Blue,Thick,Arrow[{r[Pi/2],r[Pi/2]+T[Pi/2]}]}]}]

4.2. Superficies y su curva de interseción

Las superficies con ecuaciones 4z = 3x2 − y2 y 2z = 10 − x2 − 3y2 , se intersectan en una
curva definida por la función vectorial r(t) = h2 cos(t), 2 sin(t), −1 + cos2 (t)i, con t ∈ [0, 2π ].
Las siguientes instrucciones representan esta situación.

Show[{ContourPlot3D[{3*x^2-y^2==4*z,2*z==10-x^2-3*y^2},{x,-3,3},{y,-3,3},
{z,-3,5},Mesh->False],ParametricPlot3D[{2*Cos[t],2*Sin[t],-1+4Cos[t]^2},
{t,0,2Pi},PlotStyle->{Red,Thickness[0.015]}]}]

5. Parametrizaciones básicas
Presentamos algunos ejemplos sencillos sobre parametrización de curvas, superficies o sóli-
dos; Las siguientes tres parametrizaciones servirán como principal insumo.

Parametrización de un segmento

Para dos puntos a y b en R n , la función −

→
r (t) = (1 − t)a + tb, t ∈ [0, 1], parametriza el
segmento rectilı́neo que une a a con b.

Intercambiando a y b en esta fórmula se logra cambiar el punto inicial con el punto final.
La identidad cos2 t + sin2 t = 1 motiva el uso de la siguiente elección

Parametrización de la circunferencia unitaria

La circunferencia unitaria con ecuación es x2 + y2 = 1, puede parametrizarse en la forma
−
→r (t) = hcos(t), sin (t)i , t ∈ [0, 2π ].

y y
1 1

x x
−1 1 −1 1
−1 −1
Figura 1: Gráficos de r1 (t) = hcos(t), sin(t)i y r2 (t) = hsin(t), cos(t)i

4
 También se puede asignar x = sin(t) y y = cos(t), de esta manera se describe la misma
curva, la diferencia está en los puntos terminales o en la orientación.
La ecuación x2 − y2 = 1 representa una hipérbola cuyas ramas abren hacia el eje x, la si-
guiente idea es motivada por la relación cosh2 t − sinh2 t = 1, para todo t.

Parametrización de la hipérbola unitaria

Una rama de la hipérbola unitaria se puede parametrizar por medio de la función −

→
r (t) =
hcosh(t), sinh(t)i , con t ∈ [−1, 1].
La otra rama se representa utilizando la función −
→
r (t) = h− cosh(t), sinh (t)i .
y y
1 1
x x
−2 −1 1 −2 −1 1
−1 −1

Figura 2: Gráficos de r1 (t) = hcosh(t), sinh(t)i y r2 (t) = hsinh(t), cosh(t)i

Ejemplos de parametrización de curvas

y
El segmento de recta, iniciando en el B
punto A(1, 1) y termiando en B(3, 2), se 2 •
A
parametriza simplificando la expresión 1 •
(1 − t) (1, 1) + t(3, 2), esto define la función x
r(t) = h1 + 2t, 1 + ti , con t ∈ [0, 1]. 1 2 3
z
El segmento rectilı́neo con punto inicial
(2, −1, 0) y punto final (−1, 2, 2), se pa- 2 b

rametriza simplificando (2, −1, 0)(1 − t) +

t(−1, 2, 2). Esto define la función −
→ -2
r (t) = -2 2 y
b 2
h2 − 2t, −1 + 3t, 2ti , t ∈ [0, 1]. 4
x
y
x2 y2
La elipse con ecuación 9 + 4 = 1 se para- 2
metriza al hacer x = 3 cos(t) y y = 2 sin(t). Es x
decir, −
→
r (t) = h3 cos(t), 2 sin(t)i , t ∈ [0, 2π ]. −3 3
−2

5.1. Parametrización de regiones planas

y
La expresión 4x2 + 9y2 = 36 represen-
ta una elipse, la región interior se pue-
de describir mediante la función u(r, θ ) = −3 3 x
h3r cos θ, 2r sin θ i , con r ∈ [0, 1], θ ∈ [0, 2π ].

5
 El triángulo con vértices A, B y C, se pa- C
rametriza a través de segmentos. La función
r(u, v) = (1 − v)C + v ((1 − u) A + uB) , des-
B
cribe esta región, con u ∈ [0, 1], v ∈ [0, 1].
A
y Con las ecuaciones y = x + 1 y y = x2 − 1,
3 para x ∈ [−1, 2], se simplifica la expresión
(1 − u) h x, 1 + x i + u x, x2 − 1 , que define la
función r( x, u) = h x, ( x + 1)(1 − 2u + ux )i ,
x con x ∈ [−1, 2], u ∈ [0, 1], que determina la
2
región mostrada.

El cuadrilátero con vértices A(−2, −1), B(−1, 4), y

C(3, 1), D (4, 3), se parametriza (adecuada- B•4
D
mente) con segmentos (no crear un corbatı́n). •
Consideremos P = (1 − t) A + tB y Q = (1 − t)C + tD, 2 •Q
simplificando (1 − u) P + uQ, se obtiene r(t, u) = P•
•C
ht + 5u − 2, 5t + 2u − 3tu − 1i , con t ∈ [0, 1], u ∈ [0, 1].
−2 2 4 x
A•

Con las siguientes instrucciones se consigue el gráfico

ParametricPlot[(1-u)*((1-t)*{-2,-1}+t*{-1,4})+u*((1-t)*{3,1}+t*{4,3}),
{t,0,1},{u,0,1}]

5.2. Parametrización de superficies

Usando coordenadas esféricas, la superficie de la esfera unitaria se puede parametrizar co-
mo r(θ, φ) = hcos θ sin φ, sin θ sin φ, cos φi , con θ ∈ [0, 2π ] y φ ∈ [0, π ].

x2 y2 2
Con esto, el elipsoide 16 + 9 + z4 = 1, se pa-
rametriza por medio de la función r(θ, φ) =
h4 cos θ sin φ, 3 sin θ sin φ, 2 cos φi . El gráfico se ob-
tiene ejecutando las siguientes instrucciones.

ParametricPlot3D[{4Cos[\[Theta]]*Sin[\[Phi]],3*Sin[\[Theta]]*Sin[\[Phi]],
2*Cos[\[Phi]]},{\[Theta],0,2*Pi},{\[Phi],0,Pi},PlotStyle->Orange] %% $

Las expresiones z = x2 − y2 y x2 + y2 + 2y = 8 representan un paraboloide hiperbólico y un

cilindro elı́ptico. La ecuación del cilindro representa una elipse en el plano xy, el interior de esta
región se parametriza como x = 3u cos(t), y = −1 + 3u sin(t), con la ecuación del paraboloide
hiperbólico se halla z, esto es z = x2 − y2 = 9u2 cos(2t) + 6u sin(t) − 1.
Entonces r (t, u) = 3u cos(t), −1 + 3u sin(t), 9u2 cos(2t) + 6u sin(t) − 1 con t ∈ [0, 2π ], u ∈
[0, 1], parametriza la región sobre el paraboloide hiperbólico acotada por el cilindro.

6
Los gráficos se consiguen con las siguientes instrucciones.

a=ContourPlot3D[{z==x^2-y^2,x^2+y^2+2*y==8},{x,-5,5},{y,-5,5},{z,-14,10},
Mesh->False,ContourStyle->{Directive[Orange,Opacity[0.6]],
Directive[Cyan,Opacity[0.2]]},PlotPoints->60];
b=ParametricPlot3D[{3uCos[t],-1+3uSin[t],-1+9u^2Cos[2t]+6uSin[t]},
{t,0,2*Pi},{u,0,1},PlotStyle->{Cyan},PlotPoints->80];
Show[a,b]

Las ecuaciones x2 + (y − 3)2 = 4 y 4z = 36 − x2 + y2 representan un cilindro y un para-

boloide invertido. Consideremos la superficie σ sobre el paraboloide exterior al cilindro, con
z ≥ 0. El gráfico se obtienen con ayuda de las siguientes instrucciones.

ContourPlot3D[{z==9-(x^2+y^2)/4,x^2+(y-3)^2==4},{x,-6,6},{y,-6,6},{z,0,9},
Mesh->False,ContourStyle->{Orange,Cyan}]
ParametricPlot3D[{2*(2*u+1)*Sin[t],2*(2*u+1)Cos[t]-3*u+3,
(1-u)(23+25u-12(1+2u)Cos[t])/4},{u,0,1},{t,0,2*Pi}]

z
9 y
6
7
3 P•
5 •Q
x
3
-6 1 -6
−6 −3 3 6
-3 -3 −3
3 -1 3
6 6
x y −6

La base del cilindro se parametriza como x = 2 sin(t), y = 3 + 2 cos(t), y la intersección del

paraboloide con el plano z = 0 se parametriza como x = 6 sin(t), y = 6 cos(t), con t ∈ [0, 2π ].
Con el fin de parametrizar la proyección de σ sobre el plano xy, se forman segmentos entre
dos puntos P (2 sin(t), 3 + 2 cos(t), 0) y Q (6 sin(t), 6 cos(t), 0) , es decir (1 − u) P + uQ. La re-
gión sobre el plano xy acotada por las dos circunferencias presentadas antes se representa por
−
→v (u, t) = (2 (2u + 1) sin t, 2 (2u + 1) cos t − 3u + 3, 0) . La coordenada z sobre la superficie del
paraboloide se consigue al remplazar z en esa ecución. Por lo anterior, la superficie σ se para-
metriza, con t ∈ [0, 2π ], u ∈ [0, 1], mediante la función
1
r(t, u) = h2(2u + 1) sin t, 2(2u + 1) cos t − 3u + 3, (1 − u)(25u − 12(2u + 1) cos t + 23)i
4

7
5.3. Una esfera y un cilindro elı́ptico
Consideraremos la esfera x2 + y2 + z2 = 16 y el cilindro elı́ptico 4( x − 1)2 + 9y2 = 36. En
el plano xy, la región interior al cilindro se parametriza en la forma x − 1 = 3u cos(t), y =
2u sin(t), con u ∈ [0, 1] y t ∈ [0, 2π ]; para la frontera usaremos x = 1 + 3 cos(t), y = 2 sin(t).
z
5

3P
•
1
-3 -3
-11 -1
3 -11
x 5 3
5 y
-3
•
Q
q
Sobre la curva de intersección se satisface que z = ± 5 sin2 (t) − 6 cos(t) + 6, por lo tanto,
q
−
→r (t) = h1 + 3 cos(t), 2 sin (t), ± 5 sin2 (t) − 6 cos(t) + 6i, t ∈ [0, 2π ], parametriza las cur-
2
vas en mención. Sobre la esfera y dentro q el cilindro, el valor de z satisface (1 + 3u cos(t)) +
(2u sin(t)) 2 + z2 = 16, luego z = ± 15 − 6u cos(t) − 9u2 + 5u2 sin2 (t), por tanto, una para-
metrización de la parte de la esfera dentro del cilindro es, para t ∈ [0,2π ], y u ∈ [0, 1],
q
r(u, t) = 1 + 3u cos(t), u sin (t), 15 − 6u cos(t) − 9u2 + 5u2 sin2 (t) .
q
Denotemos a z1 = 6 − 6 cos(t) + 5 sin2 (t), y consideremos dos puntos sobre el cilindro
P = (1 + 3 cos(t), 2 sin (t), z1 ) y Q = (1 + 3 cos(t), 2 sin(t), −z1 ), con ellos se parametriza el seg-
mento rectilı́neo que los une. Simplificando
 la expresión (1 − u) P + uQ, con u ∈ [0, 1], se obtie-

q
2
ne la expresión la función r(u, t) = 1 + 3 cos(t), 2 sin(t), (1 − 2u) 6 − 6 cos(t) + 5 sin (t) ,
con t ∈ [0, 2π ], u ∈ [0, 1]. Los gráficos se obtienen con ayuda de las siguientes instrucciones

Show[{ParametricPlot3D[{{1+3*u*Cos[t],2*u*Sin[t],Sqrt[15-6*u*Cos[t]-9*u^2
+5*u^2*(Sin[t])^2]}},{u,0,1},{t,0,2Pi},PlotStyle->Red],ContourPlot3D[{x^2+
y^2+z^2==16,(x-1)^2/9+y^2/4==1},{x,-4,4},{y,-4,4},{z,-4,4},ContourStyle->
{Directive[Orange,Opacity[0.15]],Directive[Gray,Opacity[0.15]]},Mesh->False]}]

Show[ContourPlot3D[{x^2+y^2+z^2==16,(x-1)^2/9+y^2/4==1},{x,-4,4},{y,-4,4},
{z,-4,4},ContourStyle->{Directive[Orange,Opacity[0.15]],Directive[Gray,
Opacity[0.15]]},Mesh->False],ParametricPlot3D[{{1+3*Cos[t],2*Sin[t],
(1-2*u)Sqrt[6-6*Cos[t]+5*(Sin[t])^2]}},{u,0,1},{t,0,2Pi},PlotStyle->Red]]

Sólido acotado por dos paraboloides

Los paraboloides elı́pticos con ecuaciones z = 3 − 2x2 − y2 y z = x2 + 2y2 , determinan una
región que se parametriza hallando la curva de intersección que corresponde a x2 + y2 = 1. El
interior de esta región se representa con x = u cos v, y = u sin v con u ∈ [0, 1] y v ∈ [0, 2π ].

8
 z Sean P y Q dos puntos sobre cada uno de los paraboloides,
en ellos, la coordenada z corresponde, respectivamente, a
3 2
z = 3 − 2x2 − y2
b
Q z1 = 12 (6 − 3u2 − u2 cos(2v)) y z2 = u2 (3 − cos(2v)).
2
2 Uniendo los puntos P(u cos v, u sin v, u2 (3 − cos(2t))) y
Q(u cos v, u sin v, 12 (6 − 3u2 − u2 cos(2t))) mediante seg-
1
z = x2 + 2y2 mentos rectilı́neos, la expresión (1 − w) P + wQ, se obtiene
-1 -1b la función que parametriza el sólido
1
P 1 2
2 y r(u, v, w) = hu cos v, u sin v, 3(1 − w) + u2 (6w − 3 −
x cos(2t))i, u ∈ [0, 1], v ∈ [0, 2π ], w ∈ [0, 1].
Una representación de esta región se obtiene con las siguientes instrucciones.

RegionPlot3D[x^2+2*y^2<=z<=3-2*x^2-y^2&&x^2+y^2<=1,{x,-1,1},{y,-1,1},{z,0,3},
PlotPoints->80,Mesh->False,PlotStyle->Orange]

6. Ejercicios
1. Para las siguientes funciones vectoriales, obtenga el gráfico de la curva y para tres puntos
arbitrarios, represente el vector Tangente unitario, Normal unitario y Binormal.
a) r(t) = hsin(t) + 2 sin(2t), cos(t) − 2 cos(2t), − sin (3t)i

b) r(t) = hcos(2t), sin (2t), cos(4t)i d) r(t) = hcos2 (t), sin2 (t), sin (2t) + cos(t)i
c) r(t) = he−t cos(t), e−t sin(t), e−t sin(t)i e) r(t) = hcos(t), 2 sin(t), ti

2. Con las funciones vectoriales del numeral anterior, obtener una representación (dinámica)
gráfica de los vectores T, N y B, en puntos diferentes;
Sugerencia: Las funciones Manipulate o Animate permiten hacerlo.
3. Utilizar las funciones Dot y Cross que permiten evaluar el producto punto y el producto
cruz para hallar la curvatura y la torsión de las siguientes funciones vectoriales.

a) r(t) = ht, 1t + 1, 1t − ti d) r(t) = hcosh(t), sinh (t), ti

3 3 √
b) r(t) = h 5t3 − t2 + 1, 5t3 + t2 − 1, √t + 3i e) r(t) = ht2 , 2 3t, 3 ln(t)i
5
c) r(t) = hcos(t), sin (t), cos (2t)i f ) r(t) = hsec(t), ln (tan(t) + sec(t)), tan (t)i

Parametrice y grafique en Mathematica las siguientes regiones.

4. La región del plano comprendida entre las y
x2 y2
x2 y2 y2
3
16 + 9 =1
elipses 16 + 9 = 1 y ( x − 1 )2 + 4 = 1.
x
−4 −2 2 4
y2
−3 ( x − 1 )2 + 4 =1

9
5. Parametrice la región sobre el paraboloi- 9. La porción del cilindro elı́ptico con ecua-
de z = 4 − x2 − y2 acotada por el plano ción 4x2 + 9y2 = 36 acotado por los planos
x − y + z = 1; también parametrice la re- 3z + x = 9 y 3z − y = −9.
gión sobre el plano acotada por el parabo-
loide y el sólido acotado por las superfi-
cies anteriores.

6. Los paraboloides 4y = z2 − x2 y y =
4 − x2 − z2 se intersecan en una curva C;
parametrice C y las superficies sobre cada
uno de los paraboloides acotadas por C.

10. El sólido acotado por el cilindro elı́pti-

co 2x2 + y2 = 16 y los paraboloides hi-
perbólicos x2 − y2 = 4z + 12 y x2 − 2y2 =
4z − 20.

7. La superficie obtenida cuando la región

2
x2
del plano acotada por las elipses 16 + y9 =
y2
1 y ( x − 1)2 + 4 = 1 se proyecta sobre el 11. El sólido acotado por dos cilindros de
igual radio con ejes de simetrı́a perpen-
paraboloide hiperbólico 3z = y2 − x2 − 2y.
diculares. En particular, se puede supo-
ner que poseen ecuaciones x2 + z2 = r2
z z
y x2 + y2 = r2 , con r > 0.
3 3
-3 -3
6
3 y x 3
y
x -3 -3

8. El sólido acotado por las gráficas del plano

p + x + y = 20 y el semicono z =
5z
x 2 + y2 .
Figura 3: x2 + z2 ≤ r2 & x2 + y2 ≤ r2 .

12. Considere una esfera (E) de radio a y un

cilindro circular recto (C) de radio b, cuyo
eje de simetrı́a contiene un diametro de E,
además 0 < b < a.
Asigne valores para a y b, luego parame-
trice el sólido interior a E y exterior a C.

10