I wrote a small code to solve the equations. The problem is whenever i change the value of variable 'S' from 1 to 1.0 maple produces different solution. It means that i am only right solution for S=1. I can't get right solution other than S=1. I haven't faced this kind of problem before. Can anybody take a look? Is there something wrong with the code? Thanks!

restart:with(PDEtools):declare(f(eta),theta(eta),quiet):S:=1:
f(eta):=sum((q^i)*f[i](eta),i=0..16):
eqa:=simplify(diff(f(eta),eta,eta,eta)+(q*g*((f(eta)*diff(f(eta),eta,eta))-((diff(f(eta),eta))*(diff(f(eta),eta)))-(S*diff(f(eta),eta))-(((S*eta)/2)*diff(f(eta),eta,eta))))):
f[0](eta):=(alpha/2)*eta^2+eta:
for m from 1 to 8 do eq[m]:=coeff(eqa,q^m) end do:
for m from 1 to 8 do (simplify(int(int(int(((-1)*eq[m]),eta),eta),eta)+f[m](eta))):f[m](eta):=simplify(%) end do:
a1:=simplify(sum(f[n](eta),n=0..8)):
b1:=simplify(subs(eta=1,a1))=(S/2):b2:=simplify(subs(eta=1,diff(a1,eta,eta)))=0:
sys:={b1,b2}:solution:=fsolve(sys):s1:=solution[1]:s2:=solution[2]:
alpha:=rhs(s1);g:=rhs(s2);beta:=g^(1/2);
c1:=simplify(a1):c2:=diff(c1,eta,eta):c3:=subs(eta=0,c2)

Hello,

I noted that add and 'add' can give different result, can somebody help me understand it?  Consider the code below:

>phi1 := GAMMA(-(1/2)*vst-8.333500000-(1/2)*r)*GAMMA((1/2)*vst+21/2+(1/2)*r); phi2 := GAMMA(16.66700000+r-2*vst)*GAMMA(2.166500000+vst); xi1 := -vst; xi2 := 16.66700000+r-2*vst; z := 37.52950222; K := 9.846618489*10^(-38)*33.330^(.5000000000*r+10.)

### Then I have the following sums (they were supposed to give the same results):

>f := unapply(K*add(phi1*(-1)^vst*z^(-xi1)/factorial(vst)+phi2*(-1)^vst*z^(-xi2)/(.5*factorial(vst)), vst = 0 .. 10000), r):

>g := unapply(K*'add'(phi1*(-1)^vst*z^(-xi1)/factorial(vst)+phi2*(-1)^vst*z^(-xi2)/(.5*factorial(vst)), vst = 0 .. upto), [r, upto]);

# I am interested in the following quantities:

>evalf(f(1)/f(0));

-6.265883733

>evalf(g(1,10000)/g(0,10000))

25.80438437

That seems very strange to me, why maple gives different results?

Hi !

I want to solve several summations which at first glance seems very easy to me.

But I definitely hit a wall. I have Maple 2018 on Windows 10.

I'm not very familiar with the sumtools or Sumtools packages.

 I have attached a file that gives you an example summation and the

 solution I found and especially I show you the way that I obtained this solution.

  The big question is:

  Is there a way to transform this summation to arrive at my solution?

  Or . Can the result I obtained be demonstrated?

  Maple seems unable to solve this summation. Mathematica did no better.

  I know very well that this is more a question for mathematicians than for

Maple users. But when I look at the questions in Mapleprimes I see

although there are good mathematicians behind most of these questions.

Thank you in advance ! 

sum1.mw

Hi !

I just noticed a strange behavior of Maple with Greek letters. 
I have Maple 2017 and 2018 on Windows 10 .

Suppose you write a procedure that contains Greek letters.

If you run the procedure definition, the Greek letters are not displayed correctly. 
They will be displayed in words: alpha, beta, gamma,...
The procedure will still work.

Now close the file. Then open this file. Maple now shows you the procedure with the Greek letters correctly displayed.

I added a file as an example.
Maybe it's an annoying little problem known for a long time ? I don't know .
I have an old version of Maple V Release V (1997) which does not have this problem. 

Best regards !

alphabeta.mw

This might be a trivial question, but I have not been able to find the answer. I am using Maple 2018.

I am working with differential operators acting on real valued functions of a real variable. On the first hand, I want to be able to do simple algebra with these operators. For example, I might want to compute the commutator of two such operators. I have been using the DEtools package, toghether with the 'mult' command to perform such calculations. See the atached file containing a simplified version of what I am doing.

Now, I also want to be able to act with my differential operators on a function, and get the resulting function. 

(a) What command allows me to do that within my framework? (i.e. that of DEtools with the way I have defined and used my operators)

(b) I there a better way to proceed? (i.e. is there a better way to do both algebra with differential operators and to act with them on functions to get the resulting function)

Many thanks

Example.mw

H.M.mw

Download H.M.mw

Hi experts,

I have an equation. I tried to use subs command to substitute x with complex number p+iq in this equation and then separate the Re and Im part of the equation, but it does not work.

How can I do this?

How can I plot the Re and Im part of the equation vs y for some different parameters?

Equation:

eq.mw

I want to arrange this equation in term of powers of x and then plot the  real and imagenery part of x vs y. How can I do this with Maple?
1-alpha*((1/x^2)+(1/(x-y)^2)+(1/(x+ay)^2))=0;

1.mw     (alpha and a are constant, for example alpha=1 and a=0.3)

Hello. I want to organize some parameter
F=fsolve(1+x^2=j).

If j=-5, for example, an empty value 'F=' will be returned. How can I write the expression

If F=(void), then action 1, else action 2 end if

The fact is that I just don't know how to write the result (void), that is, the result of 'F ='

Thanks

Hello 

Imagine you have an irregular polygon as follows:

restart; with(plottools):with(plots):

display(polygon([[0,0],[200,0],[200,300],[250,400],[500,300],[500,500],[0,500]]),color=yellow,linestyle=solid,thickness=2);

with vertices that can be easily defined, and that you have an arbitrary point P(x,y) with coordinates of (x) and (y). I want Maple to check whether the point P(x,y) is inside or outside the polygon. In Matlab, the command (inpolygon) can be used for such a task. Is there a similar readily available command in Maple. 

Thank you very much for your help and support. 

From
with(numtheory):
msolve(x^2=-1,5);                                I get
                        {x = 2}, {x = 3}
But I need the integers 2 and 3. Who can help me?

Good day!

I would welcome any help to explain how to plot a simple function, z, of two variables, (x, y) so that I can show the relationships between variables.

Basically, I need to plot a family of z-curves for varying x and y values and also would like to understand how to display the z-y and z-x plots. 

Thanks for reading!

MaplePrimes_Nov_29.mw

Hello.

I wrote some codes for solving equations in a special way. My codes work very fine when my differential system has only one equation.

But when I extend it to a system that has two or more equations, sometimes does not work and I face with error. However, in one system (2 eq) it also works. But it is only a special case. If you want, I can upload it here as well.

I am attaching the maple file containing the error. You can observe it. 

By examining many cases, I think that maybe the root of the problem is that I must change my codes in a way that they are for a "system". Indeed, instead of writing codes for each differential equation, I must define a system and write these codes for a system containing some differential equations.
However, I think that the boundary conditions sometimes also may be the root of the error.

I change the attached codes for this purpose, but I failed. I prefer to not present my wrong codes for this purpose. 

If it is possible, help me to change my codes.

Thanks a lot.

ERROR.mw

I need help on maple code for solving both linear and non linear boudary condition for fractional order partial differential equation