Maple 17 Questions and Posts

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


eq1 := a[n+2]+4*a[n+1]+10*a[n]-5*n-1;

applyrule(a[n::anything] = n*p+q+b[n], eq1);

eq2 := a(n+2)+4*a(n+1)+10*a(n)-5*n-1;

applyrule(a(n::anything) = b(n)+n*p+q, eq2);

eq3 := a(n+2)+4*a(n+1)+10*a(n) = 5*n+1;

applyrule(a(n::anything) = b(n)+n*p+q, eq3)


So my question is: Why does the first one not work? the second one works as expected, but the third one again does not?!

Good morning,

I've been having some trouble entering this expression in Maple 17.

My teacher suggested using this format, however, this format is from 2004 and certainly isn't recognized by the current version of Maple, as it keeps spitting back 0=0 no matter how many different ways I try to enter it.

Could anyone provide feedback on the correct syntax of how to enter this expression? Thank you very much for your time.

how may i encounter this problem,

Please help me on this ode solve (dsolve) wit



ps := 8933:

Eq1 := 8*(11-10*d)*(eta*(diff(f(eta), `$`(eta, 3)))+diff(f(eta), `$`(eta, 2)))/((1-phi)^2.5*(1-phi+phi*ps/pf))+4*f(eta)*(diff(f(eta), `$`(eta, 2)))+m*(1-4*(diff(f(eta), eta))*(diff(f(eta), eta))) = 0;

8*(11-10*d)*(eta*(diff(diff(diff(f(eta), eta), eta), eta))+diff(diff(f(eta), eta), eta))/((1-phi)^2.5*(1+7.958981045*phi))+4*f(eta)*(diff(diff(f(eta), eta), eta)) = 0


Eq2 := 2*knf*(eta*(diff(theta(eta), `$`(eta, 2)))+diff(theta(eta), eta))/(Pr*kf*(1-phi(eta)+phi*ps*cps/(pf*cpf)))+f(eta)*(diff(theta(eta), eta))-n*(diff(f(eta), eta))*theta(eta) = 0;

186.4367280*(eta*(diff(diff(theta(eta), eta), eta))+diff(theta(eta), eta))/(1-phi(eta)+.8253667629*phi)+f(eta)*(diff(theta(eta), eta)) = 0


Va := [.1, .1, .1]:

`Vφ` := [0.5e-1, 0.5e-1, 0.5e-1]:

etainf := 8:

bcs := f(a) = 0, (D(f))(a) = 0, theta(a) = 1, (D(f))(etainf) = 1/2, theta(etainf) = 0;

f(a) = 0, (D(f))(a) = 0, theta(a) = 1, (D(f))(8) = 1/2, theta(8) = 0


dsys := {Eq1, Eq2, bcs}:

for i to 3 do for h to 3 do phi := `Vφ`[i]; a := Va[h]; dsol[h][i] := dsolve(dsys, numeric, continuation = d); print(a); print(phi); print(dsol[h][i](0)) end do end do

Error, (in fproc) unable to store 'YP[3]' when datatype=float[8]





h bcs

Is is possible to tell the latex output what to do if it finds a symbol?

Suppose I have the following expression:

O3^2 + k_O3NO

The normal latex output just takes it as it is and also translates k_O3NO to k_O3NO with the underscore.

I however want it to read (latex notation in math mode)

\left[ \text{O}_3 \right]^2 + k_{\text{O}_3\text{NO}}

So is that possible without changing everything manually?


edit: Just as a remark, the output at the moment for this expression reads

{{\it O3}}^{2}+{\it k\_O3NO}

and is kinda awful if compiled in latex...

Is it also possible for maple to latex realizing where a line-break is necessary?


I'd like to plot two coupled functions x=x(t) and y=y(t) given by expression R0 and R1 (see below) for the interval t=0.001..0.999 and then after that extract tha data from the graphs for functions x(t) and y(t).

How can I do that in easiest way?

I will appreciate you.

Dear all,


I've just discover something quite surprising for me... It turns out like if Maple didn't take in account the linearity of the integral.


In the middle of a calculation, I need to compute the following integral

int(3628800 / (y * (1 / 2 + y)^11) - 3628800 / (y * (39 / 2 + y)^11), y = 39 / 2..infinity);

Maple 17 is not able to instantly answer me, and I'm not so patient... (I've waited more than 30s)

Nevertheless, if I ask Maple 17 to compute

3628800 * int(1 / (y * (1 / 2 + y)^11) - 1 / (y * (39 / 2 + y)^11), y = 39 / 2..infinity);

as well as

int(3628800/(y*(1/2+y)^11), y = 39/2..infinity)-int(3628800/(y*(39/2+y)^11), y = 39/2..infinity);

Maple 17 answer me instantly... which turn out to be quite surprising!

Is there someone, here, who has some ideas of what is happening in this calculation? Why Maple17 has this strange behaviour?

Many thanks for all your answers,


Why does



not simplify to




Hello everybody.

I'd like to share an observation about the integration of the Jacobi elliptic functions, in particularly the elliptic sine sn(x,k).

It's correct answer. But when I make the integration of the expression


according to the Handbook of Elliptic Integrals for Engineers and Scientists we have slightly another result

What is the difference?

