MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

How to solve the linear system by using maple

2y' + xy =2  with y(0)=1


I want to solve the following equation:

equ := 4*arctan(x)-4*x/(1+x^2) = 0

so I use fsolve(equ,x) to get all solutions of equ, but only get one of the roots:zero. Actually, equ has three roots, as the picture which can be generated by the following commands illustrated.

f1 := arctan(x):

f2 := 4*x/(x^2+1):

p1 := plot(f1, legend = ["arctan(x)"], color = red):

p2 := plot(f2, legend = ["4x/(x^2+1)"], color = blue):

as you know the number of solutions for order 3 to 5 magic squares is as follow : 
(not counting rotations and reflections)

order 3: 1 solution

order 4: 880 solutions

order 5: 275305224 solutions

higher order: unknown ...


in the first look at 880 , simply appeared it can constructed by
2x5x8x11 (every factor increased by 3)

Just a question about Maple's formatting.  When I type parenthesis, brackets, or absolute value bars, how do I control what size they are displayed?  I'm typing a bunch of expressions as in-line comments for an assignment, and there are multiple parts with parenthesis within brackets, etc.  On certain lines, ever time I put a space after typing, the size of brackets on a previous line will alternate between large and small.  I can insert spaces into random...

I cannot activate Maple 11 on my new MacBook.  I get the error message: Not enough activations left on this account.   As best I can recall, I installed Maple 11 twice before: on my desktop and on my laptop; the laptop has since crashed.

Might you be able to help me? 



I have the equation to calculate the legendre functions , but i am trying to find a way to form a loop so it would be less time consuming to find for a series of values.

Pn = sum((-1)^i*factorial(2*n-2*i)*x^(n-2*i)/(2^n*factorial(i)*factorial(n-i)*factorial(n-2*i)), i = 0 .. (1/2)*n)


This is the general equation , as you can see there are different values of P like p0,p1,p2,p3...depending...


I am rather new to maple and have a few questions considering one of my programs.

1. In my program I use 9 nested for loops. Despite trying to figure something out, I couldnt come up with a way of using less for loops.

2. In some of the for loops I want to start the loop from a variable of the for loop above it. Maple doesnt seem to be very happy about this.

3. I want to make a list of all the results produced by the for loops. However...

Hi- I'm having trouble figuring out how to combine two separate DEplot graphs into one graph that represents the SUM of the two.  Not two separate DEPlot graphs in one figure, but the sum of two DEplot graphs in one figure.  I'm using Maple 14- anyone know how to do this?  It's probably simple, and I'm just missing it.  Thanks in advance.

Having overcome a Java VM problem, I have a further problem which is now
limiting my use of MAPLE 14 on Vista:  plotting is limited to plots with
relatively few points or graphical content.

Here is an example in MAPLE 14 (in classic and the normal modes):

> with(RandomTools):with(Statistics):
> N:=400;
> rlist:=(Generate(list(float(range=-1..1),N))):
> plot(rlist,color=black,thickness=1);

This works fine up to N=170,...

m1:= [[1, 2, 3], [ 4, 5, 6], [7, 8, 9]]:

m2:= [[2, 2, 3], [ 4, 5, 6], [7, 8, 9]]:


Given my commands (above), I would like to save the output of




in the directory "c:\users\desktop\work"


and the ouput of


as m2.tex

in the directory "c:\users\desktop\work",

How could I do that?


I would like to be able to extract data from a plot strucutre. So say I've done:

p:=plot([sin(x), cos(x)],x=0..1, legend=["sin","cos"]):


now I can see the raw strucuture with


but not quite sure how to extract two lists of tuples that contain the data. I would also like to extract the two items in the legend.

Does anyone know how to do that??


eqn1:=(1-x)*exp(c*d*M*z + a*(1-d)*M*y+(M/2))-x*exp(-c*d*M*z - a*(1-d)*M*y - (M/2))=0;

eqn2:=(1-y)*exp(a*d*M*x + b*(1-d)*M*z+(M/2))- y*exp(-a*d*M*x - b*(1-d)*M*z - (M/2))=0;

eqn3:=(1-z)*exp(b*d*M*y + c*(1-d)*M*x+(M/2))- z*exp(-b*d*M*y - c*(1-d)*M*x - (M/2))=0;


When I try to solve this system of equations as:

solve( {eqn1,eqn2,eqn3},{x,y,z} );


solve( {eqn1,eqn2,eqn3,0<x,x<1,0<y,y<1,0<z,z<...


I am trying to be able to solve my plots into a file programatically (not via right click/export).

when I do this:


plotsetup(default):  #optional - just to be sure start off at default.
plotsetup(ps, plotoutput="test.eps", plotoptions=`portrait,noborder`);

things work ok

BUT.. doint this:

Hello all,

I'm sorry that I'm just not very good at working with Maple's optimization tools, but I appreciate your willingness to work with me. So, this is all related to the same problem I had been working with a few months ago. I used the "operator form" of NLPSolve to minimize the squared residuals of a temperature distribution with an external dataset by changing a few...


I do want to export as  LATEX input (displaying only 3 digits) my matrix "tableau".

After using

latex(  evalf(tableau),   3 )

I get too many decimals. What should I do to send as latex input only 3 decimals. My commands are below and my worksheet attached.

Thanks in advance for any help,





f1 :=  (w1, w2,n) -> add(add(

First 777 778 779 780 781 782 783 Last Page 779 of 1304