From: Jaime V. <vi...@fe...> - 2008-02-07 12:50:40
|
On Thu, 2008-02-07 at 12:31 +0100, Javier Arantegui wrote: > Realmente el problema no es resolver la ecuación diferencial, ... > El problema es que cuando utilizo una función para poner el seno, como > Q1(t), Maxima no calcula la transformada de Laplace. Ah! entonces ahora si entendí tu pregunta y te puedo responder. Observa lo siguiente: (%i2) eq:A*'diff(H(t),t)=Q1(t)-H(t)/R$ (%i3) Q1(t):=sin(t)$ (%i4) atvalue(H(t), t=0,0)$ (%i5) desolve(eq,H(t)); (%o5) H(t) = ?%ilt(?%laplace(sin(t),t,lvar)*R/(lvar*A*R+1),lvar,t) (%i6) desolve(''eq,H(t)); (%o6) H(t) = A*R^2*%e^-(t/(A*R))/(A^2*R^2+1)-cos(t)*A*R^2/(A^2*R^2+1) +sin(t)*R/(A^2*R^2+1) Como Q1 era desconocida en el momento en que definiste eq, cuando le pasas eq a desolve le estás pasando una función Q1 que desolve desconoce. En el momento de mostrar los resultados, Maxima, que ya sabe lo que es Q1, te muestra el valor de Q1. En %i6 he usado el operador '' para forzar a que las funciones conocidas sean substituidas en eq, pasando sin en vez de Q1 como entrada para desolve. Espero ahora si haber conseguido ayudarte. Jaime |