Ecuaciones en maple

Solo disponible en BuenasTareas
  • Páginas : 13 (3049 palabras )
  • Descarga(s) : 0
  • Publicado : 25 de septiembre de 2010
Leer documento completo
Vista previa del texto
Universidad Diego Portales Facultad de Ciencias de la Ingeniería Enero 2003

Introducción a las Ecuaciones Diferenciales en Maple
Jonathan Makuc (Jonathan.Makuc@al.udp.cl)
Esta guía tiene como objetivo mostrar algunas aplicaciones básicas del programa a las Ecuaciones Diferenciales. Los tópicos contenidos son: • • • • • • Como ingresar una ecuación diferencial Obtener el tipo de Ec. DiffResolución de Ec. Diff Ordinarias (ODE – Ordinary Diferential Ecuation) Resolución básica de Ec. Diff Parciales Resolución de Ec. Diff mediante aproximaciones numéricas. Graficación de Soluciones

La librería de interés es DEtools.

¡ CARGUE ESTA LIBRERÍA ANTES DE USAR ALGUN COMANDO DE ESTA GUÍA !

Como Ingresar una Ecuación Diferencial
Para ingresar una Ec. Diff, debemos entender que elprograma no asume nada, por lo tanto si tenemos una función f que depende de x ( f(x) ) debemos explicitar esta dependencia a Maple. Ej: dx (y-1) + dy (x+2) = 0 y −1 dy =− x+2 dx

Lo primero que debemos hacer es llevarlo a la forma

Nosotros sabemos que la solución que estamos buscando es una expresión y(x), por lo tanto escribimos en Maple: > ec1:=(y(x)-1)/(x+2)=-diff(y(x),x); y( x ) − 1 ∂ ec1 := =− y( x )   ∂x    x+2  
2 && + g sen (θ ) = 0 donde θ = d θ , una forma de && Para una ecuación de segundo grado como θ l dt 2 ingresarla sería:

> ec2:=diff(theta(t), t$2)+g/l*sin(theta(t))=0;  ∂2  g sin ( θ ( t ) ) ec2 :=  2 θ( t )  + =0  ∂t  l  

Obtener el tipo de Ecuación Diferencial
Nos introducimos ahora en el la librería DEtools de Maple, por lo tanto la cargamos. >with(DEtools): Se ha utilizado dos punto en lugar de punto y coma para ahora espacio y no presentar diferencias con las variadas versiones de Maple. Si Ud. Desea ver los comandos que incluye esta librería en su versión, solo reemplace los dos puntos por punto y coma. Con el simple uso del comando ‘odeadvisor’ podemos averiguar que tipo de Ec. Diff estamos trabajando. Usando los mismo ejemploanteriores: > odeadvisor(ec1); [ _separable ] Nos entrega que la ecuación ‘ec1’, es decir, correcto. Analizando la segunda ecuación. > odeadvisor(ec2); [ [ _2nd_order , _missing_x ], [ _2nd_order , _reducible, _mu_x_y1] ] y −1 dy = − es de variable separable, lo cual es x+2 dx

Obtenemos que es de segundo orden con un término ‘x’ no presente. También es considerada como una ecuación de segundo ordenreducible. El tercer término se hace relación al factor integrante a utilizar. Este comando ha evolucionado con las versiones de Maple. Aquí (y en toda la guía) hemos presentado la salida da la versión 7.00. Veamos que ocurre con la versión 5.00: > odeadvisor(ec2); [[_2nd_order, _missing_x], [_2nd_order, _with_linear_symmetries]] Obtenemos resultados diferentes. También identifica a la ecuacióncomo de segundo grado con termino ‘x’ ausente, pero a hora encuentra ‘simetrías lineales’ en ella. Analicemos otras ecuaciones conocidas. Aquí introducimos otras formas de ingresar las ecuaciones a Maple. > ec3:=D(y)(x)=(x^2-y(x)^2)/(x*y(x)); odeadvisor(ec3); x 2 − y( x ) 2 ec3 := D( y )( x ) = x y( x ) [ [ _homogeneous, class A ], _rational, _Bernoulli ] dy x 2 − y 2 = dx xy

>ec4:=Diff(y(x),x)=(x+y(x)-3)/(x-y(x)-1); odeadvisor(ec4); ∂ x + y( x ) − 3 ec4 := y( x ) = ∂x x − y( x ) − 1

dy x + y − 3 = dx x − y − 1

[ [ _homogeneous, class C ], _rational, [ _Abel, 2nd type, class A ] ]

> ec5:=D(y)(x)+p*y(x)=q; odeadvisor(ec5);

ec5 := D( y )( x ) + p y( x ) = q [ _quadrature ]

dy + y⋅ p =q dx

> ec6:=D(y)(x)+P(x)*y(x)=Q(x); odeadvisor(ec6); ec6 := D( y )( x ) + P( x ) y( x )= Q( )
x

dy + y ⋅ P( x ) = Q ( x ) dx

[ _linear ]

Truco de teclado: hemos realizado dos comandos por grupo de ejecución (el ‘>’ que presenta Maple). Para hacer esto, presione SHIFT – ENTER al final de la línea y una nueva línea se agregará sin ejecutar el comando que antecede. Cabe mencionar que Maple conoce la forma de muchas otras Ec. Diff, como Riccati, Bernoulli, exacta, etc. Si...