Mathematica Cheat Sheet

Algebra...

% math Use ssh to get to the % prompt In[1]:= Expand[(x+y)^2] Mathematica can expand an algebraic
Out[1]= x2 + 2 x y + y2 expression... or factor it back to a compact form.
In[1]:= Exit or hit control-d Leave Mathematica (when you’re ready to!)
In[2]:= Factor[%]
% math < sample.m > sample.lst Run Mathematica commands from sample.m Out[2]= (x + y)2
% more sample.lst (non-interactively) with output to sample.lst
In[3]:= Solve[x^2==81,x] Find the roots of an equation; note use of ==
Using Mathematica like a calculator... Out[3]={{x -> -9}, {x -> 9}} (rather than just =) in writing the equation.

In[2]:= 27.50-11.92 Mathematica as a good old calculator... hit In[4]:= Solve[x^2==-4,x] Imaginary numbers? No problem...
Out[2]= 15.58 ENTER (or shift-ENTER) after each command Out[4]= {{x -> -2I},{x -> 2I}}

In[3]:= 15! Large values are no problem; you could even In[5]:=Solve[{x+y==1,3x+y==2}] Mathematica can also solve systems of
Out[3]= 1307674368000 compute 1500 factorial if you wanted to 1 1 algebraic equations in multiple variables.
Out[5]= {{x -> -, y -> -}}
In[4]:= ?Log Need help with a function? Enter a ? followed 2 2
Log[z] gives the natural by the name of a Mathematica function. Not sure
logarithm of z (logarithm to of a function’s name? You can use a * to see Calculus...
base e). Log[b, z] gives the possible matches, e.g., ?L*
logarithm to base b. In[1]:= Limit[x/(Sqrt[x+1]-1),x->0] Evaluate a limit
Out[1]= 2
In[5]:= Log[10,3453.8] Note that Mathematica functions are case
Out[5]= 3.538 sensitive and begin with a capital letter. In[2]:= Dt[x^3+2x,x] Compute a total derivative
In[6]:= (4000/23)^3 Operations done on whole numbers are Out[2]= 2 + 3 x2
64000000000 always represented exactly when possible.
In[3]:= D[(x^2)(y^3)+4y+x+2,x] Partial derivatives work the same way
Out[6]= -----------
12167 % means “recall the last result” and //N Out[3]= 1 + 2 x y3
In[7]:= %//N means “provide an approximate numerical
In[4]:= D[x^3+2x,x,x] Take the 2nd derivative with respect to x
Out[7]= 5.26013 106 result”
Out[4]= 6 x
In[8]:= Sin[60 Degree] Function args must be put in square brackets. In[5]:= Integrate[3x^2+2x,x] Mathematica can also do integrals, just as
Sqrt[3] Trig functions are in radians by default. you’d expect.
Out[5]= x2 + x3
Out[8]= ------- Want a numeric value? Remember //N
2 Inverse functions? ArcSin[ ]/Degree In[6]:= Integrate[E^x,{x,0,1}] Definite integral are also easy to evaluate.
Out[6]= -1 + E
In[9]:= Sum[i/(i^i),{i,1,\ Numerically evaluate an infinite sum.
Infinity}]//N You can continue long Mathematica In[7]:= <<Calculus`VectorAnalysis` Cartesian space is the default, but not our only
Out[9]= 1.62847 commands lines with a \ at the end of a line In[8]:= SetCoordinates[\ option. For example, let’s find the surface area
Cylindrical]
In[10]:= BaseForm[223,2] Convert the value 223 (decimal) to base 2 of the parabola z=1+x2+y2 where x2+y2 <=1.
Out[8]= Cylindrical[Rr,Ttheta,Zz]
Out[10]//BaseForm= 110111112 (binary). In[9]:= Integrate[Sqrt[1+4Rr^2]\ Because of the nature of that restriction, it is
Rr,{Rr,0,1},{Ttheta,0,2Pi}]//N easier to work in cylindrical coordinates. We do
In[11]:= 16^^FAE7 + 16^^2C3E
Add FAE7 (hex) to 2C2E (hex); output by Out[9]= 5.33041 so via the vector analysis package (note the
Out[11]= 75557
default is in decimal, but you can then force backtick marks, not apostrophes, used when
In[12]:= BaseForm[%,16]
that output into hex, too, if you like. loading a package!). Package info is at
Out12//BaseForm= 1272516
http://documents.wolfram.com/v4/index20.html
Linear Algebra... Page 2

