**1.5.7 Solving Equations**

An expression like x^2+2x-7==0 represents an equation in Mathematica. You will often need to solve equations like this, to find out for what values of x they are true.

This gives the two solutions to the quadratic equation . The solutions are given as replacements for x

.
In[1]:= **Solve[x^2 + 2x - 7 == 0, x]**

Out[1]=

Here are the numerical values of the solutions.
In[2]:= **N[ % ]**

Out[2]=

You can get a list of the actual solutions for x by applying the rules generated by Solve to x using the replacement operator.
In[3]:= **x /. %**

Out[3]=

You can equally well apply the rules to any other expression involving x.
In[4]:= **x^2 + 3 x /. %%**

Out[4]=

Finding and using solutions to equations.

Solve always tries to give you explicit formulas for the solutions to equations. However, it is a basic mathematical result that, for sufficiently complicated equations, explicit algebraic formulas cannot be given. If you have an algebraic equation in one variable, and the highest power of the variable is at most four, then Mathematica can always give you formulas for the solutions. However, if the highest power is five or more, it may be mathematically impossible to give explicit algebraic formulas for all the solutions.

Mathematica can always solve algebraic equations in one variable when the highest power is less than five.
In[5]:= **Solve[x^4 - 5 x^2 - 3 == 0, x]**

Out[5]=

It can solve some equations that involve higher powers.
In[6]:= **Solve[x^6 == 1, x]**

Out[6]=

There are some equations, however, for which it is mathematically impossible to find explicit formulas for the solutions. Mathematica uses Root objects to represent the solutions in this case.
In[7]:= **Solve[2 - 4 x + x^5 == 0, x]**

Out[7]=

Even though you cannot get explicit formulas, you can still find the solutions numerically.
In[8]:= **N[ % ]**

Out[8]=

In addition to being able to solve purely algebraic equations, Mathematica can also solve some equations involving other functions.

After printing a warning, Mathematica returns one solution to this equation.
In[9]:= **Solve[ Sin[x] == a, x ]**

Solve::ifun: Inverse functions are being used by Solve, so some solutions may not be found.

Out[9]=

It is important to realize that an equation such as actually has an infinite number of possible solutions, in this case differing by multiples of . However, Solve by default returns just one solution, but prints a message telling you that other solutions may exist. Section 3.4.5

discusses this in more detail.

There is no explicit "closed form" solution for a transcendental equation like this.
In[10]:= **Solve[ Cos[x] == x, x ]**

Solve::tdep: The equations appear to involve transcendental functions of the variables in an essentially non-algebraic way.

Out[10]=

You can find an approximate numerical solution using FindRoot, and giving a starting value for x.
In[11]:= **FindRoot[ Cos[x] == x, {x, 1} ]**

Out[11]=

Solve can also handle equations involving symbolic functions. In such cases, it again prints a warning, then gives results in terms of formal inverse functions.

Mathematica returns a result in terms of the formal inverse function of f.
In[12]:= **Solve[ f[x^2] == a, x ]**

InverseFunction::ifun: Warning: Inverse functions are being used. Values may be lost for multivalued inverses.

Out[12]=

Solving sets of simultaneous equations.

You can also use Mathematica to solve sets of simultaneous equations. You simply give the list of equations, and specify the list of variables to solve for.

Here is a list of two simultaneous equations, to be solved for the variables and

.
In[13]:= **Solve[{a x + y == 0, 2 x + (1-a) y == 1}, {x, y}]**

Out[13]=

Here are some more complicated simultaneous equations. The two solutions are given as two lists of replacements for x and y.
In[14]:= **Solve[{x^2 + y^2 == 1, x + 3 y == 0}, {x, y}]**

Out[14]=

This uses the solutions to evaluate the expression x+y.
In[15]:= **x + y /. %**

Out[15]=

Mathematica can solve any set of simultaneous linear equations. It can also solve a large class of simultaneous polynomial equations. Even when it does not manage to solve the equations explicitly, Mathematica will still usually reduce them to a much simpler form.

When you are working with sets of equations in several variables, it is often convenient to reorganize the equations by eliminating some variables between them.

This eliminates y between the two equations, giving a single equation for x.
In[16]:= **Eliminate[{a x + y == 0, 2 x + (1-a) y == 1}, y]**

Out[16]=

If you have several equations, there is no guarantee that there exists any consistent solution for a particular variable.

There is no consistent solution to these equations, so Mathematica returns {}, indicating that the set of solutions is empty.
In[17]:= **Solve[{x==1, x==2}, x]**

Out[17]=

There is also no consistent solution to these equations for almost all values of a.
In[18]:= **Solve[{x==1, x==a}, x]**

Out[18]=

The general question of whether a set of equations has any consistent solution is quite a subtle one. For example, for most values of a, the equations {x==1,x==a} are inconsistent, so there is no possible solution for x. However, if a is equal to 1, then the equations do have a solution. Solve is set up to give you generic solutions to equations. It discards any solutions that exist only when special constraints between parameters are satisfied.

If you use Reduce instead of Solve, Mathematica will however keep all the possible solutions to a set of equations, including those that require special conditions on parameters.

This shows that the equations have a solution only when a==1. The notation x==1&&a==1 represents the requirement that bothx==1anda==1 should be True.
In[19]:= **Reduce[{x==1, x==a}, x]**

Out[19]=

This gives the complete set of possible solutions to the equation. The answer is stated in terms of a combination of simpler equations. && indicates equations that must simultaneously be true; || indicates alternatives.
In[20]:= **Reduce[a x - b == 0, x]**

Out[20]=

This gives a more complicated combination of equations.
In[21]:= **Reduce[a x^2 - b == 0, x]**

Out[21]=

Functions for solving and manipulating equations.