Hi all,
I solved this nonlinear DE
de := diff(y(r), r, r)+2*(diff(y(r), r))/r+9*(16*43)*Pi^2*sqrt(Pi/(2*(1/43)^3))*polylog(3/2, -exp(1/43*(43-y(r))))/(16*Pi^2*sqrt(43)^3) = 0;
with initial conditions y(0) = 0, D(y))(0) = 0
by using truncated series method in 'solve symbolically' option (Dsolve[interactive]).
I am getting a solution like this:
y(r) = -(3/4)*sqrt(86)*sqrt(Pi)*polylog(3/2, -exp(1))*sqrt(43)*r^2-(1161/80)*Pi*polylog(1/2, -exp(1))*polylog(3/2, -exp(1))*r^4+O(r^6).
How can I get an asymptotic potential y(r) for r-> infinity?
Thank you
MS

Dear All,

i want to used optimizatiom method in Maple ( Paritcal Swarm optimization) , i have the algorithm but i can to write the codes in Maple, vuold any one help for that.

Thanks and best regards.

Ahmed

Hi all,
I have this same equation in another active thread also with boundary conditions. But now I have INITIAL CONDITIONS INSTEAD. I thought of just starting a fresh post. Here is my non-linear differential equation
x := diff(y(r), r, r)+2*(diff(y(r), r))/r+ C*polylog(3/2, -exp(C2*(C3-y(r)))) = 0
with boundary values D(y)(0)=0 and y(0)=0. C, C1 and C2 are some numbers.
The hint says to use procedure to solve polylog and then integrate the differential equation.
I was wondering if I can get some advice on solving this equation.

Isn't there a way to get Maple to document each step in its execution of a command? E,g,, I have a complicated expression that I want to SIMPLIFY & I would like to see how Maple does it.

Alla

I want to evaluate an integral over a numerically calculated function. Specifically:

Hello everybody,

Im trying to plot a f (x) function but maple keeps giving me an empty plot error.

From my maple worksheet:

**> f:=x-> (7+4*sin((2*pi/43200)*x+3*(pi/2)));**

**>plot({f(x),8},x=-1000..360000,y=-1..11);**

Error:

I recently realized that maple/standard has some plotting capabilities not available in maple/classic. I nevertheless plan to stick with classic except for the occasional printing of graphs. I don't want to write multiple worksheets, so I want my worksheet to know when I'm printing in standard and when I'm printing in classic. I don't want a maple/classic plot to overwrite a maple/standard plot (and vice versa). I want to set some classic-specific and standard-specific options.

I have an equation of a line that intersects the ellipse. I can get Maple to find one intersection. How do I get the other?

E1:=x^2/300+y^2/1000=1;

E2:=y=2*x+40;

solve({E1,E2},{x,y});

evalf(%);

gives me {x=-4.810723698, y=30.37855260}

or I could use fsolve({E1,E2},{x,y}); and get the numeric result right away.

There is a second intersection point I can find graphically around x=-5 y=30 but I don't know how to get maple to find it numerically.

Dear Sir

Can you please help me to **minimize** the following cost function with maple 10.I want that **" t and T must be greater than zero(not equal to zero) , t < T and C > 0." **I tried a lot but I am not getting the values of t and T as mentioned above **" ". **Also why the solution obtained by the 2 methods below is different ? Which one of them is correct ?

Priyaasha

The birthday problem is a statistical problem where a randomly chosen set of people will have the same birthday.

The calculated probability is actually lower than reality. Why? Most people tend to "celebrate" on special days. The weekend of Canada day, Valentines Day, Victoria day, etc ... and so you will find higher groupings of people with birthdays 9 months following these celebration days. Hence probability is actually a little higher that 2 people will have the same birthday.

Hi all,
I have a non-linear differential equation
x := diff(y(r), r, r)+2*(diff(y(r), r))/r+ C*polylog(3/2, -exp(C2*(C3-y(r)))) = 0
with boundary values D(y)(0)=0 and y(infinity)=0.
The hint says to use a subroutine to solve polylog and then integrate the differential equation with help of a Runge-Kutta like scheme.
I was wondering if anyone could help me with this. I have no clue how to use subroutines and specially when y(r) is within polylog.
Thanks in advance,
MS

Is this the best way to produce a random 2-d array?

L:= Array(1..10,1..10,(rand(1..50));

What's unusual about the operation is that I seem to get an unusually high number of doubles in the rows. All of the time I'm getting 50% of the rows with a matching pair of numbers (ie. every time I run the line, I get 5 rows with a matching pair of numbers)

Does the probability of that seem a little high? Also are there other ways to produce a 2-d random array?

Thanks.

Hey, I have this complex ecuation on maple: (x+Iy)^3-x-Iy-A=0,

And I have to do an animation for -2<=A<=2, I tried using implicit plot and giving A that range, but it only returns a blank plot, any help would be appreciated!

Oh! I also have to search for what A values the ecuation has 1, 2, and 3 real solutions, any help in this would be appreciated to!

Thanks,

SolracQ

These functions give different indices for the x variables:

> test := j->sum(x[j-i mod 24], i = 0 .. 8); test2 := j->seq(x[j-i mod 24)], i = 0 .. 8);

> test(3); test2(3);

x[3]+x[26]+x[49]+x[72]+x[95]+x[118]+x[141]+x[164]+x[187]

x[3], x[2], x[1], x[0], x[23], x[22], x[21], x[20], x[19]

The second line uses the indices I would like. What am I doing wrong?

Hello!

Is there a function that allows finding address of array/matrix's element equal or closest to given number. Or I have to write a loop.

en example

array

[1 2 5 6 9 10]

given value

5.1

the closest is 5 and it is in 3th column.

Regards,

Smagliczka