Jun
4
A 10 minute tutorial for solving Math problems with Maxima
By Antonio Cangiano. Filed under Essential Math, Software
About 50,000 people read my article 3 awesome free Math programs. Chances are that at least some of them downloaded and installed Maxima. If you are one of them but are not acquainted with CAS (Computer Algebra System) software, Maxima may appear very complicated and difficult to use, even for the resolution of simple high school or calculus problems. This doesn’t have to be the case though, Maxima is very friendly and this 10 minute tutorial will get you started right away. Once you’ve got the first steps down, you can always look up the specific function that you need, or learn more from Maxima’s official manual. Alternatively, you can use the question mark followed by a string to obtain in-line documentation (e.g. ? integrate). This tutorial takes a practical approach, where simple examples are given to show you how to compute common tasks. Of course this is just the tip of the iceberg. Maxima is so much more than this, but scratching even just the surface should be enough to get you going. In the end you are only investing 10 minutes.
Maxima as a calculator
You can use Maxima as a fast and reliable calculator whose precision is arbitrary within the limits of your PC’s hardware. Maxima expects you to enter one or more commands and expressions separated by a semicolon character (;), just like you would do in many programming languages.
(%i1) 9+7;
(%o1)
(%i2) -17*19;
(%o2)
(%i3) 10/2;
(%o3)
Maxima allows you to refer to the latest result through the % character, and to any previous input or output by its respective prompted %i (input) or %o (output). For example:
(%i4) % - 10;
(%o4)
(%i5) %o1 * 3;
(%o5) 
For the sake of simplicity, from now on we will omit the numbered input and output prompts produced by Maxima’s console, and indicate the output with a => sign. When the numerator and denominator are both integers, a reduced fraction or an integer value is returned. These can be evaluated in floating point by using the float function (or bfloat for big floating point numbers):
8/2;
=>
8/2.0;
=>
2/6;
=>
float(1/3);
=>
1/3.0;
=>
26/4;
=>
float(26/4);
=>
As mentioned above, big numbers are not an issue:
13^26;
=>
13.0^26
=>
30!;
=>
float((7/3)^35);
=>
Constants and common functions
Here is a list of common constants in Maxima, which you should be aware of:
- %e - Euler’s Number
- %pi -

- %phi - the golden mean (
) - %i - the imaginary unit (
) - inf - real positive infinity (
) - minf - real minus infinity (
) - infinity - complex infinity
We can use some of these along with common functions:
sin(%pi/2) + cos(%pi/3);
=>
tan(%pi/3) * cot(%pi/3);
=>
float(sec(%pi/3) + csc(%pi/3));
=>
sqrt(81);
=>
log(%e);
=>
Defining functions and variables
Variables can be assigned through a colon ‘:’ and functions through ‘:=’. The following code shows how to use them:
a:7; b:8;
=>
=>
sqrt(a^2+b^2);
=>
f(x):= x^2 -x + 1;
=>
f(3);
=>
f(a);
=>
f(b);
=> 
Please note that Maxima only offers the natural logarithm function log. log10 is not available by default but you can define it yourself as shown below:
log10(x):= log(x)/log(10);
=>
log10(10)
=> 
Symbolic Calculations
factor enables us to find the prime factorization of a number:
factor(30!);
=>
We can also factor polynomials:
factor(x^2 + x -6);
=>
And expand them:
expand((x+3)^4);
=>
Simplify rational expressions:
ratsimp((x^2-1)/(x+1));
=>
And simplify trigonometric expressions:
trigsimp(2*cos(x)^2 + sin(x)^2);
=>
Similarly, we can expand trigonometric expressions:
trigexpand(sin(2*x)+cos(2*x));
=> 
Please note that Maxima won’t accept 2x as a product, it requires you to explicitly specify 2*x. If you wish to obtain the TeX representation of a given expression, you can use the tex function:
tex(%);
=> $$-\sin ^2x+2\,\cos x\,\sin x+\cos ^2x$$
Solving Equations and Systems
We can easily solve equations and systems of equations through the function solve:
solve(x^2-4,x);
=>
%[2]
=>
solve(x^3=1,x);
=>
trigsimp(solve([cos(x)^2-x=2-sin(x)^2], [x]));
=>
solve([x - 2*y = 14, x + 3*y = 9],[x,y]);
=>
2D and 3D Plotting
Maxima enables us to plot 2D and 3D graphics, and even multiple functions in the same chart. The functions plot2d and plot3d are quite straightforward as you can see below. The second (and in the case of plot3d, the third) parameter, is just the range of values for x (and y) that define what portion of the chart gets plotted.
plot2d(x^2-x+3,[x,-10,10]);