how to calculate potential energy in terms of gauss curvature?



I have been trying to create a code to solve a structure by using FEM, but I keep getting an error that stops me from moving on.


Everytime I run my code, Maple gives me the error in the Title of my question:
"Error, number of indices exceeds rank"


Unfortunately I can't find any help in Maplesoft support for such error, and I need immediate help for that.

Can someone help me with that? Should I post my entire code?


I want to solve this equation with assumptions!!!

assume(d::real, d>0):
assume(a::real, -0.01 < a, a < 0):
sys:={-800*Pi*a*cos(6.557377048*Pi*(3.470797713+d))/(a+1)^3 = -.9396060697, 800*Pi*a*sin(6.557377048*Pi*(3.470797713+d))/(a+1)^3 = -.3238482794};
solve(sys, {a,d},useassumptions=true, AllSolutions=true);

one of the solutions has true "a" but "d" is wrong, I want one true solution!

I solve a mechanical exercise but i had a problem.

I know M (mass) and K (stifness) matrices (4x4).

I want to solve the (λ2M+K)v=0  eigenvalue problem, where λ are the eigenvalues and v eigenvectors.

How can i solve this problem.  I tried with the Eigenvectors() command but it didn't give the right solution.

The Eigenvalues are okay, but the eigenvectors not

K := Matrix([[4*10^7,-1.50*10^7,2*10^7,0],[-1.50*10^7,1.50*10^7,0,1.50*10^7],[2*10^7,0,8*10^7,2*10^7],[0,1.50*10^7,2*10^7,4*10^7]]);



Acoording with the book the right eigenvectors(shape mode) are:

[0.013 991,  0.034 233,  0.073 683,  0.090 573]
[0.035 637, 0, -0.032 213, 0]
[0 ,-0.034 233, 0, 0.090 573]
[-0.013 991, 0.034 233, -0.073 683, 0.090 573]

Thank you

Regarding my recent question I would like to share an interesting observation. Here the code of the program:

R0 := ln(y)+Re(Psi(1/2+(2*(p^2+(1/2)*sqrt(2*I+4*ksi_fs^2*p^2)*tanh(sqrt(2*I+4*ksi_fs^2*p^2)*x)/(tau+0.5e-2*a)))/y))+gamma+2*ln(2)
tau:= 10.000:ksi_fs:=10:p:=0.037:
R0p:= unapply(R0, [a,x]):
R0f:= proc(a,x)
local r:= fsolve(R0p(a,x), y= 0..1);
   `if`(r::float, r, Float(undefined))
end proc:
M:= Matrix(
   (i,j)-> R0f(i, 1 + (j-1)*(0.5-0)/(100-1)),
   datatype= float[8]

After approximately 2 hours of calculations I get a message window

But I repeat this calculations on another computer with the same Windows 7 64 bit and Maple 17 I don't get such error and I obtain desired data.

So can Maple be sensitive to the hardware? 

I'm trying to numerically solve the differential equation: y' = -2xy + 1. Naturally, I come across the non-elementary integral of e^(x^2). By hand, I used a 2nd degree MacLaurin polynomial to get y = xe^(-x^2) + x^3/3e^(-x^2)+x^3/6e(-x^2). 

How do I use Maple to numerically solve this, with step sizes of h=0.1 and h=0.05 and plot them?

Hi all, I'm using text box and push button to create my user interface for my system. I  wish to use a if..then function to check some condition, if the condition is true then do statement and print out in the text box. When is key coding in the push button and click the push button, then if ...then function not working and the result that not suitable for the condition also have been printed out in text box.So,can someone help my to solve the problem?

Here in the coding i did

use DocumentTools in 

Do (C = %txtC);
Do (p = %txtpk1);
Do (q = %txtpk2);
Do (d = %txtpk3);
Do (A1 = %txtA1);
Do (A2 = %txtA2);
Do (w = `mod`((C*d), A2));
Do (mp = `mod`(w^((p+1)*(1/4)), p));
Do (mq = `mod`(w^((q+1)*(1/4)), q));
Do (hl = `mod`((mp-mq)/p, q));
Do (h2 = `mod`((-mp-mq)/p, q));
Do (m1 = mp+(hl*p));
Do (m2 = mp+(h2*p));
Do (m3 = (p*q)-m1);
Do (m4 = (p*q)-m2);
if m1 <= 2^(2*k-1) then t1 := (C+(-A1*(m1^2)))/A2 end if;
Do (t1);
if m2 <= 2^(2*k-1) then t2 := (C+(-A1*(m2^2)))/A2 end if;
Do (t2);
if m3 <= 2^(2*k-1) then t3 := (C+(-A1*(m3^2)))/A2 end if;
Do (t3);
if m4 <= 2^(2*k-1) then t4 := (C+(-A1*(m4^2)))/A2 end if;
Do (t4);

Do (%txtw = w);
Do (%txtm1 = m1);
Do (%txtm2 = m2);
Do (%txtm3 = m3);
Do (%txtm4 = m4);
Do (%txtt1 = t1);
Do (%txtt2 = t2);
Do (%txtt3 = t3);
Do (%txtt4 = t4);
end use; 

Thabk you.

First 15 16 17 18 19 20 21 Last Page 17 of 61