Maple 2016 Questions and Posts

These are Posts and Questions associated with the product, Maple 2016

Here's a short tensor manipulation which goes totally bananas. Basically I have a metric, and I define a vector k, and right at the end I calcualte the covariant derivative of it. In the metric and elsewhere I have a constant epsilon, a function of time a(t), and 2 functions of all coords, Phi and psi. Then at the end it gives the covariant derivative of k, but now with epsilon and a as functions of all coordinates. 

Any idea what's going on??





Setup(metric = ds2);





the final output is Matrix(1, 1, [[(-(2*((diff(epsilon(X), r))*(psi(X))(X))+2*epsilon(X)*(diff((psi(X))(X), r)))*(epsilon(X)*(k1(X))(X)*((a(t))(X)^2)-1)*a(t)-(2*epsilon(X)*(psi(X))(X)-1)*((diff(epsilon(X), r))*(k1(X))(X)*((a(t))(X)^2)+epsilon(X)*((diff((k1(X))(X), r))*((a(t))(X)^2))+2*epsilon(X)*(k1(X))(X)*(a(t))(X)*(diff((a(t))(X), r)))*a(t)+epsilon*psi[r]*(epsilon*k1(X)*(a(t)^2)-1)*a(t)-epsilon^2*psi[theta]*(a(t)^3)*k2(X)-epsilon^2*psi[phi]*(a(t)^3)*k3(X)-(2*(diff(a(t), t))*psi(X)*epsilon+psi[t]*a(t)*epsilon-(diff(a(t), t)))*(epsilon*k0(X)*(a(t)^2)+1))*(1/a(t))]])

where epsilon is now a function!




The help text for dsolve,numeric,events describes many kinds of triggers and actions and other event specifications, but shows only a very limited selection of examples of these possibilities.

Please tell me of any sources containing extensive examples of the wide variety of available specifications and the problems they are useful in solving.

I was learning about Maple, but I do not understand this code below. This does not look correct.

When I type it on my worksheet, I get these errors



Sometimes I wonder if any one at Maplesoft actually looks at their own help web pages and try to use them, or is just the poor users who do that.



I noticed Maple gives different order of term (but still correct) when variable y or z is used, vs. other letters in this computation.

In this example below, integral results are given, using one letter of the alphabet at a time in the integrand. 

Only when the variable is y or z, Maple shows expression in different form from all the others (negative sign pulled out). The result are correct ofcourse, but was wondering why this hapens only for these two variables? I would have expected that variable name (letter used) should not make any difference to the final form.

chars:=[seq(parse(StringTools:-Char(i)),i=97..122)]: #generate a..z variables, may be there is better way?
Matrix([seq([chars[i0],simplify(int((c0^2 - c1^2)/(chars[i0] - x0+I*c3), x0))],i0=1..nops(chars))]);

Gives (this below. Notice the very last 2 entries below)

Is this something common in Maple? Does it depend on the computation being performed? Can a user do something to tell Maple not to change order/form of expression depending on what variable letter is used?



please help me for dsolve differential equations...after much time dont answer!!!