plot2d([x^2, x^3, x^4 -x +1] ,[x,-10,10]);

f(x,y):= sin(x) + cos(y);
plot3d(f(x,y), [x,-5,5], [y,-5,5]);

Limits
limit((1+1/x)^x,x,inf);
=> %
limit(sin(x)/x,x,0);
=>
limit(2*(x^2-4)/(x-2),x,2);
=>
limit(log(x),x,0,plus);
=>
limit(sqrt(-x)/x,x,0,minus);
=>
Differentiation
diff(sin(x), x);
=>
diff(x^x, x);
=>
We can calculate higher order derivatives by passing the order as an optional number to the diff function:
diff(tan(x), x, 4);
=>
Integration
Maxima offers several types of integration. To symbolically solve indefinite integrals use integrate:
integrate(1/x, x);
=>
For definite integration, just specify the limits of integrations as the two last parameters:
integrate(x+2/(x -3), x, 0,1);
=>
integrate(%e^(-x^2),x,minf,inf);
=>
If the function integrate is unable to calculate an integral, you can do a numerical approximation through one of the methods available (e.g. romberg):
romberg(cos(sin(x+1)), x, 0, 1);
=> 0.57591750059682
Sums and Products
sum and product are two functions for summation and product calculation. The simpsum option simplifies the sum whenever possible. Notice how the product can be use to define your own version of the factorial function as well.
sum(k, k, 1, n);
=>
sum(k, k, 1, n), simpsum;
=>
sum(1/k^4, k, 1, inf), simpsum;
=>
fact(n):=product(k, k, 1, n);
=>
fact(10);
=>
Series Expansions
Series expansions can be calculated through the taylor method (the last parameter specifies the depth), or through the method powerseries:
niceindices(powerseries(%e^x, x, 0));
=>
taylor(%e^x, x, 0, 5);
=>
The trunc method along with plot2d is used when taylor’s output needs to be plotted (to deal with the
in taylor’s output):
plot2d([trunc(%), %e^x], [x,-5,5]);

