import numpy as np
from scipy.linalg import logm, expm

A = np.array([[1, 0, 0, 30],
              [0, 1, 0,  0],
              [0, 0, 1,  0],
              [0, 0, 0,  1]])
logA = logm(A)
print("log(A):")
print(logA)
print("From I to A in steps of 0.1:")
for t in range(0, 11):
  print(expm(t * 0.1 * logA))

B = np.array([[-4.52e-7, 0, -1      , 19.1],
              [ 0      , 1,  0      ,  0  ],
              [ 1      , 0, -4.52e-7, 19.1],
              [ 0      , 0,  0      ,  1  ]])
logB = logm(B)
print("log(B):")
print(logB)
print("From I to B in steps of 0.1:")
for t in range(0, 11):
  print(expm(t * 0.1 * logB))