Métodos de primer orden · Tema 17

Ecuaciones exactas

Cómo reconocer una diferencial total, reconstruir su función potencial e interpretar las soluciones como curvas de nivel.

01 · Una forma diferente de mirar la EDO

La ecuación puede ser la diferencial de una función

Una ecuación de primer orden puede escribirse como

M(x,y)dx+N(x,y)dy=0

Es exacta si existe una función potencial F(x,y) cuya diferencial total coincide con el miembro izquierdo:

dF=Fxdx+Fydy=Mdx+Ndy

Por tanto, necesitamos Fx=M y Fy=N. Si la diferencial de F se anula a lo largo de una trayectoria, F permanece constante:

F(x,y)=C

La solución aparece naturalmente en forma implícita como una familia de curvas de nivel.

02 · Regla de la cadena

Por qué F(x,y)=C describe soluciones

Si una curva se escribe como y=y(x), entonces la regla de la cadena da:

d/dx F(x,y(x))=Fx+Fyy′

La ecuación diferencial dividida por dx es

M(x,y)+N(x,y)y′=0

Cuando Fx=M y Fy=N, ambas expresiones coinciden:

d/dx F(x,y(x))=0

Así, F es constante a lo largo de cada curva solución. Cuando N≠0, también podemos despejar la pendiente:

y′=−M/N

La forma diferencial sigue teniendo sentido geométrico incluso donde la tangente es vertical y no conviene describir la curva como y(x).

03 · Prueba de exactitud

Las derivadas cruzadas deben coincidir

Si M=Fx y N=Fy, y F tiene segundas derivadas continuas, entonces

My=Fxy=Fyx=Nx

Por eso My=Nx es una condición necesaria local. En una región abierta y simplemente conexa, con derivadas continuas, también es suficiente.

Coinciden

My=Nx

En una región adecuada, existe un potencial y se aplica el método de ecuaciones exactas.

No coinciden

My≠Nx

La expresión no es una diferencial total tal como está escrita.

04 · Reconstrucción sistemática

Integrar una componente y completar la otra

Para resolver Mdx+Ndy=0:

  1. Declarar el dominio.

    Identificar denominadores, logaritmos, raíces y posibles agujeros.

  2. Comprobar exactitud.

    Calcular My y Nx.

  3. Integrar M respecto de x.

    F(x,y)=∫M(x,y)dx+g(y).

  4. Derivar respecto de y.

    Calcular Fy y compararlo con N.

  5. Determinar g.

    Resolver g′(y)=N−∂y∫Mdx.

  6. Escribir la familia.

    F(x,y)=C.

  7. Aplicar datos y comprobar.

    Evaluar C y verificar Fx=M, Fy=N.

05 · Ejemplo completo

Reconstruir un potencial cuadrático

Consideremos

(2x+y)dx+(x+2y)dy=0

Identificamos M=2x+y y N=x+2y. Las derivadas cruzadas son:

My=1,    Nx=1

La ecuación es exacta en todo ℝ². Integramos M respecto de x manteniendo y constante:

F=∫(2x+y)dx=x²+xy+g(y)

Derivamos respecto de y y comparamos:

Fy=x+g′(y)=x+2y=N

Así, g′(y)=2y y podemos elegir g(y)=y². La solución general es:

x²+xy+y²=C

Las curvas de nivel para C>0 son elipses rotadas.

06 · Seleccionar una curva de nivel

La condición inicial determina C

Para la ecuación anterior, si y(1)=2, el punto (1,2) debe pertenecer a la curva:

C=1²+1·2+2²=7

La solución implícita del PVI es:

x²+xy+y²=7

Si despejamos y, obtenemos dos ramas:

y(x)=[−x±√(28−3x²)]/2

La condición y(1)=2 selecciona el signo positivo cerca de x=1. La rama solo puede describirse como función mientras N=x+2y≠0; cuando N=0, la curva tiene tangente vertical.

07 · Camino alternativo

También se puede comenzar integrando N

Para el mismo ejemplo:

F=∫(x+2y)dy=xy+y²+h(x)

Ahora derivamos respecto de x:

Fx=y+h′(x)=2x+y=M

Entonces h′(x)=2x y h=x². Llegamos al mismo potencial, salvo una constante aditiva irrelevante.

