   Next: About this document ... Up: lab_template Previous: lab_template

Subsections

# Solving Equations and Differentiation

## Solving Equations

The best method for solving equations is to use Maple's solving capabilities. First, a plot of the function or expression is useful then you can use the Maple solve command. The following illustrates how to find the roots of a function. > f := x-> 2*x^3-5*x^2-2*x+5;
> solve(f(x)=0,x);

If you forget to type in an equation and only type in an expression without setting it equal to something, Maple automatically sets the expression to zero. In the examples below, you can see some of the solving capabilities of Maple.
> plot({x^2+2*x-1,x^2+1},x=-2..2);
> solve(x^2+2*x-1=x^2+1,x);

Unfortuately, many equations cannot be solved analytically. For example, even the relatively simple equation sin(x) = x/2 has no analytical solution. In this case, the only possibility is to solve it numerically. In Maple, the command to use is fsolve. When the solve command is used the output looks like:
> solve(sin(x)=x/2,x);

This is not incorrect, as some of the zeros of a function may be imaginary and others may be real. However, it is much better to solve numerically as shown below:
> fsolve(sin(x) = x/2, x);

A plot of both equations on the same graph will show that this solution is not complete. There are two other intersection points that the fsolve command did not output. The fsolve command allows us to solve the equation in a range of values by replacing the second argument with what looks like a plot range. This method only works for the fsolve command and not for the solve command. To find the corresponding values of the intersection points, you only need to plug back into one of the functions.
> f:=x->sin(x);
> g:=x->x/2;
> plot({f(x),g(x)},x=-2*Pi..2*Pi);
> a:=fsolve(f(x)=g(x),x=-3..-1);
> b:=fsolve(f(x)=g(x),x=-1..1);
> c:=fsolve(f(x)=g(x),x=1..3);
> evalf(f(a));
> evalf(f(b));
> evalf(f(c));

Once you have solved an equation, you may want to use the output or the solution later. In order to label the output to a solution, you need to assign a label in the same line as the solve or fsolve command as seen above. It is different, however, if there is more than one output the the solve command. For example:
> f:=x->x^2+2*x-5;
> answer:=solve(f(x)=0,x);
> f(answer);
> f(answer);

Here an expression was defined first and then the solution was assigned to the label answer''. Note that there was more than one solution. In order to substitute the answer that was listed first back into the expression, the subs command was used and  was added onto the variable name to distinguish the first solution from the second.

## The Limit Definition of the Derivative

The limit definition of the derivative of often written as is defined as: The derivative of a function at a point , often written , can be interpreted algebraically as the following limit It can be interpreted geometrically as the slope of the tangent line to the graph of at and functionally as the instantaneous rate of change of at . Probably the second and third interpretations are the most important; they are certainly closer to what makes the derivative useful. In this lab, we will use Maple to explore each of these different aspects of the derivative. You can use the definition and the Maple limit command to compute derivatives directly, as shown below. You can also compute derivatives using Maple's diff or D command. The example below shows how to use the limit definition of derivative to find the derivative and then to find with Maple.
> f := x -> x^2+3*x+5;
> limit((f(x+h)-f(x))/h,h=0);
> limit((f(1+h)-f(1))/h,h=0);


## The Maple D and diff commands

These commands can be summarized as follows.
• The D operator acts on a function to produce the derivative of that function.
• The diff command acts on an expression and differentiates that expression with respect to a variable specified by the user.

To use the D operator to compute the derivative of a function:

> f:=x->x^2;
> D(f)(x);


If you want to evaluate the derivative at a specific value of you can use the following forms of the D operator.

> D(f)(2);

To use the plotting:
> plot(D(f)(x),x=-2..2);


The D operator cannot be used on expressions, for example trying to use it to differentiate the expresssion we defined above results in an error.

> p:=3*x+2;
> D(p);

The following is also an error.
> D(f(x));


To differentiate expressions, you need to use the diff command. Here is an example.

> diff(p,x);

The diff command can also be applied to functions. However, the result of the diff command is an expression, not a function. This means that computing the value of the derivative at a specific value of requires you to use the subs command. The following example shows how to calculate using the diff command.
> subs(x=2,diff(f(x),x));


## Finding Horizontal Tangents

Combining Maple's capabilities for solving equations and finding derivatives, it should be very simple to find points on the graph of a functions where the tangent line is horizontal. For instance, the solve command finds the values where the horizontal tangents occur and plugging the values back into the function gives the corresponding values as seen in the example below. A plot of the derivative may be necessary to insure that you have found all zeros of the derivative. Remember, if the solve command fails, try fsolve.
> f := x-> 2*x^3-3*x^2-36*x+12;
> plot(D(f)(x),x=-10..10);
> a:=evalf(solve(D(f)(x)=0,x));
> f(a);
> f(a);


## Exercises

1. For the functions and , plot both functions on the same graph using an -domain that clearly shows all the intersection points and then find the and coordinates of the intersection points.

2. Find the dervative of the function using the limit definition of the derivative, the diff command and then the D command and then use all three methods to find the slope of at .

3. Find the equation of the line tangent to the graph of the function at . Include a plot of the function and the tangent line on the same graph over the interval .   Next: About this document ... Up: lab_template Previous: lab_template
Jane E Bouchard
2011-01-13