4.1. Trayectorias ortogonales.

La familia uniparamétrica de curvas f(t,y,a) = 0 es es la solución general de la ecuación diferencial F(t,y,y') =

∑f ∑f
0, que se obtiene eliminando el parámetro a en el sistema de ecuaciones f(t,y,a) = 0, ÅÅÅÅ∑ÅÅÅÅt Å + ÅÅÅÅ
∑ ÅyÅÅÅÅ y' =0,
Las trayectorias ortogonales a esta familia de curvas se obtienen hallando la solución de la ecuación diferencial
F(t,y,- ÅÅÅÅy'1Å ) = 0.

Ejemplo 4.1.1.
Determinar las trayectorias ortogonales a la familia de curvas que satisface la ecuación diferencial
1
y −  
2 xy = 0

Solución
Definimos E.D.O. de la familia de curvas.
1
ecu411 = y@xD −  x y @xD  0
2
1
y@xD −  x y @xD  0
2

La ecuación diferencial que verifican las trayectorias ortogonales a la familia de curvas dadas se obtiene sustituy-
endo y' por 1/y' en la ecuación diferencial anterior.
ecuorto411 = ecu411 ê. 8y '@xD → −1 ê y '@xD<
x
y@xD +   0
2 y @xD

Resolvemos esta E.D.O.

DSolve@ecuorto411, y@xD, xD
è!!!!!!!!!!!!!!!!!!!!!!!!!!!!! è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
−x2 + 4 C@1D −x2 + 4 C@1D
::y@xD → − 

è!!!!  
 >, :y@xD →  >>

è!!!!
2 2

Dado que Mathematica, siempre que puede, escribe las soluciones en forma explícita, han salido dos soluciones,
que podemos englobar en la expresión implícita x2 + 2 y2 = b, que engloba a todas las soluciones reales y en la que
hemos denotado 4 C[1] = b.
Esta es la familia de trayectorias ortogonales.

Ejemplo 4.1.2.
a) Hallar la ecuación diferencial de la familia uniparamétrica de curvas y 2 + 2 a x = a2 .
b) Hallar la ecuación diferencial de la familia de trayectorias ortogonales a las curvas dadas.
c) Dibujar unas pocas curvas de la familia, para valores de a entre -4 y 4.

Solución
apartado a.
Definimos la familia con y función de x, borrando las posibles asignaciones.
Clear@yD
2 Untitled-6

fam412 = y@xD ^ 2 + 2 ∗ a ∗ x  a ^ 2

2 a x + y@xD2  a2

Hallamos la derivada
D@fam412, xD
2 a + 2 y@xD y @xD  0

Despejamos la constante a de la última expresión

constante = Solve@D@fam412, xD, aD
88a → −y@xD y @xD<<

Sustituimos dicho valor en la familia. Obtendremos así la E.D.O. asociada a la familia.

ecu412 = Simplify@fam412 ê. constante@@1DDD

y@xD I2 x y @xD + y@xD I−1 + y @xD2 MM  0

apartado b.
La ecuación diferencial que verifican las trayectorias ortogonales a la familia de curvas dadas se halla sustituyendo
y' por 1/y'
ecuorto412 = ecu412 ê. 8y '@xD → −1 ê y '@xD<

