Tutorial: using the notebook and navigating the help system (long version) -- Sage Reference Manual v4.6.1 system:sage

Tutorial: using the notebook and navigating the help system (long version)

This worksheet is based on William Stein’s JPL09__intro_to_sage.sws worksheet and the Sage days 20.5_demo worksheet and aims to be an interactive introduction to Sage through exercises. You will learn how to use the notebook and call the help.

Making this help page into a worksheet

Go into the File menu at the top left of this window and click on Copy worksheet. Then, you can clear all output in the menu Action by clicking on Delete All Output.

Entering, Editing and Evaluating Input

To evaluate code in the Sage Notebook, type the code into an input cell and press shift-enter or click the evaluate link. Try it now with a simple expression (e.g., 2 + 3). The first time you evaluate a cell takes longer than subsequent times since a new Sage process is started:

{{{id=0| 2 + 3 /// 5 }}} {{{id=1| /// }}} {{{id=2| /// }}}

Create new input cells by clicking blue line that appears between cells when you move your mouse around. Try it now:

sage:

You can go back and edit any cell by clicking in it (or using the keyboard to move up or down). Go back and change your 2+2 above to 3 + 3 and re-evaluate it.

You can also edit this text right here by double clicking on it, which will bring up the TinyMCE Javascript text editor. You can even put embedded mathematics like this $sin(x) - y^3$ just like with LaTeX.

Help systems

There are various ways of getting help in Sage.

  • navigate through the documentation (there is a link Help at the top right of the worksheet),
  • tab completion,
  • contextual help.

We details in what follows the two last methods through examples.

Completion and contextual domcumentation

Start typing something and press the tab key. The interface tries to complete it with a command name. If there is more than one completion, then they are all presented to you. Remember that Sage is case sensitive, e.g. it differentiates lower case from lower case. Hence the tab completion of klein won’t show you the KleinFourGroup command that build the group \ZZ/2 \times \ZZ/2 as a permutation group. Try it on the next cells

{{{id=3| klein /// }}} {{{id=4| Klein /// }}}

To see documentation and examples for a command, type a question mark ? at the end of the command name and press the tab key as

KleinFourGroup?<tab>
sage:

Exercise A: What is the largest prime factor of 600851475143?

factor?<tab>
sage:

In the above exercise we do not affect any mathematical data to variables. The affectation in Sage is realizer through the = symbol as in:

{{{id=5| a = 3 b = 2 print a+b /// 5 }}}

Once an object is created, some methods apply to it. This means functions but instead of writing f(my_object) you write my_object.f().:

{{{id=6| p = 17 p.is_prime() /// True }}}

To know all methods of an object you can still use tab-completion. Write the name of the object followed by a dot and then press tab.:

a.<tab>
sage:

sage:

Exercise B: Create the Permutation 51324 and assign it to the variable p.

Permutation?<tab>
sage:

What is the inverse of p ?

p.inv<tab>
sage:

Does p have the pattern 123 ? What about 1234 ? And 312 ? (even if you don’t know what a pattern is, you should be able to find a command that does this).

p.pat<tab>
sage:

Some linear algebra

Exercise C: Use the matrix command to create the following matrix.

M = \left(\begin{array}{rrrr}
10 & 4 & 1 & 1 \\
4 & 6 & 5 & 1 \\
1 & 5 & 6 & 4 \\
1 & 1 & 4 & 10
\end{array}\right)

matrix?<tab>
sage:

Then using methods of the matrix

  1. Find the determinant of the matrix.
  2. Find the echelon form of the matrix.
  3. Find the eigenvalues of the matrix.
  4. Find the kernel of the matrix.
  5. Find the LLL decomposition of the matrix.
sage:

sage:

Now you know how to access the different methods of matrices:

  1. Create the vector v = (1,-1,-1,1).
  2. Compute the products: M*v and v*M.
sage: vector?<tab>
sage:

Note

Vectors in Sage are row vectors. A command such as eigenspaces might not return what you expect, so it is best to specify eigenspaces_left or eigenspaces_right instead. Same thing for kernel (left_kernel or right_kernel), and so on.

Some Plotting

The plot command allows you to draw plots of functions. Recall that you can access the documentation by pressing the tab key after writing plot? in a cell.:

sage: plot?<tab>
sage:

Here is a simple example:

{{{id=7| var('x') # make sure x is a symbolic variable plot(sin(x^2), (x,0,10)) /// }}}

Here is a more complicated plot. Try to change every single input to the plot command in some way, evaluating to see what happens:

{{{id=8| P = plot(sin(x^2), (x,-2,2), rgbcolor=(0.8,0,0.2), thickness=3, linestyle='--', fill='axis') show(P, gridlines=True) /// }}}

Above we used the show command to show a plot after it was created. You can also use P.show instead:

{{{id=9| P.show(gridlines=True) /// }}}

Try putting the cursor right after P.show( and pressing tab to get a list of the options for how you can change the values of the given inputs.:

{{{id=10| P.show( /// }}}

Plotting multiple functions at once is as easy as adding them together:

{{{id=11| P1 = plot(sin(x), (x,0,2*pi)) P2 = plot(cos(x), (x,0,2*pi), rgbcolor='red') P1 + P2 /// }}}

Symbolic Expressions

Here is an example of a symbolic function:

{{{id=12| f(x) = x^4 - 8*x^2 - 3*x + 2 f(x) /// x^4 - 8*x^2 - 3*x + 2 }}} {{{id=13| f(-3) /// 20 }}}

Exercise D: Define the symbolic function f(x) = x \sin(x^2). Plot f on the domain [-3,3] and colour it red. Use the find_root method to numerically approximate the root of f on the interval [1,2].:

sage:

Compute the tangent line to f at x=1.:

sage:

Plot f and the tangent line to f at x=1 in one image:

sage:

Exercise E (Advanced): Solve the following equation for y

y = 1 + x y^2

There are two solutions, take the one for which \lim_{x\to0}y(x) = 1. (Don’t forget to create the variables x and y!).:

sage:

Expand $y$ as a truncated Taylor series around $0$ and containing $n=10$ terms.:

sage:

Do you recognize the coefficients of the Taylor series expansion? You might want to use, or better yet, Sage’s command sloane_find which uses the online encyclopedia of integers:

sloane_find?<tab>
sage:

sage: