0% found this document useful (0 votes)
71 views

Intro Ducci On Mathematica

This document provides an introductory tutorial to Mathematica 8. It is designed as interactive course material to be worked through by students at PCs over approximately 90 minutes. The tutorial covers getting started with Mathematica, basics of the Mathematica syntax including data types, the Wolfram|Alpha integration, symbolic calculus basics, input/output, definitions and assignments, functions and patterns, shortcuts, and syntax highlighting. Comments and the kernel functionality are also explained.

Uploaded by

Abi
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
71 views

Intro Ducci On Mathematica

This document provides an introductory tutorial to Mathematica 8. It is designed as interactive course material to be worked through by students at PCs over approximately 90 minutes. The tutorial covers getting started with Mathematica, basics of the Mathematica syntax including data types, the Wolfram|Alpha integration, symbolic calculus basics, input/output, definitions and assignments, functions and patterns, shortcuts, and syntax highlighting. Comments and the kernel functionality are also explained.

Uploaded by

Abi
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 37

Introductory Tutorial

to Mathematica 8
Author: Harald Hller
version: English version v1.0
last modified: 12.01.2012
Licence: Creative Commons Licence by-nc-sa 3.0 at
This introductory tutorial to Mathematica 8 is designed as interactive course material and should ideally be worked
through together with students at PCs. Each student should progress sequentially through all Input/Output. Duration: approx. 90min.

Getting Started with Mathematica


Mathematica Basics
Some Basic Remarks
Mathematica allows analytical, numerical and grafical calculations. Mathematica offers a relatively
simple and instructional logical syntax. Most of the known functions do have a syntactically similar
correspondent in Mathematica. The argument of a function always is to set inside square brackets
(e.g. Sin[x]). If you operate on some mathematical expression, also the operation is executed onto
brackets and the arguments are separated by commas (e.g. Integrate[f[x],x]).

Data types
The standard data type, the Notebook (.nb) provides an interactive interface in which calculations are
possible. Under File -> Save As there are some further data formats to which the notebooks can be
exported (.tex, .html, .txt, .pdf). Also single expressions can be copied as LATEX code and be pasted
into a .tex-file. Notebooks can be opened by Mathematica (interactive) or the MathReader (read
only) respectively the Wolfram CDF Player which is available as free internet browser plugin.
The forward- and backward-compatibility is given partially, which means that notebooks created in
Mathematica 7 will usually run in Mathematica 8, the opposite direction will rather cause problems.

introductory_mathematica_8_tutorial_expand.nb

Wolfram|Alpha
Wolfram|Alpha is a net-based frontend for Mathematica usage with a huge database of knowledge.
Describing as scientific Google does not really cover its whole power. Either you access it via the
website www.wolframalpha.com or there is even an integration to Mathematica. Beginning an Input
with the = sign is interpreted as Wolfram|Alpha query. E.g: = cosmological redshift z=3
cosmological redshift z = 3
Result
Assuming redshift-wavelength formula Use cosmological redshift instead
Calculate emitted wavelength
observed wavelength:

575 nm

Assuming observed wavelength and emitted wavelength

Use observed frequency and emitted frequency instead

Also include: dark energy density, matter density and radiation density

universe model Hubble parameter

Input information:

redshift-wavelength formula
redshift

observed wavelength

575 nm nanometers

Result:

More units

emitted wavelength

143.8 nm nanometers
5.659 106 inches
0.1438 m micrometers

Equation:

1z

introductory_mathematica_8_tutorial_expand.nb

1z

emitted wavelength

redshift

observed wavelength

Cosmological results:

More

time ago lookback time

11.5 billion years

time since big bang

2.19 billion years

21.1 billion ly light years


6480 Mpc megaparsecs

distance comoving

2 1023 km kilometers
1.24 1023 miles

fraction of total observable radius

0.453

scale factor

0.25 current value

epoch

matter dominated, post-recombination

radiation temperature

10.9 K kelvins

based on 5year WMAP data and LambdaCDM model; current universe age: 13.7
billion years

Timeline:
Big Bang

5 billion years

first galaxies form

emitted wavelength

10 billion years

present

Earth forms

143.8 nm nanometers
5.659 106 inches
0.1438 m micrometers

Input/Output I
The input of commands in Mathematica is done via the Enter key (not Return), which is either
located at the numeric keypad or it can be typed in as Shift + Return. Each Output and Input is numbered by Mathematica.

introductory_mathematica_8_tutorial_expand.nb

The input of commands in Mathematica is done via the Enter key (not Return), which is either
located at the numeric keypad or it can be typed in as Shift + Return. Each Output and Input is numbered by Mathematica.
3
3
In this way one can refer to previous In- and Output
In2 Out2
6
respectively one can refer to the last output via %.

Help / Documentation Center


Mathematica provides a compresehsive help and documentation center which gets completed by the
online-documentation. By entering the F1 key, the documentation center is opened; one can even
mark an expression in the notebook and will get to the help entry for this term via F1.
For a brief summary of syntax usage and functionality of a function, one can also just enter question
mark + function to get an overwiev.
? Integrate
Integrate f , x gives the indefinite integral f d x.

