Definitionen der Impulsoperatoren
px := psi -> hbar/I*diff(psi,x);
py := psi -> hbar/I*diff(psi,y);
pz := psi -> hbar/I*diff(psi,z);
Ortsoperatoren
XX := psi -> x*psi;
YY := psi -> y*psi;
ZZ := psi -> z*psi;
Drehimpulsoperatoren
lx := psi->YY(pz(psi))-ZZ(py(psi));
ly := psi->ZZ(px(psi))-XX(pz(psi));
lz := psi->XX(py(psi))-YY(px(psi));
lx(psi(x,y,z));ly(psi(x,y,z));lz(psi(x,y,z));
Drehimpulsquadrat
l2 := psi->lx(lx(psi))+ly(ly(psi))+lz(lz(psi));
collect(simplify(l2(psi(x,y,z))),[x,y,z],distributed);
Kommutatoren
Kom := (a,b,psi)->a(b(psi))-b(a(psi));
simplify(Kom(lx,ly,psi(x,y,z)));simplify(Kom(ly,lz,psi(x,y,z)));simplify(Kom(lz,lx,psi(x,y,z)));
simplify(Kom(l2,lx,psi(x,y,z)));simplify(Kom(l2,ly,psi(x,y,z)));simplify(Kom(l2,lz,psi(x,y,z)));
Zentralsymmetrisch
RR := psi-> (x^2+y^2+z^2)^(n/2)*psi;
simplify(Kom(RR,lx,psi(x,y,z)));
simplify(Kom(RR,ly,psi(x,y,z)));
simplify(Kom(RR,lz,psi(x,y,z)));
simplify(Kom(RR,l2,psi(x,y,z)));