MaplePrimes Questions



I tried to solve some linear differential equations,

but the result is incomplete.

The complete solutions have been derived by hand.

How can I get the right result using maple?

The example is from determining symmetres of differential equations.

Thank you!


How to solve this problem? I want to display plot of differential equation system

this download link my problem

Thank you :)

  1. I have a set of equations described on different lines and I want them to be executed together and their answers to be displayed in seperate lines after that. How do I do that?
  2. This is particulalry important because, at times I will have to document multiple equations sequentially and later on have their answers reffere to. This way its more structured and I can group equations falling into certain category.

Please find the attached Minimum working example. Any advise is appreciated. File link: []






i look for an Maple animation about similar triangles? Thank you

HI MaplePrimes,

The input -


returns a large expression.

My had calculations reduce this to

f(k) = [(41-19*sqrt(41))/820]*[((1-sqrt(41))/2)^k+((1+sqrt(41))/2)^k)].

There may be an error.

We let f(1)=1 and f(2)=2.

The sequence, starting with 1 should read -


What is the correct expression for f(k)?




Consider a standard initial/boundary value problem for the heat equation on the interval x ∈ [0,1]:

pde := diff(u(x,t),t) = diff(u(x,t),x,x);
ic := u(x,0) = f(x);
bc := u(0,t)=0,  u(1,t)=0;

pdsolve({pde, ic, bc});
produces the expected Fourier series solution.

However, if we change the interval to x ∈ [-1,1], as in:
bc := u(-1,t)=0,  u(1,t)=0;
pdsolve({pde, ic, bc});

then Maple fails to return a solution.  Why?

Dear all

I need a help how can i solve for example the following PDEs with Initial condition and boundary condition given at x=-1, and x=1.





Many thanks




In using Isolate in RootFinding to compute roots of a real polynomial, the output contains, say, z= some number.  How to get rid of the "z =" so that I can declare that "some number" to be some variable?

I am writing a matrix to file where each matrix element is placed on a new line. An example of such code:


nRow,nCol :=LinearAlgebra[Dimension](HH);
for i from 1 to nRow do
   for j from 1 to nCol do

When this is printed to file, it will print as the power sign ^ regardless of the input it is given. This particular output file will be read in by another language, and due to this the power symbol "^" is not desired, but instead " ** ". By converting the matrix into a string format and applying string operations, this can be done; but is there a simpler way to do it?



I am solving a simple system of algebraic equations in a "for" loop. I would like to write all the values of solutions to a matrix, but I only get the symbols. Here is my code:

for i from 1 to m do
 rozw||i:=fsolve({r1||i, r2||i}, {H||i,V_A||i},H||i=0..1);Force(i,1):=H||i;Force(i,2):=V_A||i;subs(rozw||i,Force(i,..));
end do:

Thanks in advance,


I need to create randomly several {1,2,3}-sequences of length 15. I execute:

> with(RandomTools):

> w:=Generate(list(posint(range = 3), 15)):

The problem is that when I close the program and then open and run it again, it always produces the same sequences. 

What should I do in order to make it work really randomly? That is, to start at a "random" sequence as well?


Can any one help me correct my doc wherein I want to record 5 minutes with microphone and play it for the same duration through my speaker with a click of button component. I find problem in file path specification.This document in my PC is  in E:\A_PhD2017\MAPLEPhDFiles\

Components added are only three: microphone0, Speaker0 and Button0





rate := DocumentTools:-GetProperty("Microphone0", samplerate); 1; R := DocumentTools:-GetProperty("Microphone0", value)

Array(%id = 18446745571883247742)


Edit Click Action for Button contains the following commands.

R := DocumentTools:-GetProperty("Microphone0", value); 1; srate := DocumentTools:-GetProperty("Microphone0", samplerate); 1; DocumentTools:-SetProperty("Speaker0", samplerate, srate); 1; DocumentTools:-SetProperty("Speaker0", value, R)




No errors displayed for both the above doc tools use``

Also nonstop continuous evaluation takes place! Where is the loop formation if at all there is any?

Thanks for answering

Ramakrishnan V





I'm trying to write and simplify expressions involving partial derivatives of an arbitrary function, say f(x,y).

Specifically, I would like diff(f(t,y),t) to evaluate to D[1](f)(t,y), instead of 

And eval(diff(f(x,y),x),x=t) gives the same result. 

Interestingly, diff(f(2*t,y),t) does evaluate to 2*D[1](f)(2*t,y), as expected.

I could get some results by using a custom differentiation function

`diff/f` := proc(x,y,v) 
    if v = x then return D[1](f)(v,y):
    elif v = y then return D[2](f)(x,v):
    else return 0: # is there a fallback I can use without risking infinite recursion?
end proc:

and it works for the simple case, but now diff(f(2*t,y),t) does not work.

All this makes me feel that there must be a way to get what I want by default, but I can't figure it out.

Is there? Thank you in advance!

Hi Mapleprimes people and robots,

My question is regarding a recursive sequence.  It can be defined non-recursively as - 

a(r) :=  0.8*3^r + 0.2*(-2)^r.

The first few terms are - 

1,2,8,20,68,188, and so on.

Here is my Maple Worksheet.      recursive_sequence_A133467.pdf

I want some Maple code that will produce 30 terms of this sequence.  It is defined as


for n>2 we let s[n] = s[n-1] + 6*s[n-2].

Let me know if my question does not make sense.



I have a second problem which is related to converting horizontal axis to degrees ( or in terms of Pi :  0 - Pi/2 ).

Thank you for your helps


s1 := (1/3*(4*(1.66+10^(-1))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-1)))-sin(2*a);
s2 := (1/3*(4*(1.66+10^(-2))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-2)))-sin(2*a);
s3 := (1/3*(4*(1.66+10^(-3))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-3)))-sin(2*a);
k1 := series(s1, x);
       4.769067180 x     + (-0.5773502688 - 1. sin(2. a))

                           2    / 4\
          + 0.01164914507 x  + O\x /
k2 := series(s2, x);
       4.293792441 x     + (-0.5773502692 - 1. sin(2. a))

                           2    / 4\
          + 0.01293857500 x  + O\x /
k3 := series(s3, x);
       4.247636750 x     + (-0.5773502693 - 1. sin(2. a))

                           2    / 4\
          + 0.01307916821 x  + O\x /
p1 := convert(k1, polynom);
  4.769067180                                                2
  ----------- - 0.5773502688 - 1. sin(2. a) + 0.01164914507 x 
p2 := convert(k2, polynom);
  4.293792441                                                2
  ----------- - 0.5773502692 - 1. sin(2. a) + 0.01293857500 x 
p3 := convert(k3, polynom);
  4.247636750                                                2
  ----------- - 0.5773502693 - 1. sin(2. a) + 0.01307916821 x 
with(plots, implicitplot);
implicitplot([p1, p2, p3], a = 0 .. (1/2)*Pi, x = 0 .. 3, style = point, symbol = [box, cross, circle], labeldirections = [horizontal, vertical], labels = ["&theta; (radians)", "Switching  field (T)"], color = [blue, red, black], symbolsize = 13);

2 3 4 5 6 7 8 Last Page 4 of 1480