View difference between Paste ID: nVYArB8n and YB5pGeSs

-def eqS((A, B, C, D), (An, Bn, Cn, Dn)):
+def analitInitial((A, B, C, D), x0, T):
-    A = mt(A); B = mt(B); C = mt(C); D = mt(D)
+    T = arange(T[0], T[1], T[2])
-    An = mt(An); Bn = mt(Bn); Cn = mt(Cn); Dn = mt(Dn)
+    x0 = array(x0)
-    b = B[:, 0]
+
-    bn = Bn[:, 0]
+    J, V = eig(A)
-    n = len(A)
+    J = mt(J * eye(3, 3))
-    Sy = b
+    S = mt(V)
-    Sny = bn
+    y = ([], [])
-    for i in range(1, n):
+
-        Apow = copy(A)
+    for i in range(len(T)):
-        Anpow = copy(An)
+        y[0].append(T[i])
-        if not i == 1:
+        yy = C * S * mt(expm(J * T[i])) * mt(inv(V)) * mt(x0).transpose()
-            for j in range(i-1):
+        y[1].append(array(yy[0])[0][0])
-                Apow = dot(Apow, A)
+    plot(y[0], y[1])
-                Anpow = dot(Anpow, An)
+    title('Analytical construction')
-        Sy = hstack((Sy, dot(Apow, b)))
+    xlabel('t, s.')
-        Sny = hstack((Sny, dot(Anpow, bn)))
+    ylabel('y(t)')
-    return dot(Sy, inv(Sny))
+    grid(color='g', linestyle='--')
     show()