Integrate f , x, xmin , xmax gives the definite integral

xmax

f d x.

xmin

Integrate f , x, xmin , xmax , y, ymin , ymax , gives the multiple integral

xmax

xmin

With ?xyz* one gets a list of functions starting with xyz.

d x

ymax

ymin

dy f .

introductory_mathematica_8_tutorial_expand.nb

? Integ*
System`

Integer
IntegerDigits

IntegerExponent

IntegerPart

IntegerQ

IntegerString

IntegerLength

IntegerPartitions

Integers

Integral

Integrate

The Kernel
Mathematica "remembers" all definitions, In- and Outputs inside a Notebook. As described some
lines before, one can refer to previous terms easily. As long as the Kernel is running, all these definitions are available.
? Random
Random gives a uniformly distributed pseudorandom Real in the range 0 to 1.
Randomtype, range gives a pseudorandom number of the specified type, lying in the
specified range. Possible types are: Integer, Real and Complex. The default range is 0 to 1. You

can give the range min, max explicitly; a range specification of max is equivalent to 0, max.

a RandomInteger, 100
b RandomInteger, 100
c RandomInteger, 100
22
53
3
abc
78
If you want to delete entires in this memory of definitions, the Clear[Arguments] command can be
used. If you want to clear all kernel memory, i.e. kill and restart the kernel, then go to Evaluation ->
Quit Kernel -> Local. This also terminates running calculations which comes in handy, when you
realize that a calculation seems to have got stuck and you dont want to wait any longer.

introductory_mathematica_8_tutorial_expand.nb

? Clear
Clearsymbol1 , symbol2 , clears values and definitions for the symboli .
Clear" form1 ", " form2 ", clears values and
definitions for all symbols whose names match any of the string patterns formi .

Cleara, b, c
abc
abc
If one wishes not to quit ther kernel but delete all previous definitions, (of all currently opened Mathematica-Notebooks) the syntax is as follows.
Clear"Global`"

Comments
Comments (i.e. non to be executed, read-only text parts) inside Input-lines can by typed inside round
brackets and stars.
Sqrt2

Square root of 2 can also be typed via "Ctrl 2":

True

Shortcuts
There is a number of shortcuts and hotkeys in order to simplify typing formulas to the notebook.
Some examples:
ai Subscript via "Ctrl " ;
1
Fraction lines via "Ctrl " ;
2

"Esc int Esc" fx x "Esc dd Esc" ;


"Esc inf Esc" ;
Also greek letters via "Esc Letter Esc" ;

introductory_mathematica_8_tutorial_expand.nb

Syntax Highlighting
When typing formulae in a Mathematica notebook, syntax highlighting is very convenient. Known
Functions, trailing brackets, arguments etc. are recognoized and highlighted; either as positive feedback that everything is alright or - usually even more helpful - as negative feedback that there is some
syntactic problem. In case the syntax is faulty, klick on the plus-sign on the right hand side of the
input line to get more information about what might be wrong.
SolveSinx_ 3 Cosx_ 1, y

Symbolic Calculus in Mathematica - Basics


Input/Output II
The decimal separator is the dot (.) verwendet, not the comma. A space character is interpreted as
scalar multiplicator (times) and replaced by a cross.
3.14 3.14
9.8596

Definitions and Assignments


The simplest form of assigning a right hand side to a left hand side was already presented, namely the
equal sign.
?=
lhs rhs evaluates rhs and assigns the result
to be the value of lhs. From then on, lhs is replaced by rhs whenever it appears.

l1 , l2 , r1 , r2 , evaluates the ri , and assigns the results to be the values of the corresponding li .

d1
1
One can surpress the otput of an assignment by adding a colon. This assigment is evaluated not until
the definition is being used somewhere later in the Notebook. This assignment is especially useful
when the right hand side is somewhat elongate and evaluation costs computation time.

introductory_mathematica_8_tutorial_expand.nb

? :=
lhs : rhs assigns rhs to be the delayed value of lhs. rhs is maintained in
an unevaluated form. When lhs appears, it is replaced by rhs, evaluated afresh each time.

f : 1

Functions + Pattern
One main powerfulness of CAS such as Mathematica is symbolic calculus, i.e. working with functions. There is one little syntax peculiarity of Mathematica concerning this topic. Objects with arguments in the form f[x] are static; if you want the argument to be variable, then definition of the function must add an underscore (see Mathematica documentation center under keyword Pattern) to
each argument that should be a variable.
DONT :

hx : Sinx
h1
h1
DO :

gx_ : Sinx
g1
Sin1
As you can see, Mathematica does not necessarily output numerical values. In this case, the result is
an irrational number and one needs to specify further if and what numerical outoput is wanted.

Input/Output III
In order to force numerical output, one needs to type N[Argument,Digits].

introductory_mathematica_8_tutorial_expand.nb

