Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi MaplePrimes,

This YouTube video has a nice puzzle. 

It is titled "Can you solve the locker riddle". 
My first blush was to consider modular arithmatic
https://www.youtube.com/watch?v=c18GjbnZXMw

Here is a maple page -
divisors_excercise.mw

divisors_excercise.pdf

Have a very fine rest of the day.

Regards,
Matt

 

I uploaded a test document to my private space in Maple Cloud and I now want to delete it. How does one do that?  Also I just put it in the public space and want to delete it from there.


 

diff(Q(t), t) = k*A*(T__1s(t)-T__2s(t))/d

diff(Q(t), t) = h__1(t)*A*(T__1(t)-T__1s(t))

diff(Q(t), t) = h__2(t)*A*(T__2s(t)-T__2(t))

Q(t) = m__1*c__p*(T__1i-T__1(t))

Q(t) = m__2*c__p*(T__2(t)-T__2i)

h__1(t) = k(T__1(t), T__1s(t))*(.825+.387*(g*h^3*c__p*beta(T__1(t), T__1s(t))*rho(T__1(t), T__1s(t))^2*(T__1(t)-T__1s(t))/(k(T__1(t), T__1s(t))*mu(T__1(t), T__1s(t))))^(1/6)/(1+(.492*k(T__1(t), T__1s(t))/(c__p*mu(T__1(t), T__1s(t))))^(9/16))^(8/27))^2/h

h__2(t) = k(T__2(t), T__2s(t))*(.825+.387*(g*h^3*c__p*beta(T__2(t), T__2s(t))*rho(T__2(t), T__2s(t))^2*(T__2s(t)-T__2(t))/(k(T__2(t), T__2s(t))*mu(T__2(t), T__2s(t))))^(1/6)/(1+(.492*k(T__2(t), T__2s(t))/(c__p*mu(T__2(t), T__2s(t))))^(9/16))^(8/27))^2/h

 

 

rho(T__1(t), T__1s(t)) = 999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4

