MaplePrimes Questions

how to invert power series beginning with a square term?

I'm sorry if I'm duplicating the information here. I tried to post a comment but it doesn't seem to be appearing.

 

When I change = to := it does not work. Maple says

 Error, invalid input: MultiSeries:-series expects its 2nd argument, eqn, to be of type {name, name = algebraic}, but received 1.*I

I'm using the pdesolve function for a linear pde system of two equations in x and y, with indeterminate functions V(x,y), W(x,y). I'm trying to incorporate Neumann boundary conditions such as:
eval(V[x],x=0) =1;

but when I do so in the pdsolve function I get the error message that the bcs are nonlinear in V and W. Is there a different way that I should be declaring these bcs for use in pdsolve? Or should I in fact be trying to incorporate the bcs some other way?

Hello,

 

I have a huge expression involving psi and phi called A_1 which I paste below.  I am trying to integrate it using the residue theorem in the variable phi. I make a substitution in u as below for sin(k*phi)  and ask for the singularities of each operand at u. But when I get my final result it still has the variable u in it. If its evaluated the integral in u then it should not return u in the result. Evidently it’s done something...

I have a system of 3 ODEs in 3 variables x(t),y(t),z(t) to be solved numerically. I would like to set up events to halt integration as soon as one of the derivatives changes sign.

I've tried things like:

'events' = [ [ diff(x(t),t) >0, 'halt'] ] :

'events' = [ [ diff(x(t),t) = 0 .. 0.1, 'halt'] ] :

I might have tried something like increasing(x(t)) (how nice and intuitive), but if I understand correctly the help page on this, increasing(x(t)...

Hi,

 

I have been defining a function:

 

f1:=(m,r0,en)->evalf(Int(....),y=0..infinity)

 

for some reason when I call this function and pass it values say f1(0.1,2,3) Maple is unable to return me an answer, it just spits back the unevaluated inert integral with the values of the parameters inserted.

 

However, I know that Maple can do this integral if I just input :

 

Really enjoying Maple 12 in document mode, but am finding the results to be unreliable sometimes. It doesnt seem like I can attach the maple file here (which is annoying), but I often find that when, for example, solving for a variable, the results are not shown properly as an equation, but rather include text which looks as though maple is showing the equation with some kind of internal code, for example d ...

Hi guys I am trying to evaluate  an integral of the form

 

A_0 := (1/(2*Pi)+Q*cos(k*phi))^2/(1+(1/(2*Pi)+Q*cos(k*phi))^2);

In this case I specify a value for Q=15 as given below and since maple cannot evaluate the integral directly I have written code to evaluate the integral using the residue theorem. So for each singularity I calculate the residue at that point and multiply the sum of the residues by 2*Pi*I as in the code below. However the...

Hi,

I have a problem with the solve comment. I have two equations with two variable. I try to use solve or fsolve but there is a message appears as follow

" Warning, solution may have been lost"

Thus could please help me ho to find the solutions?

 

Best regards

Hi everyone, 

 

How I can split loops than in openMP library 

 

#pragma omp for

 

using threads in maple ... please a one example

 

thanks by your attention.

If I have 

F = x + 2π_Z1~

Is there a function/operator I can apply to F which returns only x, discarding 2π_Z1~?  That is, I want to return only the base angle, and dump the added integer multiple of 2π term.

Consider a function f depending on a parameter k and a variable y

f:=(k,y)->k*s(y);

where s:=y->add(sin(1/y^i),i=1..1000000) is just some nasty function.

Now, why does it take a long time to evaluate f(0,y)? And how could I speed this up?

 

More generally, I would like to plug (varying) parameters into a function, simplify and then do analysis on the resulting function. How could I do that?

for instance: s1:=implicitplot3d(z=x*y, x=-8..8,y=-8..8,z=-8..8,axes=normal): s2:=implicitplot(x^2+y^2=1, x=-8..8,y=-8..8,axes=normal): if I type display(s1,s2); it gives an error.
Pi^2*(c+(1+Pi)/Pi)^2+Pi^2)/((c-(-1+Pi)/Pi)^2*Pi^2+Pi^2);

                         2 /    1 + Pi\2     2
                       Pi  |c + ------|  + Pi
     ...

I have a equation which I want to solve with HPM method, it works this way that you put a series instead of the dependent variable, and make infinite linear differential equation out of your nonlinear, then after solving some of them that meet your accuracy and sum them, that makes the answer. Now I have many problems with this which I hope you may help me.

1. How can I programatically asssign result of a dsolve to the variable? for example when I have dsolve(diff(f(x...

First 1828 1829 1830 1831 1832 1833 1834 Last Page 1830 of 2430