A=matrix([[0,0,1,0,0],[0,0,0,1,0],[0,0,0,0,1],[0,0,0,0,0],[0,0,0,0,0]])
print A
A.jordan_form(transformation=True)
B=matrix([[1,-1,0,-1],[0,2,0,1],[-2,1,-1,1],[2,-1,2,0]])
print B
B.jordan_form(transformation=True)
u=var('u')
t=var('t')
Au=var('Au')
def f(x):
    return Au
integral(f(x),t)

Kreyszig-AppendixB-5-U (last edited 2010-12-18 11:01:39 by sri valli)