beta(T__1(t), T__1s(t)) = -(4.216485*10^(-2)-7.097451*10^(-3)*(T__1(t)+T__1s(t))+2.63217825*10^(-5)*(T__1(t)+T__1s(t))^2-4.9518879*10^(-8)*(T__1(t)+T__1s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4)
mu(T__1(t), T__1s(t)) = 2.414*10^(247.8/(.5*(T__1(t)+T__1s(t))+133)-5)

k(T__1(t), T__1s(t)) = -9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949

 

 

rho(T__2(t), T__2s(t)) = 999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4

beta(T__2(t), T__2s(t)) = -(4.216485*10^(-2)-7.097451*10^(-3)*(T__2(t)+T__2s(t))+2.63217825*10^(-5)*(T__2(t)+T__2s(t))^2-4.9518879*10^(-8)*(T__2(t)+T__2s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4)
mu(T__2(t), T__2s(t)) = 2.414*10^(247.8/(.5*(T__2(t)+T__2s(t))+133)-5)

k(T__2(t), T__2s(t)) = -9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949

 

"`h__1`(t)=(-9.481411*10^(-6) (`T__1`(t)+`T__1s`(t))^(2)+2.1356735*10^(-3) (`T__1`(t)+`T__1s`(t))+0.5599920949)/(h) (0.825+(0.387 ((g h^(3) `c__p` (-(4.216485*10^(-2)-7.097451*10^(-3) (`T__1`(t)+`T__1s`(t))+2.63217825*10^(-5) (`T__1`(t)+`T__1s`(t))^(2)-4.9518879*10^(-8) (`T__1`(t)+`T__1s`(t))^(3))/(999.9399+2.1082425*10^(-2) (`T__1`(t)+`T__1s`(t))-1.77436275*10^(-3) (`T__1`(t)+`T__1s`(t))^(2)+0.438696375*10^(-5) (`T__1`(t)+`T__1s`(t))^(3)  -0.6189861563*10^(-8) (`T__1`(t)+`T__1s`(t))^(4))) (999.9399+2.1082425*10^(-2) (`T__1`(t)+`T__1s`(t))-1.77436275*10^(-3) (`T__1`(t)+`T__1s`(t))^(2)+0.438696375*10^(-5) (`T__1`(t)+`T__1s`(t))^(3)  -0.6189861563*10^(-8) (`T__1`(t)+`T__1s`(t))^(4))^(2) (`T__1`(t)-`T__1s`(t)))/((-9.481411*10^(-6) (`T__1`(t)+`T__1s`(t))^(2)+2.1356735*10^(-3) (`T__1`(t)+`T__1s`(t))+0.5599920949) 2.414*10^((247.8)/(0.5 (`T__1`(t)+`T__1s`(t))+133)-5)))^((1)/(6)))/((1+((0.492 (-9.481411*10^(-6) (`T__1`(t)+`T__1s`(t))^(2)+2.1356735*10^(-3) (`T__1`(t)+`T__1s`(t))+0.5599920949))/(`c__p` 2.414*10^((247.8)/(0.5 (`T__1`(t)+`T__1s`(t))+133)-5)))^((9)/(16)))^((8)/(27))))^(2)"

"`h__2`(t)=(-9.481411*10^(-6) (`T__2`(t)+`T__2s`(t))^(2)+2.1356735*10^(-3) (`T__2`(t)+`T__2s`(t))+0.5599920949)/(h) (0.825+(0.387 ((g h^(3) `c__p` (-(4.216485*10^(-2)-7.097451*10^(-3) (`T__2`(t)+`T__2s`(t))+2.63217825*10^(-5) (`T__2`(t)+`T__2s`(t))^(2)-4.9518879*10^(-8) (`T__2`(t)+`T__2s`(t))^(3))/(999.9399+2.1082425*10^(-2) (`T__2`(t)+`T__2s`(t))-1.77436275*10^(-3) (`T__2`(t)+`T__2s`(t))^(2)+0.438696375*10^(-5) (`T__2`(t)+`T__2s`(t))^(3)  -0.6189861563*10^(-8) (`T__2`(t)+`T__2s`(t))^(4))) (999.9399+2.1082425*10^(-2) (`T__2`(t)+`T__2s`(t))-1.77436275*10^(-3) (`T__2`(t)+`T__2s`(t))^(2)+0.438696375*10^(-5) (`T__2`(t)+`T__2s`(t))^(3)  -0.6189861563*10^(-8) (`T__2`(t)+`T__2s`(t))^(4))^(2) (`T__2s`(t)-`T__2`(t)))/((-9.481411*10^(-6) (`T__2`(t)+`T__2s`(t))^(2)+2.1356735*10^(-3) (`T__2`(t)+`T__2s`(t))+0.5599920949) 2.414*10^((247.8)/(0.5 (`T__2`(t)+`T__2s`(t))+133)-5)))^((1)/(6)))/((1+((0.492 (-9.481411*10^(-6) (`T__2`(t)+`T__2s`(t))^(2)+2.1356735*10^(-3) (`T__2`(t)+`T__2s`(t))+0.5599920949))/(`c__p` 2.414*10^((247.8)/(0.5 (`T__2`(t)+`T__2s`(t))+133)-5)))^((9)/(16)))^((8)/(27))))^(2)"

diff(Q(t), t) = k*A*(T__1s(t)-T__2s(t))/d, diff(Q(t), t) = A*(T__1(t)-T__1s(t))*(-9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949)*(.825+.387*(((-g*h^3*c__p*(4.216485*10^(-2)-7.097451*10^(-3)*(T__1(t)+T__1s(t))+2.63217825*10^(-5)*(T__1(t)+T__1s(t))^2-4.9518879*10^(-8)*(T__1(t)+T__1s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4))*((999.9399+2.1082425*10^(-2)*(T__1(t)+T__1s(t))-1.77436275*10^(-3)*(T__1(t)+T__1s(t))^2+.438696375*10^(-5)*(T__1(t)+T__1s(t))^3-.6189861563*10^(-8)*(T__1(t)+T__1s(t))^4)^2))*(T__1(t)-T__1s(t))/(2.414*10^(247.8/(.5*(T__1(t)+T__1s(t))+133)-5)*(-9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949)))^(1/6)/(1+((.492*(-9.481411*10^(-6)*(T__1(t)+T__1s(t))^2+2.1356735*10^(-3)*(T__1(t)+T__1s(t))+.5599920949))/(2.414*c__p*10^(247.8/(.5*(T__1(t)+T__1s(t))+133)-5)))^(9/16))^(8/27))^2/h, diff(Q(t), t) = A*(T__2s(t)-T__2(t))*(-9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949)*(.825+.387*(((-g*h^3*c__p*(4.216485*10^(-2)-7.097451*10^(-3)*(T__2(t)+T__2s(t))+2.63217825*10^(-5)*(T__2(t)+T__2s(t))^2-4.9518879*10^(-8)*(T__2(t)+T__2s(t))^3)/(999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4))*((999.9399+2.1082425*10^(-2)*(T__2(t)+T__2s(t))-1.77436275*10^(-3)*(T__2(t)+T__2s(t))^2+.438696375*10^(-5)*(T__2(t)+T__2s(t))^3-.6189861563*10^(-8)*(T__2(t)+T__2s(t))^4)^2))*(T__2s(t)-T__2(t))/(2.414*10^(247.8/(.5*(T__2(t)+T__2s(t))+133)-5)*(-9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949)))^(1/6)/(1+((.492*(-9.481411*10^(-6)*(T__2(t)+T__2s(t))^2+2.1356735*10^(-3)*(T__2(t)+T__2s(t))+.5599920949))/(2.414*c__p*10^(247.8/(.5*(T__2(t)+T__2s(t))+133)-5)))^(9/16))^(8/27))^2/h, Q(t) = m__1*c__p*(T__1i-T__1(t)), Q(t) = m__2*c__p*(T__2(t)-T__2i)

diff(Q(t), t) = k*A*(T__1s(t)-T__2s(t))/d, diff(Q(t), t) = A*(T__1(t)-T__1s(t))*(-0.9481411000e-5*(T__1(t)+T__1s(t))^2+0.2135673500e-2*T__1(t)+0.2135673500e-2*T__1s(t)+.5599920949)*(.825+.387*(-.4142502071*g*h^3*c__p*(0.4216485000e-1-0.7097451000e-2*T__1(t)-0.7097451000e-2*T__1s(t)+0.2632178250e-4*(T__1(t)+T__1s(t))^2-0.4951887900e-7*(T__1(t)+T__1s(t))^3)*(999.9399+0.2108242500e-1*T__1(t)+0.2108242500e-1*T__1s(t)-0.1774362750e-2*(T__1(t)+T__1s(t))^2+0.4386963750e-5*(T__1(t)+T__1s(t))^3-0.6189861563e-8*(T__1(t)+T__1s(t))^4)*(T__1(t)-T__1s(t))/(10^(247.8/(.5*T__1(t)+.5*T__1s(t)+133)-5)*(-0.9481411000e-5*(T__1(t)+T__1s(t))^2+0.2135673500e-2*T__1(t)+0.2135673500e-2*T__1s(t)+.5599920949)))^(1/6)/(1+.4087338992*((-0.9481411000e-5*(T__1(t)+T__1s(t))^2+0.2135673500e-2*T__1(t)+0.2135673500e-2*T__1s(t)+.5599920949)/(c__p*10^(247.8/(.5*T__1(t)+.5*T__1s(t)+133)-5)))^(9/16))^(8/27))^2/h, diff(Q(t), t) = A*(T__2s(t)-T__2(t))*(-0.9481411000e-5*(T__2(t)+T__2s(t))^2+0.2135673500e-2*T__2(t)+0.2135673500e-2*T__2s(t)+.5599920949)*(.825+.387*(-.4142502071*g*h^3*c__p*(0.4216485000e-1-0.7097451000e-2*T__2(t)-0.7097451000e-2*T__2s(t)+0.2632178250e-4*(T__2(t)+T__2s(t))^2-0.4951887900e-7*(T__2(t)+T__2s(t))^3)*(999.9399+0.2108242500e-1*T__2(t)+0.2108242500e-1*T__2s(t)-0.1774362750e-2*(T__2(t)+T__2s(t))^2+0.4386963750e-5*(T__2(t)+T__2s(t))^3-0.6189861563e-8*(T__2(t)+T__2s(t))^4)*(T__2s(t)-T__2(t))/(10^(247.8/(.5*T__2(t)+.5*T__2s(t)+133)-5)*(-0.9481411000e-5*(T__2(t)+T__2s(t))^2+0.2135673500e-2*T__2(t)+0.2135673500e-2*T__2s(t)+.5599920949)))^(1/6)/(1+.4087338992*((-0.9481411000e-5*(T__2(t)+T__2s(t))^2+0.2135673500e-2*T__2(t)+0.2135673500e-2*T__2s(t)+.5599920949)/(c__p*10^(247.8/(.5*T__2(t)+.5*T__2s(t)+133)-5)))^(9/16))^(8/27))^2/h, Q(t) = m__1*c__p*(T__1i-T__1(t)), Q(t) = m__2*c__p*(T__2(t)-T__2i)

(1)

"(->)"

``

``

``

i have a system with 5 dif equations and five unknows. i have told to maple to solve it numerically with interactively solve comand (right cilck button). the window open like it normally does and i put values to my parameters, with an initial condition for the system (Q(0)=0). then i press numerically solve and that's all, the program just keep evaluating with no answer. i wait for 15 min, which i think is too much time, and got any answer yet.

hope you can help with this

thanks.. 
 

Download propuesta_transfer.mw

Hello

I have a procedure that builds an ideal from a specific set of polynomials and then calls the Groebner basis package to eliminate some of the variables.  Even though the procedure is running on a machine with 2 processors, 24 cores and 72 GB of ram, only one core has been used (and is always on a 100% usage).  Would the Grid Computing Toolbox be of some hope in this case?  If so,  how to insert the Grid commands so that Maple sends the calculations to the other cores (I find the document rather confusing)?   If I am talking non sense,  please let me know.

Many thanks

Ed

 

 

 

Hi guys,

I would like to compute the complex roots of the following equations

u*(BesselJ(0,u)^2 + BesselJ(1,u)^2) = 2 BesselJ(0,u)*BesselJ(1,u)

The function fsolve in Maple gives only 0. I was wondering whether other complex solutions could be obtained as well.

Your help is highly appreciated.

Thanks

Federiko

with(DifferentialGeometry):with(JetCalculus):
DGsetup([x],[u],E,5);
vars≔x,u,u[1],u[1,1],u[1,1,1];
PDEtools[declare](Q(vars));
TotalDiff(Q(vars),x);
TotalDiff(u[1,1],x);

 

Hi everyone,

Recently I came across the total differentiation command in the PDEtools. For its

documentation, I used the following link

http://www.maplesoft.com/support/help/Maple/view.aspx?path=DifferentialGeometry/JetCalculus/TotalDiff

Unfortunately, when I try to replicate this it did not work as expected. I am getting the total derivative of the expression to be zero. I do not understand where I am going wrong.

You can find my code above. I am also attaching the screen shot of my maple file.

I would really appreciate if someone could help me out. Thanks for your help.


 

is there a possibility to display the partial derivative of f(x,y) wrt. y as

and second order / cross-partials accordingly?

 

thx

jo

 

 

limit((x^2-1)*sin(1/(x-1)), x = infinity, complex);
infinity-infinity*I
MultiSeries:-limit((x^2-1)*sin(1/(x-1)), x = infinity, complex);
infinity

whereas the same outputs are expected. The help http://www.maplesoft.com/support/help/Maple/view.aspx?path=infinity&term=infinity does not shed light on the problem. Here are few pearls:

  • infinity is used to denote a mathematical infinity, and hence it is usually used as a symbol by itself or as -infinity.
  • The quantities infinity, -infinity, infinity*I, -infinity*I, infinity + y*I, -infinity + y*I, x + infinity*I and x - infinity*I, where x and y are finite, are all considered to be distinct in Maple. However, all 2-component complex numerics in which both components are infinity are considered to be the same (representing the single point at the "north pole" of the Riemann sphere).
  • The type cx_infinity can be used to recognize this "north pole" infinity.

I have created several animated 3-d figures using Maple.  Exporting an animation as a gif file I find that with, for example, Quicktime, I can run the animation. However I cannot rotate the 3-d figures using Quicktime.  Is there some other way to make such animations available to people not having Maple, so they can rotate the figure as well as run the animation? I would appreciate any ideas about this.

The following program hangs on the last command and a hard restart is required. The computation of a 2 x 2 matrix times a 2-vector is not that hard. Any ideas as to what is happening?

Another question: if v is a vector that depends on x and y say why does
>solve(v=0,{x,y})
not work?

It should only take a few lines of code to change v=0 to the system {components of v = 0}

I have been using Maple for many years and have created many Maple procedures for my own use in different projects. My usual practice is to place the procedures in a .txt file using WordPad or NotePad on my Windows machines. Then to load, for example, a file named C:\Users\Ed\Desktop\Top Quandles\Procedures.txt I use the command 

read("C:\\Users\\Ed\\Desktop\\Top Quandles\\Procedures.txt");
or
read("C:/Users/Ed/Desktop/Top Quandles/Procedures.txt");


Or I can leave the path off if the file is in the same folder as my worksheet.

This has always worked well for me. But have two questions:

1. Is there is a better way to save and load a bunch of procedures? Also I might mention that I frequently have to go into the file and change procedures or add new procedures, data, etc...

2. Recently I had occasion to send a worksheet and a txt file to a Mac user. We were not able to get the worksheet I sent him to read the .txt file on his Mac. Any suggestions on how to do that on a Mac? We tried putting them both in the same folder and just using a simple read("Procedures.txt"); command but the Mac didn't like that.

Any suggestions about either of these two issues would be appreciated.

Thanks,

Edwin

Does anyone know if there is a simple way to write a falling factorial (pochhammer symbol) in Maple. It seems pochhammer is the rising factorial with no option to change to falling.

Thanks in advance.

I assigned

before an algebraic calculation so I would like to get  or have the program print the 70 digits of the answer and not just 10 digits. Because when I press ENTER, I get only 10 digits.

 

fieldplot is a wonderful tool for plotting vector fields. The option 'fieldstrength' is very useful to scale the arrows so that one can better visualize the field. I often use fieldstrength=log.

However, if one includes the RealDomain library, the fieldstrength=log option fails. I don't see why it should. Can someone enlighten me before I report this as a bug.

Hi all, 
I was only wondering if there is a way how to trunk a solution using a maple command.

I want to use it to give a truncation error which is equal to :(14/45)(D^5)(y)(0) h^5

The result I am getting is as followed: (it is correct, I just want to cut off the bit with the power of 6)
........
......
>error:=expand(Yx[i+1]-Yx[i-3]-(4*h/3)*(2*f[i-2]-f[i-1]+2*f[i]));
                          
          error :=  (14/45)(D^5)(y)(0) h^5  + 7/10 (D^6)(y)(0) h^6
                  
thank you in advance

First 1024 1025 1026 1027 1028 1029 1030 Last Page 1026 of 2216