from help, it says 

The expression assuming property calling sequence evaluates the expression under the assumption property on all names in expression.

Then why 

restart;
expr := Sum((-1)^n - 1, n = 1 .. infinity):
simplify(op(1,expr)) assuming n::even;

does not simplify expr to zero, while

restart;
expr := Sum((-1)^n - 1, n = 1 .. infinity):
assume(n::even):
simplify(op(1,expr))

does simplify expr to zero.

I would have expected both to give zero. 

If I have a tensor  in 3-dimensions which is symmetric on  and anti-symmetric on , then the number of independent components should be zero. However, if I put this into Maple, using  followed by  it returns .

Any insight into why it does this would be great, thanks.

I am trying to customize the Venn Diagram generated by the VennDiagram command in the Statistics package. I would like to be able to shade just any of the regions. For example, if only two sets (A and B) are being displayed, I would like to be able to shade just A, just B, both A and B, or both A and/or B. Also, is it possible to have a rectangle surrounding the circles which represents the universal set? Could the circles be labeled with the letter "A" or "B" to be able to identify each set?

venn.mw

Here is a strange one...

1> /Library/Frameworks/Maple.framework/Versions/2019/bin/maple ; exit;
    |\^/|     Maple 2019 (APPLE UNIVERSAL OSX)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2019
 \  MAPLE  /  All rights reserved. Maple is a trademark of
 <____ ____>  Waterloo Maple Inc.
      |       Type ? for help.
> version();
 User Interface: 1435526
         Kernel: 1435526
        Library: 1435526
                                    1435526

> DE := (28*x + 44)*u(x) + (336*x^2 +
> 726*x - 12)*diff(u(x), x) + (144*x^3 + 396*x^2 - 9*x)*diff(u(x), x, x);
                               2               /d      \
DE := (28 x + 44) u(x) + (336 x  + 726 x - 12) |-- u(x)|
                                               \dx     /

                               / 2      \
             3        2        |d       |
     + (144 x  + 396 x  - 9 x) |--- u(x)|
                               |  2     |
                               \dx      /

> dsolve({DE,u(0)=2},u(x));
memory used=21.5MB, alloc=44.3MB, time=0.37
memory used=53.3MB, alloc=84.3MB, time=0.94
Error, (in dsolve) when calling 'property/ConvertRelation'. Received: 'numeric
exception: division by zero'

Presumably, the solution should be

u0:=2*HeunG((11 - 5*sqrt(5))/(11 + 5*sqrt(5)), 352/(9*(11 + 5*sqrt(5))^3*(-11 + 5*sqrt(5))^2), 1/6, 7/6, 4/3, 1/2, -8*x/(11 + 5*sqrt(5)));

(I get that by replacing coefficient 44 in DE with variable e44, solve, then substitute back e44 = 44.)

But maybe the problem is that this solution turns out to be an algebraic function:

u1:=2^(7/6)/(1 - 22*x + sqrt(-16*x^2 - 44*x + 1))^(1/6);

Dear friends, please I would like to ask for your help with an odd problem I have using the remove command. 

I have an array 

A:=Array([1,4,1,7]);

and I need to remove its first element A[1]. 

A:= remove[flatten](x -> x = A[1], A);

Instead of getting the result  A:= [4 1 7] I get  A:=[4 7], and I can't understand why. 

Could you please help me with a solution to the problem? Many thanks for the help.  

Hi,

I have been  working on a Maple code written almost 17 years back. The code generates a 3D model input for modelling analysis in ABAQUS software. The input file generated is different from the usual ABAQUS input files. I am seeking help with how to make changes to the input file generated to import into the  ABAQUS. I am uploading the notepad version of the input file as the Maple does not allow .inp files. Please do find the below attachment.

wucell.txt

How can I find the corresponding group for the Lie algebra given in the picture (using Maplesoft software)?

Also, the command Lies Third Theorem works only for Solvable representations. What to for unsolvable representation.

Given two sets of lie algebra data, How to check, using maple software, that these lie algebras are isomorphic?

for example : 

The two sets of lie algebras are given as : L1 := [[e1, e4] = e1, [e2, e3] = e1, [e2, e4] = e2]

and

L2 := [[e1, e2] = e1].

Hi there,

Could you help me with Harley's norm computation algorithm that is based on the Fast Extended Euclidean Algorithm that was suggested by Harley in an email to NMBRTHRY list in 2002 and that described in Vercauteren's thesis pp 87-90:

https://pdfs.semanticscholar.org/c945/c98267db064b272c87a885fc5eeb764b0b2d.pdf

My implementation working correctly and fast for low degree polynomials without modulo and for high degree polynomials with modulo M, where M is a prime number greater than 2^N. But all I need - it's a resultant modulo 2^N (or 2^(Nc) due to Vercauteren's Remark 3.10.3) of two large polynomials. So I should include in routine mod 2^N (or mod 2^(Nc)...) instructions to avoid exponential coefficients' growing. But since the 2^N is not prime it's a problem - polynomials contain even coefficients and this leads to some even denominators - and for example multiplicative inverse 1/2 mod 2^N doesn't exist. Please tell me how to solve this problem?

How to adapt XGCD routine for correct mod 2^N calculation of resultant (norm)?

Thank you.

mod prime version of XGCD:

XGCD.mw

Hi everyone. I'm using Maple 2020. I encountered an error as "Error, (in SumTools:-DefiniteSum:-ClosedForm) summand is singular in the interval of summation". I saw this first time. Can you help me? I added source file.

Space_Fractional.mw

 

Hi,
I want to plot an equation, but I couldn't. Who can guide me?

Ra.mw

Suppose that a function f  has derivatives of all orders at a.  The the series

∑=(f(k)(a)/(k!))*(x−a)^k (limits are infinity and k=0, i donot how to insert that)

is called the Taylor series for f  about  a, where  f(n) is the n th order derivative of  f.

Suppose that the Taylor series for e2 x sin(5 x) about 0 is

a0+a1x+a2x2+⋯+a8x8+⋯

Enter the exact values of a0 and a8  in the boxes below.

      a0=      

     a8=      

Use Maple to find the solution of the initial value problem

y*(d^(2)*y/d*x^2)+(dy/dx)^2=0 0 with initial conditions y(0)=5and y'(0)=8.

Using Maple syntax, type in your answer in the box below, or copy (Ctrl-C) from your Maple worksheet and paste (Ctrl-V) in the answer box the solution. Do NOT enter the y(x)= part of the Maple output.

Why is pdsolve's 'generalsolution' option giving the particular solution u(x, y) = 0 instead of the general solution u(x, y) = A sin(x) sin(2 y) + B sin(2 x) sin(y) for the attached problem?

Problem.mw

Is there a way to convert a Fourier series (from the OrthogonalExpansions package) automatically into the sum of odd/even terms if the even/odd terms are 0 respectively?

Download FourierSeries.mw