Después de leer http://math.tut.fi/~piche/pde/pde.pdf , no sé cómo calcular el vector propio
Cómo encontrar el vector propio de un segundo orden de la ecuación diferencial?
¿por qué algunas personas usan el pecado como autovector? es sólo el pecado puede ser vector propio?
El problema es que eigenfunction de expansión, el primer paso es encontrar autovalor y autovector, pero no sé cómo calcular el vector propio de la ecuación diferencial
por ejemplo Arce código
x*diff(f(x), x$2) + 2*x*diff(f(x),x) + f(x) = 0
x*diff(f(x), x$2) + 2*x*diff(f(x),x) + x = 0
Actualizado
sol := dsolve(t*diff(phi(x),x$2)-x*diff(phi(x),x)+n*phi(x),phi(x));
phi := unapply(rhs(sol), x);
BC := [phi(0)=0,phi(1)=0];
with(linalg):
Ccoef := genmatrix(BC, [_C1,_C2]);
CharEqn := det(Ccoef) = 0;
restart;
sol := dsolve(t*diff(phi(x,t,n),x$2)-x*diff(phi(x,t,n),x)+n*phi(x,t,n),phi(x,t,n));
phi := unapply(rhs(sol), x);
BC := [phi(0,0,0)=0,phi(1,1,1)=0];
with(linalg):
Ccoef := genmatrix(BC, [_C1,_C2]);
CharEqn := det(Ccoef) = 0;
**lo siento, solo los domingos momento oportuno para leer este archivo, me parece que el pecado de la función que viene de la etapa de cálculo de la ecuación característica uso del archivo pdf método para calcular la anterior ecuación diferencial para eignvector,
esta ecuación es Ermitaño después de probado, la ecuación característica es cero, es decir no autovector supongo que este cálculo de arce código tiene algo de malo
cómo calcular esto?**
Actualizado el 2 de
Originalmente espero encontrar Ermitaño H(x) y, a continuación, utilizar la suma(H*z^m/m!, m=0..infinity) encontrar un A*exp(B) donde B es en términos de z y t y es sólo una simple fórmula ahora siguiendo los pasos, supongo que la H es la solución de la función de green acerca de la expansión
se vuelven más compicated para H(x), y me parece que hay un D[2] pero no sabes de dónde vienen. entonces no sabes lo que paso es H(x), me imagino vterm o vv
sol := dsolve(t*diff(phi(x),x$2)-x*diff(phi(x),x)+n*phi(x),phi(x));
phi := unapply(rhs(sol),x);
odetest(sol,ode);
eq1:=limit(rhs(sol),x=0,right)=0;
eq2:=eval(rhs(sol),x=1)=0;
Ccoef := LinearAlgebra:-GenerateMatrix([eq1,eq2],[_C1,_C2]);
CharEqn:=LinearAlgebra:-Determinant(%[1])=0;
solve(CharEqn,t);
step1 := map(xi->simplify(subs(t=RootOf(KummerM(1/2-(1/2)*n, 3/2, 1/(2*_Z))),xi)),Ccoef);
with(linalg):
NN := nullspace(step1);
subs(_C1=NN[1][1],_C2=NN[1][2],t=RootOf(KummerM(1/2-(1/2)*n, 3/2, 1/(2*_Z))),phi(x));
phi := (n,t,x) -> KummerM(1/2-(1/2)*n, 3/2, (1/2)*x^2/RootOf(KummerM(1/2-(1/2)*n, 3/2, 1/(2*_Z))))*x;
assume(j,posint):
interface(showassumed=0):
Gterm := unapply(-phi(n,t,x)*phi(n,t,x)*exp(-lambda(j)*t)/int(phi(n,t,x)^2,x=0..1),(j,n,x,y,t)):
G:=Sum(Gterm(j,n,x,y,t),j=1..infinity);
vterm := int(D[2](Gterm)(n,1,x,t-tau),tau=0..t);
vv := sum(Sum(op(n,vterm),j=1..infinity),n=1..2);