My procedure produces multiple different values for the same variable which is fine but I need a way of finding the lowest common multiple of all of the values. Is there a way to do this? Every time I try to do it, Maple gives me the last value it calculated.

Hello,

I used the folowing command, but the problem is that tehre is a blank space in the name of a directory.

!curl -m 3 -o /Users/John/Documents/P Admin/a.dat  --insecure https://www.google.com/?gws_rd=ssl;

I used macosx, and in a Terminal window the correct path must be written as :

/Users/John/Documents/P\ Admin/a.dat

But, impossible to set the right expression with Maple.

I used for instance:

!curl -m 3 -o /Users/John/Documents/P\\ Admin/a.dat  --insecure https://www.google.com/?gws_rd=ssl;

But it failed.

What is the right expression?

Best regards

Dears,

Realted to this post I am trying to solve the following inequalities in the variable "t"

Cosenos := proc (m, x) local C, x1, x2, y, x0, i, j, j0, k, K0, Kaux, K1, L, R, S, t1, t2, t, r1, r2, d; C := NULL; K0 := NULL; Kaux := NULL; K1 := NULL; L := NULL; R := NULL; S := NULL; d := nops(x); for k to d-1 do if evalf(x[k]) = 1 then j0 := m else j0 := 1+floor(x[k]*m) end if; K0 := K0, [(j0-1)/m, j0/m] end do; K0 := [K0]; S := K0[1][1] <= t, t <= K0[1][2]; for k from 2 to d-1 do S := S, K0[k][1] <= 1/2-(1/2)*cos(Pi*m^(k-1)*t), 1/2-(1/2)*cos(Pi*m^(k-1)*t) <= K0[k][2] end do; K1 := solve({S}, t) end proc

However, this run very very slow....For instance, for Cosenos(50, [.3225, .25877, .325, 1])

we interrupt the opertion after more than 10 minutes...I am doing something wrong?

Many thanks for your comments and suggestions. 

Hi, I use Maple 2017 on Windows 10 Pro. When I try to save a package that was included with the book entitled
"Finite Elements Using Maple - A Symbolic Programming Approach" by A. Portela and A. Charafi, Springer Verlag, Berlin, Heidelberg, 2002, with the line: "savelib('Plotter'):", maple return the error "Error, cannot open archive, C:\Program Files\Maple 2017\lib, for writing.". Someone can help me? Tank you.

I am new to Maple and recently installed the software. Surprisingly, it attempt to load by bringing the splash screen and later open the interface but goes off without any error or warning. Can anyone help me out as I really need to use it?

Dear Users!

Hope you would be fine. In the following maple code, I want to write the derivative of psi in term of psi like it did manually in red portion. For higher M and k it very hard to do it manully. It there any command to fix my problem for any value of k and M.

restart; k := 2; M := 4;

with(linalg); with(LinearAlgebra);

printlevel := 2;

for i while i <= 2^(k-1) do

for j from 0 while j <= M-1 do

psi[M*i+j-M+1] := simplify(2^((1/2)*k)*sqrt(GAMMA(j+1)*(j+alpha)*GAMMA(alpha)^2/(Pi*2^(1-2*alpha)*GAMMA(j+2*alpha)))*(sum((-1)^i1*GAMMA(j-i1+alpha)*(2*(2^k*x-2*i1+1))^(j-2*i1)/(GAMMA(alpha)*factorial(i1)*factorial(j-2*i1)), i1 = 0 .. floor((1/2)*j))));

`&psi;&psi;`[M*i+j-M+1] := simplify(diff(psi[M*i+j-M+1], x))

end do

end do; r := 2^(k-1)*M;

VV := Vector[column](r, proc (i) options operator, arrow; psi[i] end proc);

DV := Vector[column](r, proc (i) options operator, arrow; `&psi;&psi;`[i] end proc);