y@xD i i + 1 y
j 2x y
jy@xD j
j−1 z
z −  z
z0
k k y @xD2 { y @xD {

Nota: Multiplicando esta E.D.O. por y 2 se obtiene la E.D.O. inicial.

En efecto
FullSimplify@−ecuorto412@@1DD ∗ y '@xD ^ 2 − ecu412@@1DDD
0

Luego la E.D.O. que verifica la familia de trayectorias ortogonales es la misma que verifica la familia dada. Esta
familia es una familia de curvas autoortogonales. Por cada punto pasan dos curvas de la familia, correspondientes a
diferentes valores del parámetro a, que son ortogonales.

apartado c.
Dibujamos las curvas pedidas, generadas a través de la instrucción Table, utilizando para ello el paquete Implicit-
Plot. Cargamos el paquete
<< Graphics`ImplicitPlot`

representamos las curvas

Untitled-6 3

ImplicitPlot@Table@y ^ 2 + 2 a ∗ x  a ^ 2, 8a, −4, 4<D, 8x, −3, 3<, AspectRatio → 1D;

Solve::svars : Equations may not give solutions for all "solve" variables. More…

ImplicitPlot::epfail : Equation y2  0 could not be solved for points to plot.

6
4
2

-3 -2 -1 1 2 3
-2
-4
-6

Nota: Para ver que es una familia autoortogonal representamos las curvas que pasan por un punto del plano por
ejemplo el (0,2).
Hallamos los valores del parámetro a, para que las curvas pasen por dicho punto.
Solve@fam412 ê. 8x → 0, y@xD → 2<, aD
88a → −2<, 8a → 2<<

Representamos las curvas

ImplicitPlot@8Hy ^ 2 + 2 a ∗ x  a ^ 2L ê. a → 2, Hy ^ 2 + 2 a ∗ x  a ^ 2L ê. a → −2 <,
8x, −3, 3<, AspectRatio → 1D;
4

-3 -2 -1 1 2 3
-2

-4

Son ortogonales.

Ejemplo 4.1.3.
Representar algunas curvas de de la familia de curvas ortogonales a la familia de hipérbolas x2 − a2 y 2 = a2 .

Solución
Definimos la familia con y función de x, borrando las posibles asignaciones.
Clear@yD

fam413 = x ^ 2 − a ^ 2 ∗ y@xD ^ 2  a ^ 2

x2 − a2 y@xD2  a2

Hallamos la derivada
D@fam413, xD

2 x − 2 a2 y@xD y @xD  0

Despejamos la constante a de la última expresión

4 Untitled-6

constante = Solve@D@fam413, xD, aD

è!!!! è!!!!
x x
::a → −  ! >, :a → 
è!!!!!!!!!!!!
è!!!!!!!!!!!!! è!!!!!!!!!!!! ! >>

è!!!!!!!!!!!!!
y@xD y @xD y@xD y @xD

Sustituimos dicho valor en la familia. Escogemos la primera solución pues en la ecuación la constante aparece
elevada al cuadrado. Obtendremos así la E.D.O. asociada a la familia.
ecu413 = Simplify@fam413 ê. constante@@1DDD
y@xD x
x Jx − 
 N  

y @xD y@xD y @xD

La ecuación diferencial que verifican las trayectorias ortogonales a la familia de curvas dadas se halla sustituyendo
y' por 1/y'
ecuorto413 = Simplify@ecu413 ê. 8y '@xD → −1 ê y '@xD<D
y @xD
x Jx +  + y@xD y @xDN  0
y@xD

Resolvemos esta E.D.O.

solorto4131 = DSolve@ecuorto413, y@xD, xD

InverseFunction::ifun :
Inverse functions are being used. Values may be lost for multivalued inverses. More…

Solve::ifun : Inverse functions are being used by Solve, so some

solutions may not be found; use Reduce for complete solution information. More…
è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
99y@xD → − ProductLog@e−x2 +2 C@1D D =, 9y@xD → ProductLog@e−x2 +2 C@1D D ==

Definimos las soluciones y representamos algunas de ellas

solprim4131@x_D = solorto4131@@1, 1, 2DD ê. C@1D → c
è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
− ProductLog@e2 c−x2 D

solseg4131@x_D = solorto4131@@2, 1, 2DD ê. C@1D → c

è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ProductLog@e2 c−x2 D

Table@solprim4131@xD, 8c, 0, 3<D

è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! è!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
9− ProductLog@e−x2 D , − ProductLog@e2−x2 D , − ProductLog@e4−x2 D , − ProductLog@e6−x2 D =

graf4131 = Plot@Evaluate@Table@solprim4131@xD, 8c, 0, 2, .5<DD, 8x, −2, 2<D;

-2 -1 1 2
-0.25

-0.5

-0.75

-1

-1.25

-1.5

-1.75
Untitled-6 5

graf4132 = Plot@Evaluate@Table@solseg4131@xD, 8c, 0, 2, .5<DD, 8x, −2, 2<D;

1.75

1.5

1.25

0.75

0.5

0.25

-2 -1 1 2

Y ahora las dos

Show@graf4131, graf4132D;

1.5

0.5

-2 -1 1 2
-0.5

-1

-1.5

Para ver que la familia de curvas obtenidas es la familia ortogonal representaremos algunas de las curvas de la familia
original
graf4133 = ImplicitPlot@Table@x ^ 2 − a ^ 2 ∗ y ^ 2  a ^ 2, 8a, −1, 1, .5<D,
8x, −1.5, 1.5<, PlotStyle −> RGBColor@1, 0, 0DD;

Solve::svars : Equations may not give solutions for all "solve" variables. More…

-1.5 -1 -0.5 0.5 1 1.5

-1

-2
6 Untitled-6

Y ahora las curvas de las dos familias

Show@graf4131, graf4132, graf4133, AspectRatio → AutomaticD;

-2 -1 1 2

-1

-2

Las curvas son ortogonales