ARCHIVOS TEÓRICOS DEL MUSEO MATEMÁTICO.
ÁREA- ANÁLISIS-
TEMA-CÁLCULO DE VARIACIONES.ENFOQUE GLOBAL.
////// Problemas variacionales ///////////
Autor- Leonard Echagüe-Mate UBA Museum- lechague@dm.uba.ar
Cicloide como minimal de tiempo de caída.(Brachistócrona).
Planteo por cálculo de variaciones:
Por Euler-Lagrange simplificada
f-diff(f, )* =a
Quedando para resolver:
> restart;
> z:=diff(y(x),x);
> ode5 := a=y(x)+y(x)*z^2;
>
> ans5 := dsolve(ode5);
> with(DEtools): odeadvisor(ode5);
> z:=diff(y(x),x);
> ode6 := 10=y(x)+y(x)*z^2;
> ans6 := dsolve(ode6);
Copiando, pegando y adaptando constantes....
> plot({sqrt(-y^2+10*y)-5*arcsin(-1+1/5*y),-10/2*Pi-sqrt(-y^2+10*y)+5*arcsin(-1+1/5*y)},y=0..10);
>
> y1:=0;sqrt(-y1^2+10*y1)-5*arcsin(-1+1/5*y1);evalf(%);y1:=10;sqrt(-y1^2+10*y1)-5*arcsin(-1+1/5*y1);evalf(%);
Verificación gráfica.....
> with(plots):
> p1:=plot({sqrt(-y^2+10*y)-5*arcsin(-1+1/5*y),-10/2*Pi-sqrt(-y^2+10*y)+5*arcsin(-1+1/5*y)},y=0..10):
> p2:=plot([-5*cos(t)+5,5*t-5*sin(t)-15/2*Pi,t=0..2*Pi],color=blue,style=point):
> display([p1,p2]);
>