Ng1, 100
0.8414709848078965066525023216302989996225630607983710656727517099
919104043912396689486397435430526959
There is also the possibility to append functions and objetcs that operate on the whole expression by
a double slash.
g1 N
0.841471

Basic Calculus (Scalars)


An addition of objetcs (scalars as well as vectors and matrices) is typed - as expected - via the plus (+)
key, subtraction via minus (-). Scalar multiplication can be set via the star (*) or a blank character,
division is typed viaone slash.
22
4
22
4
22
1
Nested expressions are typed as with pen and paper using round brackets.
1 4 3 2 1 1 1
5
2

10

introductory_mathematica_8_tutorial_expand.nb

Clear"Global`"
please enter this input to delete all previous definitions
which would partially conflict with the upcoming section

Basic Commands - I/O III


Until now we have seen Mathematica mainly as somewhat sophisticated calculator. However its true
power abounds when applying such a computer algebra system to complex problems such as integrals, differential equations or algebraic system of equations that are hardly or even not solvable at
all with pen and paper.
E.g.: Differentiation
Just to get used to the syntax we regard a simple example. The differentiation symbol is the capital D:
?D
D f , x gives the partial derivative f x.

D f , x, n gives the multiple derivative n f xn .


D f , x, y, differentiates f successively with respect to x, y, .
D f , x1 , x2 , for a scalar f gives the vector derivative f x1 , f x2 , .
D f , array gives a tensor derivative.

When we compute the first derivative of a general function F[x], we get the output F'[x]. As described
in the Mathematica documentation, the arguments (function and variable) of this operation are
separated by a comma.
DFx, x
F x
We have to further specify the function to be derived in order to get a calculated output of course. As
an example we define a function G in and compute ther first partial derivative with respecto to one
variable,
Gx_, y_ : x y
DGx, y, y
x

introductory_mathematica_8_tutorial_expand.nb

the second partial derivative


DGx, y, x, 2
0
or subsequent derivatives with respect to a list of variables.
DGx, y, x, y
1

Simplify und Expand


? Simplify
Simplifyexpr performs a sequence of algebraic and other transformations on expr, and returns the simplest form it finds.
Simplifyexpr, assum does simplification using assumptions.

? FullSimplify
FullSimplifyexpr tries a wide range of transformations on
expr involving elementary and special functions, and returns the simplest form it finds.
FullSimplifyexpr, assum does simplification using assumptions.

? Expand
Expandexpr expands out products and positive integer powers in expr.
Expandexpr, patt leaves unexpanded any parts of expr that are free of the pattern patt.

E.g.: Polynom
x x 1^3 x^2 3 x^5 x 1 x^2 x x^2 x 1^2^2
1 x x2 x 1 x3 1 x2 x x2 3 x x2
2

11

12

introductory_mathematica_8_tutorial_expand.nb

FullSimplify
2 x2 5 x 3 x 1 3 x5 x
Expand
2 5 x2 3 x3 x4 243 x5 405 x6 270 x7 90 x8 15 x9 x10

Plots / Graphics / Figures - Basics


