bx=zeros(M,N,P);
by=zeros(M,N,P);
bz=zeros(M,N,P);
% wx= C cos y + A sin z
% wy= A cos z + B sin x
% wz= B cos x + C sin y
i=sqrt(-1);
A=0.4; B=0.3; C=0.5;
bx(1,2,1)=C/2;
bx(1,N,1)=C/2;
bx(1,1,2)=A/2/i;
bx(1,1,P)=-A/2/i;
by(1,1,2)=A/2;
by(1,1,P)=A/2;
by(2,1,1)=B/2/i;
by(M,1,1)=-B/2/i;
bz(2,1,1)=B/2;
bz(M,1,1)=B/2;
bz(1,2,1)=C/2/i;
bz(1,N,1)=-C/2/i;

A=0.5; B=0.4; C=0.3;
bx(1,3,1)=C/2;
bx(1,N-1,1)=C/2;
bx(1,1,3)=A/2/i;
bx(1,1,P-1)=-A/2/i;
by(1,1,3)=A/2;
by(1,1,P-1)=A/2;
by(3,1,1)=B/2/i;
by(M-1,1,1)=-B/2/i;
bz(3,1,1)=B/2;
bz(M-1,1,1)=B/2;
bz(1,3,1)=C/2/i;
bz(1,N-1,1)=-C/2/i;

t=0;
deltat=0.005;
g=1;
nu=0.06;
% OTHER INITIALIZATION CAN GO HERE

[kx,ky,kz,ksquare]=makek(M,N,P);