Digits := 15; SYS := [.16783*h1(theta)-0.96238e-3*(diff(h1(theta), theta, theta))+0.61603e-1*(diff(h2(theta), theta))+0.14870e-4*(diff(h2(theta), theta, theta, theta))-.23703*h3(theta)-0.84431e-3*(diff(h3(theta), theta, theta))+3.4919*10^(-7)*(diff(h1(theta), theta, theta, theta, theta)) = 0, 2.3940*h2(theta)-.35329*(diff(h2(theta), theta, theta))-0.68260e-1*(diff(h1(theta), theta))-0.16526e-4*(diff(h1(theta), theta, theta, theta))+3.0808*(diff(h3(theta), theta))-0.17833e-2*(diff(h3(theta), theta, theta, theta)) = 0, 9.4813*10^(-7)/((1.+1.5802*10^(-8)*h3(theta))*ln(10.+1.5802*10^(-7)*h3(theta))^2)-3.1867/((1.-0.26556e-1*h3(theta))*ln(10.-.26556*h3(theta))^2)-7.6530/((1.-0.31888e-1*h3(theta))*ln(10.-.31888*h3(theta))^2)-4.2551/((1.-0.35459e-1*h3(theta))*ln(10.-.35459*h3(theta))^2)-9.0315/((1.-0.37632e-1*h3(theta))*ln(10.-.37632*h3(theta))^2)-4.6587/((1.-0.38822e-1*h3(theta))*ln(10.-.38822*h3(theta))^2)-9.4520/((1.-0.39384e-1*h3(theta))*ln(10.-.39384*h3(theta))^2)+0.74999e-1/((1.+0.12500e-2*h3(theta))*ln(10.+0.12500e-1*h3(theta))^2)-.69143/((1.-0.28810e-2*h3(theta))*ln(10.-0.28810e-1*h3(theta))^2)-.12945/(1.-0.38836e-1*h3(theta))^4-0.38260e-1/(1.-0.11478e-1*h3(theta))^4-0.24826e-1/(1.-0.37240e-2*h3(theta))^4+0.74712e-3*(diff(h3(theta), theta, theta, theta, theta))+2.6337*10^(-8)/(1.+1.5802*10^(-8)*h3(theta))^4-0.27242e-1*(diff(h3(theta), theta, theta))-3.0707*(diff(h2(theta), theta))+0.17833e-2*(diff(h2(theta), theta, theta, theta))-.23618*h1(theta)-0.84126e-3*(diff(h1(theta), theta, theta))-0.89222e-1/(1.-0.26767e-1*h3(theta))^4-.21340/(1.-0.32010e-1*h3(theta))^4-.13154/(1.-0.19732e-1*h3(theta))^4-0.88519e-1/(1.-0.26556e-1*h3(theta))^4-4.7454/((1.-0.39545e-1*h3(theta))*ln(10.-.39545*h3(theta))^2)-0.36390e-1/(1.-0.10917e-1*h3(theta))^4-1.3100/((1.-0.10917e-1*h3(theta))*ln(10.-.10917*h3(theta))^2)-.89374/((1.-0.37240e-2*h3(theta))*ln(10.-0.37240e-1*h3(theta))^2)-1.3773/((1.-0.11478e-1*h3(theta))*ln(10.-.11478*h3(theta))^2)-.11820/(1.-0.35459e-1*h3(theta))^4-.26259/(1.-0.39388e-1*h3(theta))^4-4.7356/((1.-0.19732e-1*h3(theta))*ln(10.-.19732*h3(theta))^2)+28.586*h3(theta)-.21258/(1.-0.31888e-1*h3(theta))^4-9.4531/((1.-0.39388e-1*h3(theta))*ln(10.-.39388*h3(theta))^2)-4.6603/((1.-0.38836e-1*h3(theta))*ln(10.-.38836*h3(theta))^2)-9.0393/((1.-0.37664e-1*h3(theta))*ln(10.-.37664*h3(theta))^2)-4.2631/((1.-0.35526e-1*h3(theta))*ln(10.-.35526*h3(theta))^2)-7.6824/((1.-0.32010e-1*h3(theta))*ln(10.-.32010*h3(theta))^2)-.25088/(1.-0.37632e-1*h3(theta))^4-.12919/(1.-0.19378e-1*h3(theta))^4-.11842/(1.-0.35526e-1*h3(theta))^4-.12941/(1.-0.38822e-1*h3(theta))^4+0.20833e-2/(1.+0.12500e-2*h3(theta))^4-0.19206e-1/(1.-0.28810e-2*h3(theta))^4-4.6507/((1.-0.19378e-1*h3(theta))*ln(10.-.19378*h3(theta))^2)-.26255/(1.-0.39384e-1*h3(theta))^4-.13182/(1.-0.39545e-1*h3(theta))^4-.25109/(1.-0.37664e-1*h3(theta))^4-3.2120/((1.-0.26767e-1*h3(theta))*ln(10.-.26767*h3(theta))^2) = 0]

ode2 := diff(SYS[2], theta); SYS2 := {ode2, SYS[1], SYS[3]}; bcs2 := {h1(0) = 0, h1(1) = 0, h2(0) = 0, h2(1) = 0, h3(0) = 0, h3(1) = 0, ((D@@1)(h1))(0) = 0, ((D@@1)(h1))(1) = 0, ((D@@1)(h3))(0) = 0, ((D@@1)(h3))(1) = 0}; bcs22 := eval[recurse](convert(SYS[2], D), `union`({theta = 1}, bcs2)); res2 := dsolve(`union`(`union`(SYS2, bcs2), {bcs22}), 'maxmesh' = 2024, numeric, method = bvp[middefer], range = 0 .. 1, abserr = 0.1e-3, output = listprocedure)