``&psi;&psi;`[2] := 8*sqrt((alpha+1)*(1/2))*sqrt(2)*sqrt(alpha*GAMMA(alpha)^2*4^alpha/GAMMA(2*alpha))/sqrt(Pi) = 8*sqrt((alpha+1)*(1/2))*psi[1];

`&psi;&psi;`[3] := 16*sqrt((2+alpha)*(alpha+1)/(1+2*alpha))/sqrt(2)*(2*sqrt(2)*sqrt((alpha+1)*GAMMA(alpha)^2*4^alpha/GAMMA(1+2*alpha))*alpha*(4*x+1)/sqrt(Pi)) = 16*sqrt((2+alpha)*(alpha+1)/(1+2*alpha))*psi[2]/sqrt(2)

I am waiting your response. Thanks

Hi everyone:

I want to obtain only the vertical asymptotes this function in the interval [a, b] and i used the asymptotes command but i didn't get true answer, can everybody help me? a=0.3 and b= 0.36

I attached my code. 

tnx...

Download Asymptotes.mw

Where is the mistake?
 

I'm working with tensors using the DifferentialGeometry package (not yet checked out the Physics package since I plan on working with Lie algebras), and in particular I'm computing the Riemann tensor of a metric.

Raising the first index of CurvatureTensor(g), the Riemann tensor has the symmetries of a symmetric product of two 2-forms, and this is how I would like to present it, identifying terms like dw1dw2dw1dw2 with terms like dw1dw2dw2dw1 with the appropriate sign. Ideally I would arrive at an expression that is written in terms of (dw1^dw2) (dw1^dw2) (etc.) basis elements, provided that my tensor has these symmetries.

Is there some readily available functionality to do this, or a simple, standard method in Maple?

Dear Users! 

Hope you would be fine with everything. I want to define a matrix F of M+1 by M+1 order having element of the following form:

I derived the F[r,s] but confuse who to generate matrix now.

restart; M := 5; printlevel := 3; for r from 2 while r <= M+1 do for s while s <= r-1 do if type(r+s, odd) then F[r, s] := 2^(k+1)*sqrt((2*r-1)*(2*s-1)) end if end do end do

I have got the following function:

R = ((n1-n2)/(n1+n2))^2

I would like to plot R as a function of (n1/n2). How can I do this in maple?

Hi everybody:

How can i plot this function in maple? for x=0.2..0.4 and I is complex number(I^2=-1).

tnx.

f:=(x)->-(2.50000000000000000000000000000*(-2.424176758*10^30*x^5+5.982448001*10^29*x^3+(5.336360731*10^33*I)*exp((4.0405213687327619451756171962*I)/x)+(3.144879905*10^98*I)*exp((48.6318841588794957405477169504*I)/x)+(8.489643736*10^35*I)*exp((4.04052136873276194517561719622*I)/x)))*x^6/((1.009157478*10^101*I)*exp((48.6318841588794957405477169504*I)/x)*x^6-2.760182453*10^32*x^11-(7.651751413*10^37*I)*exp((4.0405213687327619451756171962*I)/x)*x^4+1.887831524*10^32*x^9-(8.567811466*10^36*I)*exp((4.04052136873276194517561719622*I)/x)*x^4-(1.207009672*10^100*I)*exp((48.6318841588794957405477169504*I)/x)*x^4+(6.538572988*10^29*I)*x^6+(1.246166933*10^31*I)*x^10-1.837651455*10^31*x^7+(1.504191760*10^38*I)*exp((4.04052136873276194517561719622*I)/x)*x^6+1.574891043*10^34*exp((4.0405213687327619451756171962*I)/x)*x^5-(7.612867062*10^30*I)*x^8+4.556143298*10^99*exp((48.6318841588794957405477169504*I)/x)*x^5-(2.488678360*10^105*I)*exp((52.6724055276122576857233341466*I)/x)+6.791123646*10^36*exp((4.04052136873276194517561719622*I)/x)*x^5-3.454612059*10^36*exp((4.0405213687327619451756171962*I)/x)*x^3-4.710342238*10^98*exp((48.6318841588794957405477169504*I)/x)*x^3+(3.488285966*10^35*I)*exp((4.0405213687327619451756171962*I)/x)*x^6+5.512268188*10^106*exp((52.6724055276122576857233341466*I)/x)*x)

Im trying to create a basis as similar to this document. I have created the 72 equations under root of unity method (aswell have 44 equations on a seperate document for the sum shidoku product way) but cant seem to generate a basis for it such as they did in the document with these equations. If anyone can help as i may be entering the Basis() wrong. If you need the 72 equations i can post them as well. Any tips that you have would be appreciated!

http://reu.dimacs.rutgers.edu/~ericlee/Shidoku.pdf

Hello everyone,

I have a project that uses several machines to make a computation. I wrote a package (called here "package") with procedures, "proc".  I wrote several .mpl files structured as:

read "package.mpl":

W:= Read("file1.mpl"):    #I wrote a procedure to read and assign to a variable. This is not the problem.

Export("output.mpl", proc(W));

The situation is as follows: in all but one of my machines, this works fine, and proc(W) is exported correctly. However, in this ugly duckling (which is running the GUI version of Maple18), maple prints as an output, in blue:

Export("output.mpl",...)

and in place of ... it gives the answer expected. But it does not export automatically. I can, of course, manually correct this, but this is not optimal. Moreover, if the output is too big, there will be display issues.

Help? Any light is appreciated, thanks.

Marcelo

PS: I do not expect an answer for the PS; but I also do not know how to make my package "official", in the sense that I always read the .mpl file with all the procedures. Is there a way to call it like the LinearAlgebra and so on and use the package:-proc syntax? 

I have a problem with this sequence, as defined : f(x)=1-1/x+2*trunc(1/x)

The sequence starts with the number x=50100 and with 100000 steps the aim of the sequence is to see how many numbers are between 0 and 1( and 1 and 2, 2 and 3, and 3 and 4). My idea is:

 f:=x->1-2*(1/x-trunc(1/x))+1/x;
       x:=50100;
       k:=0 : a:=0 : b:=0 : c:=0 : d:=0 : while k<100000 do x:=f(x): k:=k+1 : if(x>0 and x<1)   then(a=a+1)           elif(x>1 and x<2) then(b=b+1) elif(x<2 and x>3) then(c=c+1) elif (x>3 and x<4) then (d=d+1) end             if : end do: printf("%d %d %d %d", a,b,c,d)                                                                                                                                             But i have problems :(.  Does anyone see the mistake?