Description
While the desolve_system command can fully solve a given linear system of differential equations, it is still useful to know how to compute a general solution by finding the eigenvalues and eigenvectors of the corresponding matrix. Consider the system
(1)The corresponding matrix to this system is
(2)We may use Sage to find both eigenvalues and an eigenvector for each eigenvalue:
Sage Cell
Notice that the .eigenvectors_right() method gives us the eigenvalue corresponding to each eigenvector, so the .eigenvalues() method is redundant in this case, bit it's still useful to know. From the output, we see we have two eigenvalues: $\lambda_1 = -3$ with eigenvector $\left<1, -1/6\right>$ and $\lambda_2 = 2$, with eigenvector $\left<1, -1\right>$. So, letting $c_1$ and $c_2$ be our arbitrary constants, the general solution to our system is
(3)Code
A = matrix(QQ, [[-4, -6], [1, 3]])
print(A.eigenvalues())
print(A.eigenvectors_right())
Options
none
Tags
Primary Tags:
Secondary Tags:
A list of possible tags can be found at The WeBWorK Open Problem Library. For linear algebra tags see the Curated Courses Project.
Related Cells
Any related cells go here. Provide a link to the page containing the information about the cell.
Attribute
Permalink:
Author:
Date: 13 May 2020 04:31
Submitted by: Zane Corbiere
