Here is the simplest example of the problem. If I type

assume(j, integer, j>0); assume(k, integer, k > 0);

followed by

int(cos(j*t)*cos(k*t), t = -Pi..Pi);

Maple gives the answer 0, not, as I would have expected, Pi*KroneckerDelta[j,k].

Is there a way to force Maple to alllow for both the possibilities that j = k and j <> k, and thereby to give the general result using the Kronecker delta (or equivalent)?

Bought this book way back in 2000 to learn programming in Maple for procedures.

Problem was that i could not test myself with the exercises, because there were no answers included 

Has someone answers in Maple from these exercises ? 

That is the only way to learn from this book and check my knowledge about programming.

Jan 

The following spreadsheet contains an integral to be done

in the framework of the Maxwell–Jüttner distribution and therefore date back to 1911

integrale_juttner_primes.mw
 

Download integrale_juttner_primes.mw

Hi,

I am using Maple to do a parameter estimation on systems of ODEs using linear least-squares.

Though I have solved the problem, I would like to arrive into better symbolic expressions for the parameters, as they are obtained as the solution of systems of linear equations where the coefficient matrix and the RHS vector have sums as their elements.

The difficulty I have now is that the expressions I arrive at are, for instance:

and what I would like to do, now, would be to express this as

(2*h^2*sum(x[i]^2,i=1..4))*alpha+(2*h^2*sum(x[i]^2*y[i],i=1..4))*beta+2*h*sum(x[i]^2,i=1..4)-2*h*sum(x[i]*x[i+1],i=1..4)

x and y are already defined as vectors.

Is there any way I can do it?

Regards,

Rudnei

Hello,

Is it possible to solve a differential equation, where the solution is a 8x8 matrix u(t)?

The equation to solve for, "vgl", is this one:

The values for u(0) are known:

The functions QL(t) and RL(t) are also 8x8 matrix functions. They are actually Fourier series.

Can anyone point me in a direction how to tackle this in Maple?

For completness, I add the entire Maple worksheet (Toon.mw), but my last line obviously doesn't work ;-)

Also all the matrices for the Fourier series are added, just for completeness: RL0.txtRL1.txtRL_1.txtstrQL0.txtstrQL1.txtstrQL2.txtstrQL_1.txtstrQL_2.txtustart.txt

Any hints appreciated a lot!

Stef

Before Maple 2020, it was possible to start up Maple from the Linux commandline with multiple files, as  in

maple   -x   worksheet-1.mw   worksheet-2.mw   ...   worksheet-n.mw

That started up Maple's Java interface and loaded the n worksheets in separate tabs.

That no longer seems to work in Maple 2020 — now the command exits silently when invoked with multuple files.  Why?

Aside: Loading a single file works as expected, but not multiple files.

I'm still trying to get used to the commands. How do you find the height of a cylinder and plot to a graph?

Thanks in Advance.
 

Download heightofcylinder.mw

Hello

I wanted to create a group of interactive components that, when given an eqaution in a math container of an ellepses (in the form of ((x/A)^2)+((y/B)^2)=1) would pick A and B and give back two formulas of the type

x(t)=A*cos(t)

y(t)=b*sin(t)

So, is there a way to extract A and B from the initial formula?

Probably a dumb question, but I don't manage to read complex numbers from a text file, not even when the text file is created by Maple itself. I've been looking around, but didn't find an answer (which must be rather a simple one, I suppose).

Why is A not a matrix with complex numbers, as is B?

Hi, I'm on Maple 2020 on Mac, and running this:

" />

In fact, that isn't the Chromatic number of that graph. Is not possible to color that graph with 2 colors becuase C7 (Cycle with 7 vertex) is a subgraph which is imposible to color with 2 colors. Am I doing something wrong? Can you replicate this? Is this a bug or is me?

Thanks

Dear all

I want if possible to solve the inequation   f(x,y)>0 using maple

positive_function.mw

Thanks for any help

Maple formerly accepted

rule1 := forall(string(y), TD(y) = D(y)):

but now that statement prodices an error message:

Error, (in forall) expecting bound variable(s) to be a name or list of unique names

The entire value of this forall structure is that it should apply to objects of specific types, but it is not obvious how to express those types; the Help page for 'forall' is essentially useless in relation to this question. 

Hi, I am a beginner of maple and I met a problem recently.

I extracted some equations from a .txt file like below:

> read “D:/blablabla/filename”

                       EQ_1,EQ_2,...

Then I want to collect all these equations to make a list EQlist=[EQ_1,EQ_2,...] but the following code doesn't work:

> EQlist:=read “D:/blablabla/filename”

Error, invalid assignment

Could anyone tell me how to construct such EQlist?

Thank you in advance.  P.S. I can't just copy and paste because these commands are supposed to be in a loop.

Dear maple users,

Greetings.

How to plot this function equation "An" for x=0.0001..1,0.02 with 0..1 range.

Wating for replay.

restart;
A2 := 1.107444364; A4 := 1.124502164; ad := .5; ed := 0.1e-1; pd := 21; ld := .3;
f := unapply(3*x^2-2*x^3-1.238616691*x^2*(x-1)^2-.7382714588*x^2*(x-1)^3+.1921034396*x^2*(x-1)^4+.5253305667*x^2*(x-1)^5+.7364291997*x^2*(x-1)^6+1.032724351*x^2*(x-1)^7+.8058204155*x^2*(x-1)^8+.3290860035*x^2*(x-1)^9, x);
t := unapply(.339997432+1.547096375*x^2-2.488736512*x^3+8.154594212*x^4-15.63643668*x^5+15.85832377*x^6-8.734300202*x^7+1.959461605*x^8, x);
b1 := f(x);
b2 := diff(f(x), x);
b3 := diff(f(x), x, x);
b4 := t(x);
b5 := diff(t(x), x);
An := A4*(1+(4/3)*ad)*(b5^2+b4^2/ld^2)+4*ed*pd*(b2^2/x^2+b1^2/x^4-b1*b2/x^3)/ld;
As := seq(An, x = 0.1e-2 .. 1, 0.5e-1);
L := [seq([x, As], x = 0.1e-2 .. 1, 0.5e-1)];

with(plots);
plots:-display(plot(L, style = point), plot(As, x = 0.1e-2 .. 1), color = blue, linestyle = solid, labels = ["&eta;", "f'"], thickness = 1, labeldirections = [horizontal, vertical], labelfont = ['TIMES', 'BOLDOBLIQUE', 16], size = [450, 450], axes = box);
mp.mw
 

Download mp.mw

Is this a bug?

eq:=5*x+y-10=0:
eq2:=y=solve(eq,y):

f:=(x,y)->eq2
                              
f(3,1)
                   

I think the way I applied it above is a bit unorthodox to Maple but maybe it is a bug.  Normally I think we use unapply in this situation. 

f:=unapply(eq2,x,y)
                             
f(3,1)
                     

Okay better, but I think that first way should work no?