Conviene integrar primero la componente más sencilla. El resultado debe ser compatible con ambas derivadas parciales; esa doble comprobación detecta signos o términos omitidos.

Un término no puede desaparecer

Si al comparar Fy con N aparece una expresión que todavía depende de x donde debería estar g′(y), algo salió mal o la ecuación no era exacta.

08 · Cuando falla el criterio

Un factor integrante puede crear exactitud

La forma

Mdx+Ndy=0

puede no ser exacta, pero volverse exacta al multiplicarla por una función no nula μ(x,y). Encontrar un factor general puede ser difícil; hay dos pruebas especiales útiles.

Criterios para factores integrantes de una sola variable.
CocienteSi depende solo de...Factor
(My−Nx)/Nxμ(x)=exp(∫[(My−Nx)/N]dx)
(Nx−My)/Myμ(y)=exp(∫[(Nx−My)/M]dy)

Ejemplo breve

Para y dx+x²dy=0, M=y, N=x², de modo que My=1 y Nx=2x: no es exacta. Sin embargo,

(My−Nx)/N=(1−2x)/x²

depende solo de x. En un intervalo que no cruce cero:

μ(x)=exp[∫(1/x²−2/x)dx]=e−1/x/x²

Al multiplicar, la forma se vuelve exacta y un potencial es F(x,y)=ye−1/x. Por tanto, la familia queda ye−1/x=C.

09 · Geometría y topología

Gradiente perpendicular y dominios con agujeros

El vector ∇F=(M,N) es perpendicular a las curvas F=C. Un vector tangente es (N,−M), porque

(M,N)·(N,−M)=MN−NM=0

Su pendiente es −M/N, exactamente la exigida por la ecuación.

La igualdad de derivadas cruzadas es local. En el plano sin el origen, la forma

[−y/(x²+y²)]dx+[x/(x²+y²)]dy

tiene derivadas cruzadas iguales, pero su potencial sería el ángulo polar, que no puede definirse como una función real continua y monovaluada en todo el dominio perforado. Al recorrer una vuelta alrededor del origen, el ángulo cambia en 2π.

10 · Exploración interactiva

Laboratorio: potencial, gradiente y tangente

Elegí un punto para la ecuación exacta (2x+y)dx+(x+2y)dy=0. La curva dorada es el nivel F=x²+xy+y²=C; el vector coral es el gradiente y el violeta es una dirección tangente.

EXPERIMENTO 17

Curvas de nivel de un potencial

x² + xy + y² = 7,00

Los resultados numéricos aparecen debajo.
ConstanteC = 7,00
M(x₀,y₀)M = 4,00
N(x₀,y₀)N = 5,00
Pendientey′ = −0,80

El gradiente (4,00; 5,00) es perpendicular a la curva de nivel.

Las curvas tenues son otros niveles del mismo potencial. En el origen, gradiente y dirección tangente se anulan.

Preguntas para explorar

  1. Mové el punto y comprobá visualmente que gradiente y tangente permanecen perpendiculares.
  2. Buscá un punto donde N=0. ¿Qué tipo de tangente tiene la curva?
  3. Elegí el origen. ¿Por qué el nivel C=0 deja de ser una elipse?
  4. Cambiá solo la escala. ¿Se modifica la constante o únicamente la visualización?

11 · Representación en código

Comprobar exactitud mediante derivadas conocidas

const derivadaMy = (x, y) => 1;
const derivadaNx = (x, y) => 1;

function comprobarExactitud(my, nx, puntos, tolerancia = 1e-12) {
  return puntos.every(([x, y]) =>
    Math.abs(my(x, y) - nx(x, y)) <= tolerancia
  );
}

const puntos = [[0, 0], [1, 2], [-2, 3]];
console.log(comprobarExactitud(derivadaMy, derivadaNx, puntos));

Representar el potencial y sus derivadas

const potencial = (x, y) => x ** 2 + x * y + y ** 2;
const M = (x, y) => 2 * x + y;
const N = (x, y) => x + 2 * y;

