電磁場のディラック方程式(maxima ソース)
電磁場のディラック方程式の計算に使った maxima (wxMaxima) のソースです。コメント等ありませんが,だいたい見たらわかると思います。gt 〜 gz がγ行列です。 電流項になにか計算ミスがあるので L と L2 のふたつ行列を定義しましたが,本当は同じはずです。
code:maxima
g0:;
gt: matrix(
);
gx: matrix(
);
gy: matrix(
);
gz: matrix(
);
D:(Dt*gt-Dx*gx-Dy*gy-Dz*gz);
F: (Ex* gt.gx+ Ey* gt.gy + Ez*gt.gz) + (Bz*gx.gy+By*gz.gx+Bx*gy.gz) ;
J:(Jt*gt-Jx*gx-Jy*gy-Jz*gz);
K: matrix(
);
realpart(expand(K.D.F.L));
imagpart(expand(K.D.F.L));
realpart(expand(K.J.L));
imagpart(expand(K.J.L));
/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$
==============