I hope you’ll find this useful and that it will help you get started with Maxima. CAS can be powerful tools and if you are willing to learn how to use them properly, you will soon discover that it was time well invested.
If you enjoyed this post, then make sure you subscribe to our RSS Feed.
Comments
60 Responses to “A 10 minute tutorial for solving Math problems with Maxima”
Leave a Reply
Subscribe Articles
Subscribe Comments
Subscribe via Email
This is awesome. Awesome! Thank you.
– Josh
Cool beans. Very nice write-up.
nice infos, in my opinion these kind of software should be much more supported by universities, instead of (implicitly) encouraging the use of pirated software…
thanks! I was just trying to improve my Maxima skills…
ciao ciao!
Thank you. Very useful intro.
Wow. Looks awesome. The thing that REALLY sealed the deal was the TeX output. Love it.
[…] A 10 minute tutorial for solving Math problems with Maxima (tags: math tutorial) […]
[…] A 10 minute tutorial for solving Math problems with Maxima (tags: math lisp tutorial) […]
[…] Link Original: A 10 minute tutorial for solving Math problems with Maxima […]
Great software. Great for my high school subject.
this was NOT a guide for high school calc…
it’s a great, vague, manual snapshot
..
a guide would have a scan from a high school calc book, and STEPS with screenshots
A 10 minute tutorial for solving Math problems with Maxima(10 分鐘學會使用Maxima)
之前那篇
3 Awesome Math Programs(3個超棒的數學繪圖程式)
裡面有提到這個程式Maxima
官方的教學文件在這裡
剛剛在Digg上看到一篇好文章
圖片和解釋都說得很詳盡(想當初用Mathematica用到快起肖= =R…
man, where was all this cool free math software when I was stuck using maple on mac plus? argh!
probably a good thing that it “appears” difficult as it’d a shame for highschoolers to simply use it instead of learning basic algebra.
This is awesome little tool. Wish I knew about this in my school.
[…] read more | digg story […]
I have been using Maxima as I re-learn Algebra. I use it it verify what I do by hand, but sometimes I would like to see some of the intermediate steps. Is there any way to do that?
How the hell do you translate these formulas to be used in software programs like javascript and actionscript?
Could you possibly create a tutorial on what all the symbols mean?
[…] Read the original source… […]
Excellent write-up. How did you get the LaTeX output on the fly? The presentation looks great!
The fancy calculus stuff doesn’t translate. Most of the basic functions (sin, sqrt, log, etc) are available as Math.sin() and whatnot. Google for “javascript math” to get a full list of Javascript’s math functions and constants.
@rwinston
Thank you very much. I used a LaTeX plugin for Wordpress, and I wrote the TeX code for simple cases and Maxima’s tex(expr) function for the longer ones.
great job ! do u know of a similar tutorial for R ?
@easan
Thanks a lot easan. For R, take a look at this tutorial: http://www.cyclismo.org/tutorial/R/
This is really cool, sadly the school year is now over and I can’t use it!!! fdasdafhasdkj! GRRR!
Maxima has its roots in Macsyma which existed before Mathematica and Maple.
See the following article about Macsyma:
http://en.wikipedia.org/wiki/Macsyma
I’ve recommended Maxima to many people. Most of these people rejected it and did not even bother to look it up because they never heard of the name or believe free software cannot be of good quality.
I would like to see a side-to-side comparison between Maxima and the commercial alternatives. I can imagine Mathematica and Maple having more features (never missed them though) and a more optimized solver.
[…] A 10 minute tutorial for solving Math problems with Maxima (tags: mathematics opensource tutorial) […]
great program! helps a lot at the university
[…] Link: Math Blog […]
Thanks!
I’d messed with Maxima a few months ago, and found it a painful experience trying to figure it out from the manual. Went thru your tutorial last night and now I can use this for solving real problems.
This will be a HUGE asset for me.
Alex
[…] http://math-blog.com/2007/06/04/a-10-minute-tutorial-for-solving-math-problems-with-maxima/ […]
Hi!The information was delivered in a simple way that it is very easy to understand in using the software which initially i found it very difficult to use Maxima when I explore it al by my own. But now you had help me to solve my probs.
It is indeed a fantastic piece of info for me especially at this moment in rushing my assignment. Thanks…
Got exactly what I wanted.
Thank you so much@
Hi
Thanks for the Tut. However I cant find information to help me. My problem is as follows:
I want to solve f in the following two functions, when I choose the following:
s=2400 and d = 2000. f is the required value, which is the focal point of a parabole.
h=d^2/(16*f);
s=d/2*(sqrt((4*h/d)^2+1))+2*f*ln((4*h/d)+sqrt((4*h/d)^2+1));
I substitude h in the second function and then evaluate for the said values of s and d.
All goes well and is calculated. However I then can not get a calculated value for f from Maxima. The following eauation is given.
2400=2*ln(500/f+sqrt(250000/f^2+1))*f+1000*sqrt(250000/f^2+1)
What can I do to get a calculated value for f?
Thanks
Jacobus
with Mathematica, I can do wonderful things to turn a series of calculations into a semi-professional looking document. Can I do that with maxima? (or wmaxima, or …?) I’ve played with it (on windows) and it just seems like a ‘calculator’ I can solve it all, but it’s not really made to give a nice formatted output. Seems. I claim no full knowledge, maybe I’m missing something.
Thank you. Your article is very useful to me to start with with Maxima.
Rajendra
JiK, you can use Maxima from within TeXmacs or from a SAGE notebook.
[…] read more | digg story […]
i think that knowing all this math will be very helpful in the long run. it should really help you in the future
Any suggestions on how to get Maxima to generate a surface of revolution; e.g., the surface formed between y=x^2, y=0 and x=6 is revolved about either axis?
Thanks.
Very usefull. Muito bom mesmo.
fantastic program! great work guys!
After a short tutorial, I tried a real problem:
solve([x*(1+(0.083/y)^(1/z))-41, x*(1+(0.833/y)^(1/z))-44.9, x*(1+(8.333/y)^(1/z))
-47.8], [x,y,z]);
____
But get no response, anything I was wrong?
Hr.stein, I run your system of equations and as a response I got that it has no solutions ([]).
thanks for your reply, Antonio!
In fact, I would expect an approximate solution using least square approach, with indication of deviation.( somewhat like the Mathematica does).
I still think Maxima could solve the equations, only I miss some built-in functions or there is mistake in my expression.
Regards!
log (3m+7)-log (m+4)= 2log 6 - 3log 3
6 6 6 6
This is a fantastic tutorial. Thanks to the author, I was able to discover the great power of this amazing program. Based on this page, I started to write a tutorial for Maxima in Greek.
http://lomik.wordpress.com/maxima/
Hello
How can I tell maxima that a given variable is real and greater than zero in the expression below?
integrate((a/pi)/(a^2+x^2),x,-inf,inf)
fantastic!! This is a great getting started tutorial and showes the power of the program.
I teach at University of Colorado and really push the value and quality of open source software (which I use in my own research). Next time I teach numerical analysis I will definitly use Maxima!
Excellent!! Thanks for your great tutorial for Maxima.
Wow, just from reading about the Maxima wants me to use it, it seems that its a very powerful mathematical tool
I used to lust after Mathematica but it is outrageously expensive: $2,500. We are approaching the point where open source alternatives, like Maxima, are the way to go for heavy duty computer assisted symbolic and numerical mathematics.
[…] Essendo studente di ingegneria, mi trovo spesso nella necessità di utilizzare un CAS (Computer Algebra System) ovvero un programma in grado di eseguire calcoli numerici, simbolici, grafici e altre operazioni correlate. Il software più diffuso per effettuare questo tipo di operazioni è Mathematica il quale però è un software commerciale, così mi sono messo alla ricerca di un programma open source con funzionalità simili. Tale programma esiste e si chiama Maxima. Con esso è possibile svolgere tutta una serie di operazioni come derivazione, integrazione, risoluzione di sistemi di equazioni, grafici di funzioni in 2 e 3 dimensioni e molto altro ancora. Per facilitare l’utilizzo del programma è disponibile wxMaxima un comodo front-end grafico, che permette rapido accesso alle principali funzioni del programma. Di seguito trovate alcuni screenshots con qualche esempio. Il testo degli esempi è tratto dall’ottima guida: A 10 minute tutorial for solving math problems with Maxima. […]
I like the way you are presenting this material. I have used Maxima and Scilab in teaching a course on Computation. I am a big fan of Maxima. Scilab is OK but crashes a lot. ON the other hand Maxima is extremely solid with many options for the user interface. It is promoted as a symbolic math program but it can do most of the run of the mill numerical math. Here is the link to the course material I taught.
http://www.eng.ysu.edu/~jalam/engr6924s07/
Feel free to take any part of the material if you want to write more articles about math calculation using Maxima.Also any feed back from is highly welcomed.
Very nice tutorial, Antonio. I’m signing up for your RSS feed.
Tell me, is there a way to zoom and rotate the 3d graphs? Also, my graphs are not solid planes like your screenshots but rather grids. How does one make the planes solid? Thanks.
Richard,
Use assume(a>0);
Dotan,
To make the graphs coloured and solid you could add [gnuplot_pm3d,true] to the end of the statement eg:
plot3d(f(x,y), [x,-5,5], [y,-5,5],[gnuplot_pm3d,true]);
Re: Jacobus
Better late than never?
(%i1) find_root(subst( [h=d^2/(16*f),d=2000,s=2400], s=d/2*(sqrt((4*h/d)^2+1))+2*f*log((4*h/d)+sqrt((4*h/d)^2+1))),f,1,2400);
(%o1) 421.811063853618
I’ve 2 issues here:
FIRST: Halting problem in TM
——
The below symbolic example doesn’t terminate after many minutes in maxima-5.14.0 although the human solution is f(x)-g(x)-h(x):
factor(expand((f(x)-g(x)-h(x))^500));
To fix this problem in the future: better
artificial sustitution in the solving of
factor or intelligent lazy computation as
factor(expand(Alpha)) => Alpha.
SECOND: Non-deterministic Memoization
——-
With
factor(expand((sqrt(x+a)+b)^10));
i got once the imprecise solution
10 b x^4 sqrt(x+a) + 120 b^3 x^3 sqrt(x+a) + … + 45 a^4 b^2 + a^5
but i got many times sqrt(x+a)+b)^10 .
I think that the maxima system
uses memoization that remember some
internal answers of previous computations
but it’s non-deterministic.
axima 5.14.0 (with CLISP 2.41) has 2 symbolic bugs:
* First bug: it doesn’t factor when it uses sqrt instead sqrt2 :
(%i1) factor(expand((sqrt(x)+y)^2));
2
(%o1) y + 2 sqrt(x) y + x
(%i2) factor(expand((sqrt2(x)+y)^2));
2
(%o2) (y + sqrt2(x))
* Second bug: wrong result due to the sign: (y-x)^1000 instead of (x-y)^1000 :
(%i3) factor(expand((x-y)^1000));
1000
(%o3) (y - x)
Is there any way to get nice looking latex fonts in the EuMathT notebook? For example, if one calculates >:integrate(x,x) one gets
x^2
—–
2
and this in not nice looking. How to get LaTeX fonts?
Seem like a good maths tool. Students of maths should welcome this tool. Nowadays, learners go for effective learning tools. This is a good target.