请高手帮我编程求解指数矩阵（精细积分求解）：
指数矩阵满足如下关系：T=exp(B1)=[exp(B1/m)]^m=[I+(B1/m)+(B1/m)^2/2*1+(B1/m)^3/3*2*1……]^m=[I+T^(0)]^m=[I+T^(1)]^(2^N-1)，有递推式：T^(i+1)=2*T^(i)+T^(i)*T^(i)，其中m=2^N，N=5，根据以上条件编程得到exp(B1)具体的值。
以下是我编写的程序，但是得不到结果，不知道什么原因造成的，请高手调试……十分谢谢！
function Tn=juzhen(Ta,N)
syms k w
A=[0 -i*k 0 0 0 0 0 -2/0.7
    -0.3*i*k 0 -1 0 0 0 1 0
    0 1 0 -1 0 0 0 0
    0 0 0.3*(i*k)^2 0 1 0 0 0
    0 0 0 -2*0.7*(i*k)^2 0 -12*(1/0.02)^2 0 0
    0 0 -0.91*(0.02^2/12)*(i*k)^4+7840*0.02*w^2*0.91/(2.1*10^11*0.02) 0 0.3*(i*k)^2*0.02^2/12 0 -1 0
    0 -7840*0.02*w^2*0.91/(2.1*10^11*0.02) 0 0 0 1 0 i*k
    0.91*(i*k)^2+7840*0.02*w^2*0.91/(2.1*10^11*0.02) 0 0 0 0 0 0.3*i*k 0]
F=[0 0 0;0 0 0;0 0 0;0 0 0;0 0 0;1/2*pi 1/2*pi 1/4*pi;0 0 0;0 0 0]
A1=[0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0]
A2=[0 0 0;1 0 0;0 2 0]
B=[A, F;A1, A2]
B1=pi*B
N=5,m=2^N
I=eye(size(B1))
Ta=(B1/m)+((B1/m)^2)/2+((B1/m)^3)/(3*2)+((B1/m)^4)/(4*3*2)
for j=0:1:N-1
    Ts=2*Ta+Ta*Ta
    Ta=Ts
end
Tn=Ts
T=I+Tn
还有一个问题：
series(Dirac(a)/2*Pi,a=0,3)=Dirac(0)/2*Pi+a*Dirac(1,0)/2*Pi+a^2*Dirac(2,0)/4*Pi+0(a^3)   请问怎样根据a来确定Dirac(0)、Dirac(1,0)、Dirac(2,0)的取值？