function datosGeometricos(x, y) {
  const gradiente = [M(x, y), N(x, y)];
  const tangente = [N(x, y), -M(x, y)];
  const productoPunto = gradiente[0] * tangente[0]
    + gradiente[1] * tangente[1];
  return { nivel: potencial(x, y), gradiente, tangente, productoPunto };
}

console.log(datosGeometricos(1, 2));

Parametrizar una curva de nivel

function puntoDelNivel(constante, angulo) {
  if (constante < 0) throw new RangeError("Este potencial exige C ≥ 0");
  const u = Math.sqrt(constante / 1.5) * Math.cos(angulo);
  const v = Math.sqrt(constante / 0.5) * Math.sin(angulo);
  return {
    x: (u + v) / Math.sqrt(2),
    y: (u - v) / Math.sqrt(2)
  };
}

const punto = puntoDelNivel(7, Math.PI / 3);
const residuo = punto.x ** 2 + punto.x * punto.y + punto.y ** 2 - 7;

console.log(punto);
console.log(Math.abs(residuo) < 1e-12);

12 · Comprensión

Confusiones frecuentes

«Exacta significa que tiene solución exacta»

Aquí «exacta» significa que Mdx+Ndy es la diferencial total de una función potencial.

«Basta comprobar My=Nx en cualquier dominio»

La suficiencia global requiere condiciones sobre la región, como ser simplemente conexa.

«Al integrar M, agrego una constante numérica»

Debe agregarse g(y), porque cualquier función de y actúa como constante respecto de x.

«La solución debe despejar y»

La forma natural es F(x,y)=C. Despejar puede dividir la curva en ramas y excluir tangentes verticales.

«Si no es exacta, no puede resolverse»

Puede admitir otro método o un factor integrante. Solo falla el método directo de potencial.

«Cualquier factor no nulo conserva la exactitud»

Multiplicar por una función arbitraria puede destruirla o no crearla. El factor integrante debe satisfacer una condición específica.

13 · Práctica guiada

Ejercicios de exactitud

Ejercicio 1: comprobar

Determiná si (3x²+2y)dx+(2x+4y³)dy=0 es exacta.

Ver solución

My=2 y Nx=2. Es exacta en todo ℝ².

Ejercicio 2: reconstruir el potencial

Resolvé la ecuación del ejercicio anterior.

Ver solución razonada

Integrando M respecto de x: F=x³+2xy+g(y). Como Fy=2x+g′(y)=2x+4y³, resulta g=y⁴. La solución es x³+2xy+y⁴=C.

Ejercicio 3: aplicar un punto

Para (2x+y)dx+(x+2y)dy=0, encontrá la curva que pasa por (2,−1).

Ver solución

C=2²+2(−1)+(−1)²=3. La curva es x²+xy+y²=3.

Ejercicio 4: detectar que no es exacta

Analizá (y+2x)dx+x²dy=0.

Ver solución

My=1 y Nx=2x. No coinciden en una región abierta, por lo que no es exacta tal como está.

Ejercicio 5: pendiente implícita

Para F=x²+xy+y², obtené y′ y verificá que coincide con −M/N.

Ver solución

Derivando F=C: 2x+y+(x+2y)y′=0. Entonces y′=−(2x+y)/(x+2y)=−M/N.

Ejercicio 6: dominio perforado

¿Por qué la igualdad de derivadas cruzadas fuera del origen no basta para que el diferencial angular tenga potencial global?

Ver solución

El dominio ℝ²\{0} tiene un agujero. La integral alrededor de una vuelta cerrada es 2π, no cero, y el ángulo cambia de rama; solo existen potenciales locales.

14 · Síntesis

Ideas para recordar

  • Una ecuación exacta tiene la forma dF=Mdx+Ndy=0.
  • Sus soluciones son curvas de nivel F(x,y)=C.
  • La condición cruzada es My=Nx.
  • La suficiencia global del criterio depende del dominio.
  • Al integrar M respecto de x, debe añadirse una función g(y).
  • También puede comenzarse integrando N respecto de y.
  • El gradiente del potencial es perpendicular a las curvas solución.
  • Una solución implícita conserva ramas y tangentes verticales.
  • Algunas formas no exactas pueden volverse exactas mediante un factor integrante.

En el próximo tema aplicaremos ecuaciones lineales de primer orden a balances de mezcla y circuitos eléctricos RC.