I am looking for help with basic algorithms to approximate numerical solution to the N body problem.
I have been trying to develop basic algorithms using Runge Kutta, Euler ect. Initially I specified to the 3 body case and the equations of motion. I am now trying to develop a more general algorithm. If anyone has knowledge in this are or can look through my annotated procedure to see where I have gone wrong I would appreciate it.
I am very new to maple and would appreciate detailed exposition, "for dummy's", if you will.
Thanks in advance

Maple should do this

1e123456789876 + 1e123456789876

float(infinity)

The answer is 2e123456789876.

Is there a way to get maple to actually do the calculation without doing any manipulation of the numbers before hand? Maybe there is a package out there for working with extremely large numbers?

How do I do to draw a graph from an adjacency matrix? Notice: the matrix is symmetric and when I use GraphTheory package, occurs the error: "Error, (in GraphTheory:-Graph) an undirected graph cannot contain loops (loop detected at vertex 1)".

I've tried the "graphplot3d" command, but I need to know the in-degrees of each vertex (it´s a rotine). It's impossible, in this case, using the "InDegree" command.

Help me, please! Regards,

> PDE := diff(u(x, t), t) = k*(diff(u(x, t), x))+diff(u(x, t), x, x);
> a := 1/70;
> IBC := {u(0, t) = 0, u(5, t) = 0, u(x, 0) = exp(-(x-2.5)^2/(2*a^2))/((2*Pi)^(1/2)*a)};
> pds := pdsolve(PDE, IBC, numeric, time = t, range = 0 .. 5, timestep = 1/2000, spacestep = 1/2000);
>R := pds:-value(output = listprocedure);
>uu := subs(R, u(x, t))

Hi,all,
here is my question:
alias(alpha=RootOf(x^32+x^26+x^19+x^15+x^13+x^11+x^9+x^8+x^4+x+1));
and then i want to compute:
alpha^12345678987654321;
But when I did it in maple, it failed.
I think there should be some way to do it.
Would you like to help me?
Thanks
Gepo

Is there a problem (strange) with the new algolib, or maybe an existing problem elsewhere that shows up when algolib is present? Is it only under MacOS, or more general?

The new version of algolib has been out about 2 weeks now. It is here... http://algo.inria.fr/libraries/

I downloaded files algolib.mla and algolib.hdb, but them in a directory "/Applications/Maple/algolib". It seems to work fine. But today I had a strange problem with latex(), and...

Hello,

I am wrting my thesis which is on some MapleSim models. In Ref I introduce MapleSim, I wanna know if MapleSim has Document Object Identifier (DOI) number or name that I can refer to.

Thanks so much,

Bye

I checked some formulas and realized, that my Maple version simply has no representation for some Greek symbols e.g. delta or epsilon, though it can print alpha, beta, gamma, xi omega. So if I type 'delta*alpha', Maple just prints 'alpha' though it calculates 'delta*alpha'.

Does anyone of you have any idea why this is the case, or even better, how I can change that?

Btw.: My version is 'Maple 12.0'.

Thanks.

Hi there!

Suppose I have an expression

Expr:=a*f(x[0],...) + b*f(x[1],...) + c*f(x[2],...),

where the ellipses represent (different) functions of x[0]..x[2]. I'd like to use frontend to get the coefficients a, b, c corresponding to x[0], x[1] and x[2]. If the function f only has one argument, that's peanuts:

Hi,

How I do to run the commands of the package in Maple 13?

Gracias

For some reason Maple just will not solve this...

2.160417096*x^(2.513) + 0.9921351694*x^(0.512) -15.95645167 = 0

It runs for about 10 minutes before I decide to give up.

The only thing I can think of is I had one of those 'lost kernel' messgaes the other day and I closed Maple while it was trying to restart the kernel so maybe that caused some problems with my system...

I want to factor the 34 th Fermat number.

F(n) := 2^(2^n) + 1

I can enter the 29th Fermat number: 2^(2^29) +1, and Maple will evaluate it.

I seem to be under the maxdigits limit:

Hello Friends;

How could make maple show me the steps and rules applied in the following exerted?

3^(-2*x^2)*9^(-2*x^2-3) = (1/27)^(6*x-2);

Something very similar to the command:

Student[Calculus1][ShowSolution]

Thanks

I'm getting started with this. I'm reading books about General Relativity. Their notation uses superscripts and subscripts in Greek, for examples. But in the lessons for DifferentialGeometry and in applications I find in the Maple website, I see coordinate systems of (x,y), for example. Indices, raising and lowering, and other operations involve "d" and "D_" notation. Examples use something like T12 to represent what I would want to see as "T" with two subscripts having values of 1 and 2.
I apologize if I'm not explaining the background to my question well.

Hi guys,

I am not a experienced user of Maple. So the question is probably very basic. Any help is appreciated.

I have 2 coupled first order linear equations. I may be able to solve them manually, but I thought it would be a perfect application for Maple. The problem is that I have no idea how to do it in maple. So here it goes.

d(y1)/dt = (a-y1)/T1 - k1*y1 + k2*y2;

d(y2)/dt = (b-y2)/T2 - k2*y2 + k1*y1;

Under the following constraint: a*k1 = b*k2;

How may I solve this problem in Maple?

Thanks a lot.