Plot
A major scientific tool in Mathematica is the manifold of plotting routines that help you to depict
mathematical results grafically. The basic syntax for 2D graphs (as they are of course also done by
muich simpler programs like gnuplot) is Plot[Argument, {Variables, Boundaries}]. These boundaries
define the abscissa, the ordinate is scaled with the additional argument, namely PlotRange ->
{Boundaries}.
? Plot*
System`

Plot

Plot3Matrix

PlotJoined

PlotLayout

PlotPoints

PlotRangeClipping

Plot3D

PlotDivision

PlotLabel

PlotMarkers

PlotRange

PlotRangePadding

PlotRegion
PlotStyle

introductory_mathematica_8_tutorial_expand.nb

Plot Functions
PlotExpx, x, 0, 10

8000

6000

4000

2000

10

PlotSinx x, x, 10, 10
1.0

0.8

0.6

0.4

0.2

10

5
0.2

10

13

14

introductory_mathematica_8_tutorial_expand.nb

Plot1 x^2, x, 10, 10, PlotRange 0, 5,


PlotStyle Thick, Dashed, Background White, Filling Automatic
5

10

10

Dynamic Graphics
The two major commands for interactice plotting are:
? Dynamic
Dynamicexpr represents an object that displays as the dynamically updated current value
of expr. If the displayed form of Dynamicexpr is interactively changed or edited, an assignment
expr val is done to give expr the new value val that corresponds to the displayed form.
Dynamicexpr, None does not allow interactive changing or editing.
Dynamicexpr, f continually evaluates f val, expr during interactive changing or editing of val.
Dynamicexpr, f , fend also evaluates fend val, expr when interactive changing or editing is complete.

Dynamicexpr, fstart , f , fend also evaluates fstart val, expr when interactive changing or editing begins.

? Manipulate
Manipulateexpr, u, umin , umax generates a version
of expr with controls added to allow interactive manipulation of the value of u.
Manipulateexpr, u, umin , umax , du allows the value of u to vary between umin and umax in steps du.
Manipulateexpr, u, uinit , umin , umax , takes the initial value of u to be uinit .
Manipulateexpr, u, uinit , ulbl , labels the controls for u with ulbl .
Manipulateexpr, u, u1 , u2 , allows u to take on discrete values u1 , u2 , .
Manipulateexpr, u, , v, , provides controls to manipulate each of the u, v, .
Manipulateexpr, cu u, , cv v, ,
links the controls to the specified controllers on an external device.

introductory_mathematica_8_tutorial_expand.nb

ManipulatePlota Sina x x,
x, 7, 7, PlotRange 15, 15, a, 0, 10, 1

15

10

10

15

Graphics und Show


? Graphics
Graphics primitives, options represents a twodimensional graphical image.

? Show
Showgraphics, options shows graphics with the specified options added.
Showg1 , g2 , shows several graphics combined.

15

16

introductory_mathematica_8_tutorial_expand.nb

ShowPlotSinx, x, 0, Pi,
AspectRatio Automatic, PlotStyle Thick, Gray,
Background White, PlotRange 0, 2.5,
PlotExpx x x^2, x, 0, Pi, AspectRatio Automatic,
PlotStyle Thick, Blue,
Background White, PlotRange 0, 2.5,
GraphicsThick, Gray, Arrow0.2, Sin0.2, 0.2, 2,
GraphicsThick, Gray, Arrow.5, Sin.5, .5, 2,
GraphicsThick, Gray, Arrow.8, Sin.8, .8, 2,
GraphicsThick, Gray, Arrow1.1, Sin1.1, 1.1, 2,
GraphicsThick, Gray, Arrow1.1, Sin1.1, 1.1, 2,
GraphicsThick, Blue, Arrow1.7, Exp1.7 1.7 1.7^2,
2.5, Exp1.7 1.7 1.7^2,
GraphicsThick, Blue, Arrow1.8, Exp1.8 1.8 1.8^2,
2.5, Exp1.8 1.8 1.8^2,
GraphicsThick, Blue, Arrow1.9, Exp1.9 1.9 1.9^2,
2.5, Exp1.9 1.9 1.9^2,
GraphicsThick, Blue, Arrow2, Exp2 2 2^2,
2.5, Exp2 2 2^2,
GraphicsThick, Red, Circle2, 1,
GraphicsThick, Red, Circle2, 1, 1.1, Pi 8, 0,
GraphicsThick, Red, Circle2, 1, 1.1, Pi 8, 2 Pi 8,
GraphicsThick, Red, Circle2, 1, 1.1, 3 Pi 8, 4 Pi 8,
GraphicsThick, Red, Circle2, 1, 1.1, 5 Pi 8, 6 Pi 8,
GraphicsThick, Red, Circle2, 1, .9, Pi 8, 0,
GraphicsThick, Red, Circle2, 1, .9, Pi 8, 2 Pi 8,
GraphicsThick, Red, Circle2, 1, .9, 3 Pi 8, 4 Pi 8,
GraphicsThick, Red, Circle2, 1, .9, 5 Pi 8, 6 Pi 8

introductory_mathematica_8_tutorial_expand.nb

17

2.5

2.0

1.5

1.0

0.5

0.0

0.5

1.0

1.5

2.0

2.5

3.0

Linear Algebra - Calculus on Lists,


Vectors, Matrices etc.
On The Difference Between Lists and Vectors
Handling Vectors, Matrices and Tensors
Whenever working with CAS respectively even more basic level of data processing such as using
common programming languages, one has to keep in mind that a list of entries does not necessarily
make a vector, nor does a two dimensional list make a tensor etc. From algebra we know that all
tensorial quantities are only defined properly when also defining a base. This standard base is - as
expected - the cartesian n-dimensional unit base if nothing else is specified.
A list of functions, numbers or otherwise entries is input via curly brackets and the comma as separator.
vec a, b, c, d
a, b, c, d

18

introductory_mathematica_8_tutorial_expand.nb

? Dimensions
Dimensionsexpr gives a list of the dimensions of expr.
Dimensionsexpr, n gives a list of the dimensions of expr down to level n.

Dimensionsvec
4
Of course this concept can be used to arbitrarily interlace lists of lists etc.
tensor : a, b, c, a, b, c, a, b, c, a, b, c,
a, b, c, a, b, c, a, b, c, a, b, c
Dimensionstensor
2, 4, 3
E.g.: The Identity Matrix
? IdentityMatrix
IdentityMatrixn gives the n n identity matrix.

M : IdentityMatrix10
DimensionsM
10, 10

introductory_mathematica_8_tutorial_expand.nb

E.g.: Automatic Filling of Lists


? Table
Tableexpr,
Tableexpr,
Tableexpr,
Tableexpr,
Tableexpr,

imax generates a list of imax copies of expr.


i, imax generates a list of the values of expr when i runs from 1 to imax .
i, imin , imax starts with i imin .
i, imin , imax , di uses steps di.
i, i1 , i2 , uses the successive values i1 , i2 , .

Tableexpr, i, imin , imax , j, jmin , jmax , gives a nested list. The list associated with i is outermost.

H Table1 i j 1, i, 1, 15, j, 1, 15;


MatrixForm
1
1
2
1
3
1
4
1
5
1
6
1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15

1
2
1
3
1
4
1
5
1
6
1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16

1
3
1
4
1
5
1
6
1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17

1
4
1
5
1
6
1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18

1
5
1
6
1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19

1
6
1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20

1
7
1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21

? Subscript
Subscriptx, y is an object that formats as x y .
Subscriptx, y1 , y2 , formats as x y1 ,y2 , .

1
8
1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22

1
9
1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23

1
10
1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23
1
24

1
11
1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23
1
24
1
25

1
12
1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23
1
24
1
25
1
26

1
13
1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23
1
24
1
25
1
26
1
27

1
14
1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23
1
24
1
25
1
26
1
27
1
28

1
15
1
16
1
17
1
18
1
19
1
20
1
21
1
22
1
23
1
24
1
25
1
26
1
27
1
28
1
29

19

20

introductory_mathematica_8_tutorial_expand.nb

TableSubscriptm, i, j, i, 3, j, 3 MatrixForm
m1,1 m1,2 m1,3
m2,1 m2,2 m2,3
m3,1 m3,2 m3,3

Coordinate Systems + Vector Calculus


Since the topic of coordinate systems is especially interesting in physical and astrophysical applications, we want to take a look at the implemented functionalities in Mathematica. A number of packages in Mathematica are not loaded into the Kernel by default (which makes the startup faster) but
have to be loaded on demand. The VectorAnalysis package is one of those - other exaples can be
found here: tutorial/PackagesForSymbolicMathematics.

The VectorAnalysis Package


Loading the package
Needs"VectorAnalysis`"
enables a number of additional commands. E.g. with
CoordinateSystem
Cartesian
one can inquire the momentarily used coordinate system. As mentioned before, this is cartesian by
default. When we change to a spherical coordinate system (remark: American coordinate ranges by
default)
SetCoordinatesSphericalr, ,
Sphericalr, ,
and compute the divergence of the unit vector in r-direction, (1,0,0), we get the well known result.
Div1, 0, 0
2
r

