function X=eAtfjs2(A,t,e)
[n,m]=size(A);
a=eig(A);
Q=eye(n);
syms s
a(1)=fix(a(1)*10^e)/10^e;
r=exp(a(1)*t);
S=r*Q;
for i=2:n;
    r1=inline(r);
    a(i)=fix(a(i)*10^e)/10^e;
    r=int(exp(a(i)*(t-s))*r1(s),s,0,t);
    Q=(A-a(i-1)*eye(n))*Q;
    S=S+r*Q;
end
X=S;