title: Teorema de Malus-Dupin
text: # Teorema de Malus-Dupin
  
  El teorema de Malus-Dupin es uno de los teoremas fundamentales de la óptica geométrica que la relaciona con la óptica ondulatoria.
  
  # Teorema de Malus-Dupin
  ## Enunciado
  
  Si sobre cada rayo emitido por un foco recorremos caminos ópticos iguales, entonces los puntos que los delimitan forman una superficie normal a todos los rayos. Denominamos a dicha superficie frente de ondas. Coincide con el frente de onda dado por la teoría oscilatoria. Al deducirse del principio de Fermat, es válido a pesar del número de reflexiones o refracciones que pueda sufrir el rayo antes de llegar a su destino.
  
  # Teorema de Malus-Dupin
  ## Demostración
  
  La demostración se realiza a partir del principio de Fermat.
  Tomamos dos trayectorias distintas separadas infinitesimalmente, $[AB]$ y $[AB']$, donde $A$ es el foco y $B$ y $B'$ son los puntos de llegada separados por caminos ópticos iguales. Entonces definimos los respectivos caminos ópticos como:
  $L_{AB}=\int _{A}^{B}n({\textbf {r}})ds$
  $L_{AB'}=\int _{A}^{B'}n({\textbf {r}}')ds'$
  Siendo ${\textbf {r}}$ y ${\textbf {r}}'$ los respectivos vectores de posición, $ds$ y $ds'$ los respectivos diferenciales de espacio y $n({\textbf {r}})$ el índice de refracción.
  Admitimos que el índice de refracción es derivable.
  Emplearemos:
  Desarrollo en serie de Taylor de primer orden: $f({\textbf {r}})=f({\textbf {r}}_{0})+{\vec {\nabla }}f({\textbf {r}}_{0})\cdot ({\textbf {r}}-{\textbf {r}}_{0})$
  Ecuación de la trayectoria de un rayo luminoso (deducida a partir del principio de Fermat):${d \over ds}(n{\textbf {u}})={\vec {\nabla }}n$ de lo que $(n{\textbf {u}})={\vec {\nabla }}nds$.
  La relación ${\textbf {r}}'={\textbf {r}}+\varepsilon {\textbf {b}}$ de la cual ${\textbf {dr}}'={\textbf {dr}}+\varepsilon {\textbf {db}}$.
  Admitimos que el índice de refracción admite un desarrollo en serie de Taylor de orden 1. Entonces se obtiene que
  $n({\textbf {r}}')=n({\textbf {r}})+{\vec {\nabla }}n({\textbf {r}})\cdot ({\textbf {r}}'-{\textbf {r}})=n({\textbf {r}})+{\vec {\nabla }}n({\textbf {r}})\cdot \varepsilon {\textbf {b}}$.
  Por otro haremos lo mismo con el módulo del vector posición:
  $\mid {\textbf {r}}'\mid =\mid {\textbf {r}}\mid +{\vec {\nabla }}\mid {\textbf {r}}\mid \cdot ({\textbf {r}}'-{\textbf {r}})=\mid {\textbf {r}}\mid +{{\partial } \over {\partial r}}r{\hat {u}}_{r}\cdot ({\textbf {r}}'-{\textbf {r}})=\mid {\textbf {r}}\mid +{\textbf {u}}\cdot ({\textbf {r}}'-{\textbf {r}})=\mid {\textbf {r}}\mid +{\textbf {u}}\cdot {\textbf {r}}'-{\textbf {u}}\cdot {\textbf {r}}=\mid {\textbf {r}}\mid +{\textbf {u}}\cdot {\textbf {r}}'-\mid {\textbf {r}}\mid ={\textbf {u}}\cdot {\textbf {r}}'$.
  De modo que:
  $ds'=\mid {\textbf {dr'}}\mid ={\textbf {u}}\cdot {\textbf {dr}}'={\textbf {u}}\cdot {\textbf {dr}}+{\textbf {u}}\cdot \varepsilon {\textbf {db}}=ds+\varepsilon {\textbf {u}}\cdot {\textbf {db}}$.
  Se reemplaza en el camino óptico:
  $L_{B'}=\int _{A}^{B'}n({\textbf {r}}')ds'=\int _{A}^{B'}[n({\textbf {r}})+{\vec {\nabla }}n\cdot \varepsilon {\textbf {b}}](ds+\varepsilon {\textbf {u}}\cdot {\textbf {db}})$
  $L_{B'}=\int _{A}^{B'}n({\textbf {r}})ds+\int _{A}^{B'}{\vec {\nabla }}n\cdot \varepsilon {\textbf {db}}+\int _{A}^{B'}n({\textbf {r}})\varepsilon {\textbf {u}}\cdot {\textbf {db}}+\int _{A}^{B'}{\vec {\nabla }}n\cdot \varepsilon {\textbf {b}}\varepsilon {\textbf {u}}\cdot {\textbf {db}}$
  Eliminamos los elemento de orden de $\varepsilon ^{2}$:
  $L_{B'}=\int _{A}^{B'}n({\textbf {r}})ds+\int _{A}^{B'}{\vec {\nabla }}n\cdot \varepsilon {\textbf {db}}+\int _{A}^{B'}n({\textbf {r}})\varepsilon {\textbf {u}}\cdot {\textbf {db}}$
  Calculamos $\Delta L$:
  $\Delta L=\int _{A}^{B'}n({\textbf {r}})ds+\int _{A}^{B'}{\vec {\nabla }}n\cdot \varepsilon {\textbf {b}}ds+\int _{A}^{B'}n({\textbf {r}})\varepsilon {\textbf {u}}\cdot {\textbf {db}}-\int _{A}^{B}n({\textbf {r}})ds$
  Factorizamos los términos por potencias de $\varepsilon$:
  $\Delta L=\int _{A}^{B'}n({\textbf {r}})ds-\int _{A}^{B}n({\textbf {r}})ds+\varepsilon \left[\int _{A}^{B'}{\vec {\nabla }}n\cdot {\textbf {b}}ds+\int _{A}^{B'}n({\textbf {r}}){\textbf {u}}\cdot {\textbf {db}}\right]$
  Por la ecuación de las trayectorias tenemos que $d(n{\textbf {u}}\cdot {\textbf {b}})=d(n{\textbf {u}})\cdot {\textbf {b}}+n{\textbf {u}}\cdot {\textbf {db}}={\vec {\nabla }}n\cdot {\textbf {b}}ds+n{\textbf {u}}\cdot {\textbf {db}}$ de modo que:
  $\Delta L=\int _{A}^{B'}n({\textbf {r}})ds-\int _{A}^{B}n({\textbf {r}})ds+\varepsilon \left[\int _{A}^{B'}d(n{\textbf {b}}\cdot {\textbf {u}})\right]$
  Si suponemos que hemos escogido $B\equiv B'$ entonces:
  $\Delta L=\int _{A}^{B'}n({\textbf {r}})ds-\int _{A}^{B'}n({\textbf {r}})ds+\varepsilon \left[\int _{A}^{B'}d(n{\textbf {b}}\cdot {\textbf {u}})\right]=0$
  $\Delta L=+\varepsilon \left[\int _{A}^{B'}d(n{\textbf {b}}\cdot {\textbf {u}})\right]=0$
  Por lo que:
  $\int _{A}^{B'}d(n{\textbf {b}}\cdot {\textbf {u}})=0$
  $n({\textbf {r}}_{B'}){\textbf {b}}({\textbf {r}}_{B'})\cdot {\textbf {u}}({\textbf {r}}_{B'})-n({\textbf {r}}_{A}){\textbf {b}}({\textbf {r}}_{A})\cdot {\textbf {u}}({\textbf {r}}_{A})=0$
  Como el punto A es el foco la separación es siempre nula, en consecuencia:
  $n({\textbf {r}}_{B'}){\textbf {b}}({\textbf {r}}_{B'})\cdot {\textbf {u}}({\textbf {r}}_{B'})=0$
  ${\textbf {b}}({\textbf {r}}_{B'})\cdot {\textbf {u}}({\textbf {r}}_{B'})=0$
  ${\textbf {b}}({\textbf {r}}_{B'})\perp {\textbf {u}}({\textbf {r}}_{B'})$
  Como B' es un punto arbitrario se tiene que ${\textbf {b}}({\textbf {r}})\perp {\textbf {u}}({\textbf {r}})$.
  Teníamos que ${\textbf {r}}'={\textbf {r}}+\varepsilon {\textbf {b}}$, por lo que $\varepsilon {\textbf {b}}({\textbf {r}})$ une los puntos de la superficie. De lo que la superficie formada es ortogonal a cada rayo. Podemos justificar que los puntos forman una superficie por continuidad.
  * Datos: Q3527117