499 Reputation

16 Badges

11 years, 70 days

MaplePrimes Activity

These are replies submitted by Alger


Thank you for your answer.

My objective id to define a Fourier serie of type "a0+Sum(an(n)*cos(n*Pi/ln(x2/x1)*ln(x/x1))+ bn(n)*sin(n*Pi/ln(x2/x1)*ln(x/x1)),n=1..N)" in range [x1..x2]

However, int(C(n,x)*S(m,x)/x,x=x1..x2) is not null.


C := (n,x) -> cos(n*Pi*ln(x/x1)/ln(x2/x1));

S := (n,x) -> sin(n*Pi*ln(x/x1)/ln(x2/x1));

Your above answer "Fact" is not simple for me to understand. 

I don't know if it's possible to get this type of serie.


@Rouben Rostamian  

Thanks for your answer.

My problem is int(C(n,x)*S(m,x)/x,x=x1..x2) is not null.

@Rouben Rostamian  

Thank you for for your above answer.

The basis cos(n*Pi*ln(x)), sin(n*Pi*ln(x)) is orthogonal in the period 1..exp(2)

My goal if possible to get a basis (orthogonal) as

cos(n*Pi/ln(x2/x1)*ln(x/x1)) and sin(n*Pi/ln(x2/x1)*ln(x/x1)) in the range [x1 x2]


@Rouben Rostamian  



a0+Sum(an*cos(n*Pi*ln(x))+bn*sin(n*Pi*ln(x)),n=1..N) is ok

In the period 1 to exp(2) it is biorthogonal

Thank you very mach for this post.

I'm looking for more improvement for pdsolve as:

pde[25] with f(x,y) as source and not 0.

pde[23] with f(r,theta) as source and not 0.

and other improvements for Helmoltz equations



Thank you


Yes, Tabulate do a good work

Thank you


Many thanks for your great answer. I will read it carrefully. I am very greatfull. 

Have the best

Thank you 

@Preben Alsholm 

The post is just a simple example.

I want that i is shown (print(i)) at the same time the for command is executing.

This is to know at what iterator the for command is during the calculation.

Assume for example we have a procedure that take 15 minutes and the i must vary from 1 to 20

for i from 1 to 20 do x[i]:=f_proc(i); print(i); end do:

I want to know with print(i) the stage of calculation at same time of calculation

I do'nt know if it's clear


Thank you.

Yes it is a good idea that Maple define my local variables.

Thank you for the answer.

I solved the problem when changing e1 as

e1:=seq(map(evalf,s1), j=1..Qr):


Thank you for the response.

I will think on the two propositions.

I have many local variables.

@Carl Love 

I understand the diffivulty.

Thank you


Thank you for the answer

I want to use MatrixFunction for BesselY(A,B) ?

Is it possible or not ?

@Adri van der Meer 

Thank you for the answer


MatrixFunction(A, sin(x), x); is not the same as sin~(A) ?

I understand that we can not use MatrixFunction for BesselY(A,B) ?


1 2 3 4 5 6 7 Last Page 2 of 17