introductory_mathematica_8_tutorial_expand.nb

21

JacobianMatrixSpherical MatrixForm
Cos Sin r Cos Cos r Sin Sin
Sin Sin r Cos Sin r Cos Sin
Cos
r Sin
0
Just as further example for some functionalities of that package, we calculate the unit volume of a
cylindre with radius and height one by the differential geometric rule as the integral over the Jacobian determinant.
SetCoordinatesCylindricalr, , z
Cylindricalr, , z
IntegrateJacobianDeterminantCylindrical,
r, 0, 1, , 0, 2 , z, 0, 1

Just to check, we compute the divergence of the unit vector in x-direction for cartesian coordinates
SetCoordinatesCartesian
CartesianXx, Yy, Zz
Div1, 0, 0
0
which vanishes of course.

Scalar Product, Cross Product


Even without the VectorAnalysis package activated, there are vector commands loaded by default,
such as the scalar (inner) product of two vectors, input via the dot (.).
vec a, b, c
a, b, c

22

introductory_mathematica_8_tutorial_expand.nb

vec.vec
a2 b2 c2
With the star * the entries of the list are multiplied pairwise.
vec vec
a2 , b2 , c2
The cross profuct (vector product) is set with the command Cross.
Crossvec, vec
0, 0, 0
Since a vector is parallel to itself, the cross product has to vanish.
Cross1, 0, 0, 0, 1, 0
0, 0, 1
The cross product between two unit base vectors in three dimensions yields the third unit base vector.

Solving Sets of Equations


optional: Solving Linear Sets of Equations
One of the most importand application of computer algebra and numerical codes is solving systems
of equations. Wether you want to solve a complicated switching circuit by solving the Kirchhoff rules
or you look for the numerical solution of differential equations, mathematically these problems
reduce to solving (huge) sets of linear equations or - even more basic - it means to invert big matrices.
Although CAS and numerical schemes are very powerful tools for computationally expensive problems, there are of course limitations as well. This shall be demonstrated with the help of a little example.
The main command for solving linear systems of equations is

introductory_mathematica_8_tutorial_expand.nb

23

? LinearSolve
LinearSolvem, b finds an x which solves the matrix equation m.x b.
LinearSolvem generates a LinearSolveFunction which can be applied repeatedly to different b.

which solves equations of the form m.x=b.