In[1]:= w={{a,b},{c,d}} Create a 2x2 matrix (we’re using symbols, but In[3]:=!!mydata.dat Work with (x,y) data points from an external file.
Out[1]= {{a, b}, {c, d}} you could equally easily use numeric values) 4.1 10.7 !!mydata.dat shows us the contents of the file.
[etc] Read in pairs of numbers from that file, storing
In[2]:= w.{x,y}=={k1,k2} Use a dot product to apply that matrix of In[4]:= newvals=ReadList[ \ the list of values by the name newvals. Plot the
Out[2]= {a x + b y, c x + d y} == coefficients to two variables to form a system "mydata.dat",{Number,Number}] dataset. Fit a line to the points & plot that.
{k1, k2} of two equations with constants {k1, k2} Out[4]= {{4.1,10.7},[etc]} Finally, overlay both and save as a gif
In[5]:= plot1=ListPlot[newvals]
In[3]:= Transpose[w]//MatrixForm Mathematica can easily do most standard
Out[5]= -Graphics- [not shown]
Out[3]//MatrixForm= a c linear algebra operations, for example, we
In[6]:= Fit[newvals,{1,x},{x}]
b d can easily transpose matrix w...
Out[6]= 5.14286 + 9.96429 x
In[4]:=Inverse[{{1,-1},{2,2}}] Or compute the inverse of a 2x2 numeric In[7]:= plot2=Plot[%,{x,1,8}]
1 1 1 1 matrix... Out[7]= -Graphics- [not shown]
Out[4]={{-, -},{-(-), -}} In[8]:= Show[plot1,plot2]
2 4 2 4 Out[8]= -Graphics-
In[9]:= Display["b.gif",%,"GIF"]
In[5]:= Det[{{a,b,c},{d,e,f},\ Or compute the determinant of a 3x3 Out[9]= -Graphics-
{g,h,i}}] symbolic matrix...
Out[5]= -(c e g) + b f g + c d h
Mathematica As A Programming Language...
- a f h - b d i + a e i
(* Approach No. 1 *) If Mathematica doesn’t have precisely what you
In[6]:= Table[If[EvenQ[i]||EvenQ[j]\ In addition to entering matrices on an element
w=Join[Table[0,{7}],Table[5,{7}],\ need (or what it has is overkill), you can always
,1,0],{i,3},{j,3}]//MatrixForm by element basis, Mathematica will also let us Table[10,{4}],{25}]; use Mathematica as a programming language
Out[6]= 0 1 0 construct matrices using rules, such as this
1 1 1 <<DiscreteMath`Combinatorica` and write your own code. For example, assume
example that sets elements of a 3x3 matrix to
0 1 0 x=Union[KSubsets[w,7]]; you have a pile of 5, 10 and 25 pound weights.
be 1 if the column or row is an even number. Select[x,(Plus@@##)<=45&]\ Using no more than 7 of them in any instance,
//TableForm how many combinations can you form that will
Plotting in Mathematica...
Print["\n ",Length[%]," soln’s"] total no more than 45 pounds?
In[1]:= Plot[x^2,{x,-5,5}] Plot a function over an interval. If connecting
(* Approach No. 2 *)
Out[1]= -Graphics- from a Unix workstation or an X terminal, your We can solve that problem using Mathematica’s
solns=0;
In[2]:= Display["a.gif",%,"GIF"] graph will be shown in a new window; we also Combinatorica package, or we can just write a
Do[If[((25i+10j+5k<=45)&&\
Out[2]= -Graphics- show saving graphic output in gif format. little program to solve that problem directly by
(i+j+k<=7)),\
solns++, Null],\ looping through a three way nested do loop,
Note: besides GIF format, you can also use {i,0,1},{j,0,5},{k,0,7}]; using an if statement to tally only solutions
the Display function to save Mathematica Print["\n",solns," soln’s"] that meet the specified restriction.
graphics in PDF, EPS, PCL, PBM and
other formats. Mathematica on other platforms...

X has a site license for Mathematica covering

More Information About Mathematica...

The Mathematica Book, 4th Ed., by Stephen See also http://www.wolfram.com/ and
Wolfram [ISBN 0-521-64314-7, 1470 pages] is http://documents.wolfram.com/ for online
the definitive reference. copies of many Mathematica documents.