softwarepackages.Mostpackageswillalsoreturntheeigenvectors,if
requested.UsingOx:
#include
main()
{
decl ECN437
ComputingEigenvaluesUsingSoftware
Packages
Hereareafewexamplesshowinghoweigenvaluesarecomputedinvarious a,eigs;
a = <1, 2; 3, 4>;
// Compute the eigenvalues and store them in eigs.
The & is
// needed so that Ox can use eigs to return the
eigenvalues.
eigen( &eigs a, );
print( eigs );
}
UsingGauss:
let x[2,2] = {1 2,3 4};
eig(x);
UsingMaple:
# This line activates the linear algebra package. It
produces two
# messages about norm and trace being redefined, which is
OK.
with(linalg):
#
#
#
Create a matrix
a:=matrix( [ [1,2],[3,4] ] );
#
#
#
Calculate its eigenvalues
b:=eigenvals(a);
#
#
#
Exit from maple
quit;