E.g.: Hilbert - Matrix
We ragard an arbitrary real vector b (filled with random real numbers between 0 and 100)
b TableRandomInteger, 100, n, 1, 15
31, 84, 7, 64, 98, 93, 11, 84, 30, 94, 51, 62, 75, 11, 9
and solve the system H.x=b, i.e. we want to determine the solution vector x. The matrix H was
defined above in subsection Automatic Filling of Lists .
LinearSolveH, b MatrixForm
3 734 071 933 185
751 754 904 677 280
37 714 435 829 386 680
828 864 504 280 612 800
9 984 977 804 685 904 020
74 111 982 969 427 603 200
362 869 877 239 232 683 320
1 221 043 396 114 981 604 880
2 889 102 402 497 160 559 710
4 847 877 238 595 981 990 400
5 737 637 243 526 263 875 800
4 682 450 439 923 249 994 000
2 507 202 365 724 910 608 300
792 653 569 836 038 568 000
112 131 660 922 158 101 400
On some computers this may take quite some time since the determinant of this matrix, also known
as Hilbert-Matrix, gets numerically very small with increasing problem size. From the theory of
linear algebra we know that a system of linear equations is only solvable when the determinant is
non-zero. Computationally, matrix inversion gets more and more expensive, the smaller the determinant.

24

introductory_mathematica_8_tutorial_expand.nb

? Det
Detm gives the determinant of the square matrix m.

DetH N
1.05854 10124
The determinant of the 15-dim. Hilbert-Matrix is practically zero. The higher dimensional H, the
smaller gets its determinant and the more difficult gets its inversion.
? HilbertMatrix
HilbertMatrixn gives the nn Hilbert matrix with elements of the form 1 i j 1.
HilbertMatrixm, n gives the mn Hilbert matrix.

DetHilbertMatrix50 N
1.392615568935140 101466
Clear"Global`"

Solving Coupled Systems of Equations


Clearly we are often also confronted with nonlinear sets of coupled equations, where we cannot
express the problem as sketched before. In this case, we apply the command
? Solve
Solveexpr, vars attempts to solve the system expr of equations or inequalities for the variables vars.
Solveexpr, vars, dom solves over the domain dom. Common choices of dom are Reals, Integers, and Complexes.

in which documentation entry it states already indicatively that it attempts to solve the problem
posed. Mathematica has implemented a variety of algorithms that are tried out.

introductory_mathematica_8_tutorial_expand.nb

E.g. : System of three linear equations


Solve3 a b, 2 a c, c 1, a, b, c
a

1
2

, b

3
2

, c 1

25

26

introductory_mathematica_8_tutorial_expand.nb

E.g.: System of three nonlinear equations


Solve3 a^3 b^2, 2 a^2 c^2, c^3 1, a, b, c
a

, b

2
a

, b

2
a

1
4
1
4

1
4
1
4

1
4
1
4

1
4
1
4

234
3

234

3
4

3
4

3
4

3
4

3
4

3
4

, b

, b

, b

, b

1
2

1
2

2
3

1
2

1
2

1
1

2
3

3
234

1
1

, b

, b

, b

, b

34

34

34

34

1
1

, c

, c

34

34

34

34

, c

, c

, c 1,

, c 1,

234

, c 1, a

, b

, c 1, a

, b

1
2

1
2

1
2

1
2

, c

, c

, c

, c

1
2

1
2

1
2

3 ,
3 ,

3 ,

3 ,

3 ,
3 ,

3 ,
3

Mathematica will list all found solutions also if they sometimes not usable for our purposes because
they are unphysical or violate some other conditions. In this case we frequently have to come back to
pen and paper in order to pick the right results.

introductory_mathematica_8_tutorial_expand.nb

27

Mathematica will list all found solutions also if they sometimes not usable for our purposes because
they are unphysical or violate some other conditions. In this case we frequently have to come back to
pen and paper in order to pick the right results.
Some more useful commands concerning linear algebra:
? RandomReal
RandomReal gives a pseudorandom real number in the range 0 to 1.
RandomRealxmin , xmax gives a pseudorandom real number in the range xmin to xmax .
RandomRealxmax gives a pseudorandom real number in the range 0 to xmax .
RandomRealrange, n gives a list of n pseudorandom reals.

RandomRealrange, n1 , n2 , gives an n1 n2 array of pseudorandom reals.

A RandomReal1, 1, 5, 5
0.474727, 0.388645, 0.461371, 0.938493, 0.950533,
0.830478, 0.139342, 0.801181, 0.0834233, 0.259048,
0.606089, 0.750426, 0.0652498, 0.50112, 0.928461,
0.903335, 0.207855, 0.368359, 0.051829, 0.0604615,
0.0480989, 0.359201, 0.128467, 0.236191, 0.0770996
? MatrixRank
MatrixRankm gives the rank of the matrix m.

MatrixRankA
5
? Eigenvalues
Eigenvaluesm gives a list of the eigenvalues of the square matrix m.
Eigenvaluesm, a gives the generalized eigenvalues of m with respect to a.
Eigenvaluesm, k gives the first k eigenvalues of m.
Eigenvaluesm, a, k gives the first k generalized eigenvalues.

EigenvaluesA
0.195529 1.07151 , 0.195529 1.07151 ,
0.325532 0.509743 , 0.325532 0.509743 , 0.269558

28

introductory_mathematica_8_tutorial_expand.nb

