All published worksheets from http://sagenb.org
from sympy import * from sympy.matrices import Matrix from numpy import * a,b,x=symbols("abx")
Mat_A=Matrix(( (a,b,0,0), (x,a,b,0), (x,-2*x,a,b), (-x**2,-x,x,a), )) Mat_A
h=symbols("h") eq_1=(Mat_A-eye(4)*h).det() y=symbols("y") eq_2=eq_1.subs(h,a+y) eq_2
eq_2.expand()
lambda=a+sqrt(abs(x)*abs(b)*sqrt(-1*(b+1))) lambda