Is there a way to automatically convert and paste my clipboard contents as 2D math?

Using a mouse's back/forward buttons or alt+ left arrow/ alt+right arrow did not work in Maple 2015 and they still don't work in Maple 2016. 

Also, if you select anywhere in the search bar, it automatically selects all the words.  I sometimes (most the time) want to change only one word.  

Im using x64 Windows 10 and it also did it when running Windows 8/8.1 .   Does anybody else have this problem?

Can it be fixed please if not?

hi..when i use rule [for] in maple code i encounter error'''''

Error, (in dsolve/numeric/process_input) input system must be an ODE system, found {f1(x), f2(x), f3(x), ApproximateInt(-4*cos(theta)^2, theta = 0. .. 1, method = simpson), ApproximateInt(4*cos(theta)^2, theta = 0. .. 1, method = simpson), ApproximateInt(8*cos(theta)^2, theta = 0. .. 1, method = simpson)}''''''''

please help me for remove it

i want for different beta for example beta=0, 40 and 80 this lines will computed three times  

 ''with(Student[Calculus1]); a1 := ApproximateInt(g1*g1, theta = a .. 1, method = simpson); a2 := ApproximateInt(2*(g1*g1)+3*g1*(diff(g1, theta, theta)), theta = a .. 1, method = simpson).......................''          

 by {for i from 1 by 1 to 3 do } and final gain ''ITRA_1_W[m] := eval(fy33*g3, fixedparameter)'' that have 3 amount.

when i use  rule {for i from 1 by 1 to 3 do ...  }   integral not computed and showed for example:::

a1=ApproximateInt(4*cos(theta)^2, theta = 0. .. 1, method = simpson)!!!!!!!!


I read a posting by Mr. Stephen Forrest on Thingiverse about using the exportplot command to export a 3D plot as a .stl file. I have a 3D printer and need the convert the plot from the Maple file into an .stl file.

Here is the posting I am referencing:


I was able to follow the commands successfully in the referenced paper for a hyperboloid or revolution that I plotted in Maple - it was very helpful. However, I could not figure out how to access the temporary .stl file that was created so that I could open it using my 3D printer's software.

My question is: Once I execute the "exportplot(stl, hyprev)" command that I inputted, how can I find the file in the Temporary Directory in order to open it from the software I use for my 3D printer?

I would appreciate any further details you can provide for accessing the temporary .stl file I created in order to be able to print the object.

Thank you!

If you try and scroll a page up or down, it just selects the content.  How hard would it be to fix this?

I use Maple on a Surface Pro 4 tablet and the math input panel does not work with Maple.  This would be a great feature to have.  Or does anybody know how to get it working?

Why can't I sub both values at the same time and/or why does k stay symbolic but not d?

Hello Dr/Prof?sir?madam

i have problem on running the ode bcs

there is cos n sine in ther ode

any idea to solve this ?

i have attched it

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/ .


Using Maple 2016, I created a new Document with one line:


Right-clicking on this expression I thought I could choose Assign to a Name, as is done in Clickable Calculus Series - Part 3: Multivariate Calculus, found here:

However, Assign to a Name does not appear in the context menu.

On the other hand, if I write


(with no parentheses) and right-click on this expression then Assign to a Name does appear in the context menu.

Any insights on this?


Can Maple determine the value of DthetaZero such that the solution to the ODE, for some specific value of t, simultaneously provides the two values theta(t)=Pi and diff(theta(t),t)=0?

sol:=dsolve({4*sin(theta(t))*cos(theta(t))-9.8100*sin(theta(t))-(diff(theta(t), t, t)) = 0, theta(0) = 2*Pi*(1/3), D(theta)(0) = DthetaZero}, numeric)

odeplot(sol,[t,theta(t),diff(theta(t),t)],t=0..5) for trial values of DthetaZero shows that the desired value is

1.0340*Pi <DthetaZero<1.0345*Pi.

First 26 27 28 29 30 31 32 Last Page 28 of 37