? Eigenvectors
Eigenvectorsm gives a list of the eigenvectors of the square matrix m.
Eigenvectorsm, a gives the generalized eigenvectors of m with respect to a.
Eigenvectorsm, k gives the first k eigenvectors of m.
Eigenvectorsm, a, k gives the first k generalized eigenvectors.

EigenvectorsA
0.144597 0.32199 , 0.647119 0. , 0.0223833 0.456948 ,
0.391941 0.0945935 , 0.127755 0.261649 ,
0.144597 0.32199 , 0.647119 0. , 0.0223833 0.456948 ,
0.391941 0.0945935 , 0.127755 0.261649 ,
0.0352433 0.480426 , 0.327103 0.0011468 ,
0.259849 0.380617 , 0.542821 0. , 0.276221 0.278583 ,
0.0352433 0.480426 , 0.327103 0.0011468 ,
0.259849 0.380617 , 0.542821 0. , 0.276221 0.278583 ,
0.320862, 0.369529, 0.656308, 0.272299, 0.505578
? CharacteristicPolynomial
CharacteristicPolynomialm, x gives the characteristic polynomial for the matrix m.
CharacteristicPolynomialm, a, x gives the generalized characteristic polynomial with respect to a.

CPx_ CharacteristicPolynomialA, x
0.116983 0.603626 x 0.979111 x2 1.36765 x3 0.529564 x4 x5
The eigenvalues of a matrix (resp. a linear map) are the zeros of the characteristical polynomial
which we can probe:
CPEigenvaluesA
2.22045 1016 0. , 2.22045 1016 0. ,
2.77556 1017 5.55112 1017 ,
2.77556 1017 5.55112 1017 , 1.73472 1017
... is approximately zweo.

Analysis - Working with Functions,


Derivatives, Integrals

introductory_mathematica_8_tutorial_expand.nb

29

Analysis - Working with Functions,


Derivatives, Integrals
Introductory Remarks
Mathematical Entities in Mathematica
As mentioned before, Mathematica has some syntactic idiosyncrasies. As an example we regard one
of the most beautiful formulas in mathematics:
DONT :

e^i pi 1
ei pi 1
DO :

E^I Pi 1
True
?E
E is the exponential constant e base of natural logarithms, with numerical value 2.71828.

?I
I represents the imaginary unit

1 .

? Pi
Pi is , with numerical value 3.14159.

30

introductory_mathematica_8_tutorial_expand.nb

? Infinity
Infinity or is a symbol that represents a positive infinite quantity.

optional: Natural Constants etc.


PhysicalConstants`
SpeedOfLight
299 792 458 Meter
Second
FineStructureConstant
0.00729735
AgeOfUniverse
4.7 1017 Second
MagneticFluxQuantum
2.06783 1015 Weber

Important Embedded Functions


Clearly the most important functions are implemented in Mathematica; as mentioned before all
these commands begin with a capital letter. Some Examples:
? Sin
Sinz gives the sine of z.

introductory_mathematica_8_tutorial_expand.nb

31

? Log
Logz gives the natural logarithm of z logarithm to base e.
Logb, z gives the logarithm to base b.

? Gamma
Gammaz is the Euler gamma function z.
Gammaa, z is the incomplete gamma function a, z.
Gammaa, z0 , z1 is the generalized incomplete

gamma function a, z0 a, z1 . It is also a unit of magnetic flux density.

? DiracDelta
DiracDeltax represents the Dirac delta function x.
DiracDeltax1 , x2 , represents the multidimensional Dirac delta function x1 , x2 , .

Series
In order to get used to handling series, we examine a frequent application, namely determining the
limit of a series.
Clear"Global`"
please enter this input to delete all previous definitions
which would partially conflict with the upcoming section
We can define series analogously to functions
an_ 1 1 n^n

or recursively
bn_ : bn 1 bn 2
b0 : 1
b1 : 2
b2 : 2

32

introductory_mathematica_8_tutorial_expand.nb

and determine the limit or generate a list of entries.


? Limit
Limitexpr, x x0 finds the limiting value of expr when x approaches x0 .

Limitan, n Infinity

Tablebn, n, 0, 10
1, 2, 2, 4, 8, 32, 256, 8192,
2 097 152, 17 179 869 184, 36 028 797 018 963 968
Of course also sums are implemented im Mathematica
? Sum
Sum f , i, imax evaluates the sum f .
imax

Sum f , i, imin , imax starts with i imin.


Sum f , i, imin , imax , di uses steps di.
Sum f , i, i1 , i2 , uses successive values i1 , i2 , .
i1

Sum f , i, imin , imax , j, jmin , jmax , evaluates the multiple sum f .


imax

Sum f , i gives the indefinite sum f .

jmax

iimin j jmin

Sumx^n n^2, n, 0, Infinity


BesselI0, 2

? BesselI
BesselIn, z gives the modified Bessel function of the first kind In z.

introductory_mathematica_8_tutorial_expand.nb

