Wow. Cool problem. Putzer's algorithm is not that old, so it is not that well known.
But, having researched matrix exponentials, I think I may be able to help a wee bit.
By Putzer's algorithm,
eAt=p1(t)M0+p2(t)M1
Where
M0=I, M1=[a−λ1cbd−λ1]
p1(t)=eλ1t, p2(t)=λ1−λ21(eλ1t−eλ2t)
Of course, I assume you know that the 2 lambdas represent the eigenvalues of matrix A.
The characteristic polynomial of said matrix is
λ2+4. Which means the eigenvalues are
λ1=2i, λ2=−2i
Now, proceeding. I am not going to go into every derivation and proof of why and/or how this works. You can probably find all that somewhere if you need it.
Therefore, from all of the above:
eAt=[eλ1t+λ1−λ2a−λ1(eλ1t−eλ2t)λ1−λ2c(eλ1t−eλ2t)λ1−λ2b(eλ1t−eλ2t)eλ1t+λ1−λ2d−λ1(eλ1t−eλ2t)]
Plug in a,b,c,d and the eigenvalues.
The matrix exponential is used to solve DE's. So, I suppose this is all that is needed because I see no initial values.