Quantcast
Viewing all articles
Browse latest Browse all 1087

Why doesn't solve give any indication that there can be a non-trivial solution to this system of two linear equations?

I created the following worksheet to illustrate my question.

Image may be NSFW.
Clik here to view.
restart

We have two equations. Image may be NSFW.
Clik here to view.
w
 and Image may be NSFW.
Clik here to view.
w__0
 are parameters and we wish to solve for Image may be NSFW.
Clik here to view.
A__1
 and Image may be NSFW.
Clik here to view.
A__2
.

Image may be NSFW.
Clik here to view.
eq1 := A__1*(w^2-3*w__0^2)+A__2*w__0^2 = 0
 = Image may be NSFW.
Clik here to view.
A__1*(w^2-3*w__0^2)+A__2*w__0^2 = 0
Image may be NSFW.
Clik here to view.
NULL

Image may be NSFW.
Clik here to view.
eq2 := A__1*w__0^2+A__2*(w^2-w__0^2) = 0
 = Image may be NSFW.
Clik here to view.
A__1*w__0^2+A__2*(w^2-w__0^2) = 0
Image may be NSFW.
Clik here to view.
NULL

Image may be NSFW.
Clik here to view.
solve({eq1, eq2}, {A__1, A__2})
 = Image may be NSFW.
Clik here to view.
{A__1 = 0, A__2 = 0}
Image may be NSFW.
Clik here to view.
NULL

 

Maple says the only solution is the trivial solution.

If we check the determinant of the matrix of the system we see it can be zero for certain values of Image may be NSFW.
Clik here to view.
w
 given Image may be NSFW.
Clik here to view.
w__0
.

Image may be NSFW.
Clik here to view.
M := Matrix([[w^2-3*w__0^2, w__0^2], [w__0^2, w^2-w__0^2]])
 = Image may be NSFW.
Clik here to view.
Matrix(%id = 36893488152149895276)
Image may be NSFW.
Clik here to view.
NULL

Image may be NSFW.
Clik here to view.
d := LinearAlgebra:-Determinant(M)
 = Image may be NSFW.
Clik here to view.
w^4-4*w^2*w__0^2+2*w__0^4
Image may be NSFW.
Clik here to view.
NULL

Image may be NSFW.
Clik here to view.
solutions := solve(d = 0, w)

Image may be NSFW.
Clik here to view.
(2+2^(1/2))^(1/2)*w__0, -(2+2^(1/2))^(1/2)*w__0, (2-2^(1/2))^(1/2)*w__0, -(2-2^(1/2))^(1/2)*w__0

(1)

If Image may be NSFW.
Clik here to view.
w
 is one of these values then the system of equations is singular and has non-zero solutions.Image may be NSFW.
Clik here to view.
``

 

For example, subbing the first value above into the equations and solving Maple gives us non-trivial solutions.

 

Image may be NSFW.
Clik here to view.
solve({subs(w = solutions[1], eq1), subs(w = solutions[1], eq2)}, {A__1, A__2})

Image may be NSFW.
Clik here to view.
{A__1 = -(1+2^(1/2))*A__2, A__2 = A__2}

(2)

Image may be NSFW.
Clik here to view.
``

Why didn't Maple give us any indication that there could be non-zero solutions in the first call to Image may be NSFW.
Clik here to view.
solve
?


For context, the system of equations comes from a calculation involving coupled oscillators.

I had a system of differential equations, guessed at a solution, plugged it in and got the equations shown in the worksheet. The guess isn't a correct solution in general, but it is a solution if w is one of the values computed in the worksheet (the values of the variable "solutions").

Download solve-nonzero-solutions.mw


Viewing all articles
Browse latest Browse all 1087

Trending Articles