Sumx^k k, k, 0, n
x Gamma1 n, x
n
and especially useful is also the Series command which expands a function into its Taylor series.
? Series
Series f , x, x0 , n generates a power series expansion for f about the point x x0 to order x x0 n .

Series f , x, x0 , n x , y, y0 , n y , successively finds series expansions with respect to x, then y, etc.

SeriesSinx, x, 0, 15
x

x3

x5

x7

x9

5040 362 880


x13
x15

Ox16
39 916 800 6 227 020 800 1 307 674 368 000
6
120
11
x

Differentiation and Integration


Integration
Basic differentiation syntax was presented already in section I/O III. Integration commands are
analogous
Clear"Global`"
please enter this input to delete all previous definitions
which would partially conflict with the upcoming section
? Integrate
Integrate f , x gives the indefinite integral f d x.

Integrate f , x, xmin , xmax gives the definite integral

xmax

f d x.

xmin

Integrate f , x, xmin , xmax , y, ymin , ymax , gives the multiple integral

xmax

xmin

d x

ymax

ymin

dy f .

33

34

introductory_mathematica_8_tutorial_expand.nb

E.g.: Some Integrals


IntegrateSinx, x
Cosx
IntegrateE^x^2, x
1
2

Erfix

? Erfi
Erfiz gives the imaginary error function erfiz i.

IntegrateSinx Cosx, x, Pi, Pi


0
IntegrateE^x, x, 0, Infinity
1

Differential Equations
Solving Ordinary Differential Equations (ODEs)
Solving differential equations is another major strenght of CAS. While solving an ODE can require a
sophisticated ansatz and several analytic techniques, Mathematica is a lot faster in giving results
with help of the command DSolve.
Clear"Global`"
please enter this input to delete all previous definitions
which would partially conflict with the upcoming section

introductory_mathematica_8_tutorial_expand.nb

35

E.g.: Solving ODEs without boundary condition


? DSolve
DSolveeqn, y, x solves a differential equation for the function y, with independent variable x.
DSolveeqn1 , eqn2 , , y1 , y2 , , x solves a list of differential equations.
DSolveeqn, y, x1 , x2 , solves a partial differential equation.

DSolvey'x yx 0, yx, x
yx x C1
E.g.: Solving ODEs with boundary condition
DSolvez'x zx 0, z0 1, zx, x
zx x
E.g.: Solving systems of ODEs
eqn1 : s'x tx 0
eqn2 : t'x sx 0
DSolveeqn1, eqn2, sx, tx, x
sx

1
2

tx

x 1 2 x C1
1
2

1
2

x 1 2 x C1

x 1 2 x C2,
1
2

x 1 2 x C2

E.g.: Numerical solutions of ODEs


Whenever Mathematica fails to find a closed analytic solution of an ODE, such as the following
example shows
DSolvey''x y'x yx 1 yx 0, yx, x
DSolveyx 1 yx y x y x 0, yx, x
there is the possibility to look for solutions numerically. Iterative algorithms implemented in Mathematica are capable of finding numerical solutions of ODEs.

36

introductory_mathematica_8_tutorial_expand.nb

there is the possibility to look for solutions numerically. Iterative algorithms implemented in Mathematica are capable of finding numerical solutions of ODEs.
? NDSolve
NDSolveeqns, y, x, xmin , xmax finds a numerical solution to the ordinary differential
equations eqns for the function y with the independent variable x in the range xmin to xmax .
NDSolveeqns, y, x, xmin , xmax , t, tmin , tmax finds a numerical solution to the partial differential equations eqns.
NDSolveeqns, y1 , y2 , , x, xmin , xmax finds numerical solutions for the functions yi .

sol NDSolveyx 1 yx y x y x 0, y0 1, y'0 2,


yx, x, 0, 15
yx InterpolatingFunction0., 15., x
The result is output as InterpolatingFunction meaning that the function is only known at those
points where the algorithm has provided a numerical result.
PlotEvaluateyx . sol, x, 0, 15,
PlotRange 1.5, PlotStyle Thick, Background White
1.5

1.0

0.5

10

12

14

0.5

1.0

1.5

CAS: Alternatives
Mathematica is a comprehensive computer algebra system, however it is a proprietary product.
There are several free open source alternatives:
Axiom:
Scilab:
Octave:
Maxima:
Wiris:
Geogebra:

http://axiom-wiki.newsynthesis.org/FrontPage
http://www.scilab.org/
http://www.gnu.org/software/octave/
http://maxima.sourceforge.net/
http://www.wiris.com/
http://www.geogebra.org/cms/

Mathematica is a comprehensive computer algebra system, however it is a proprietary product.


There are several free open source alternatives:
introductory_mathematica_8_tutorial_expand.nb

Axiom:
Scilab:
Octave:
Maxima:
Wiris:
Geogebra:

http://axiom-wiki.newsynthesis.org/FrontPage
http://www.scilab.org/
http://www.gnu.org/software/octave/
http://maxima.sourceforge.net/
http://www.wiris.com/
http://www.geogebra.org/cms/

37

You might also like