It doesn't take much to start building interpolants. import numpy as np import loewner_framework as lf import loewner_framework.linear_daes as ld Lambda = np.diag([0. ...