MaplePrimes Questions

Good day everyone,

I am having a problem writing a PDE code and it is given an error code

"Error, (in pdsolve/numeric) initial/boundary conditions must be defined at one or two points for each independent variable"

The link is attached below. Please, help. Thanks

I am trying to compute the determinant of a 3 x 3 matrix. It has very involved symbolic entries.

The cpu goes to 100% and after about a minute I get BSOD on windows 10. Using a 4 core i7 64 gig of ram.

Can I stop Maple using all the cores to see if this helps the BSOD.?

I bought a subscription for the Maple calculator on google play and the app still wont give me more than 5 steps a day. Does anyone know what the problem could be?

How can i get mapple code for solving Fractional Patfial Differential Equation with Laplace, Adomian and New Iteration Method


I want to draw a 2d directionfield with arrows for a differential equation system with three equations. I have tried the following: 

mu := 1/(365*80);
alpha := 20;
psi := 100;
varphi := 2;
K := 100;
theta := 0.5;
sigma := 0.001;

beta := 0.3;
des := diff(S(t), t) = (psi + mu)*K - beta*S(t)*In(t) - (mu + varphi + psi)*S(t) + (theta - psi)*V(t) - psi*In(t), diff(V(t), t) = varphi*S(t) - sigma*beta*In(t)*V(t) - (mu + theta)*V(t), diff(In(t), t) = beta*S(t)*In(t) + sigma*beta*In(t)*V(t) - (mu + alpha)*In(t);

phaseportrait([des], [S(t), V(t), In(t)], t = 0 .. 100, [[S(0) = 70, V(0) = 0, In(0) = 30]], V = 0 .. 100, In = 0 .. 100, stepsize = 0.05, arrows = medium, scene = [V(t), In(t)], linecolour = sin(t*Pi/2))

The problem is, Maple plots only the curve for the initial values and not the arrows of the direction field. How can I solve this problem?


How do I use ErrorPlots to put error bars on my scatter plot of data, AT THE LOCATION OF THE SCATTER PLOT DATA POINT? Below I write an example demonstrating how I plot the data and error bars using ErrorPlots.  When I use display to combine plots you clearly see that the errorbars are not on the scatter plot data points.  How do I give ErrorPlots, the x-y locations of the data points, so it knows to print the error bars on the location of the data points?  Is there a better approach to putting horizontal and vertical error bars onto scatter plot data?


I have some complicated calculations which result in a lot of 1e(-7) or even 1e(-12). These numbers greatly complicate the expersisions without adding any benefits. 

How do I instrcut Maple to round the calculations to the nearest decimal precision?

I tried "Digits" and "Rounding" and "interface(precisiondisplay)" but none of them get rid of something like 3.343^(-12)* expr. 

I still see such infistimal terms showing up. 

If you're wondering, I am exporting these symbolic terms to Matlab to be used as optimization functions. So, having these extra terms add unnecessary optimization steps. 


A system of differential equations

But solution I want to express in terms of other specified function

How to dsolve to get a solution that in terms of a list of specified function? For example in terms of sin , cos , Bessel etc ?

Sol := [a(t) = sin(t) + cos(t) , b(t)= sin(t) ....]

I try to eliminate m  in                    X := (a*m^2 + 2*m^2*p + 2*p)/(2*m^2)
                                                      Y := a/(2*m)

eliminate({X, Y}, m);
without results. How can I do ? Thank you.

The following code works in one worksheet but not in another (where it returns the error in the title). What could be the reason?

[> restart:
    P:=display(polygonplot([[0,0],[0,1],[1,1],[1,0]], color=red,axes=none,style = polygon)):
[> exportplot("Test.jpg",P, plotoptions="height=30,width=30,quality=90");

In the worksheet (.mw) where it does work, I can't put the code in one command line unless the file is already there, in which case it overwrites.

- I create a function of 2 variables based on Procedures as follows:

Test := proc (a, b)

local z;

if not is({args}, set(numeric)) then return ('procname')('args') end if;

z := a*x^2+b*y^2-3*x+3*y+3;

minimize(z, location)[1]

end proc:

- The function calculates its value if given the specific values of the variables:

Test(2, 2.5) --> 0.975

- But somehow I can't graph this function with the following syntax:

plot3d(Test(a, b), a = 1.0 .. 5.0, b = 2.0 .. 3.0)

- MAPLe 17 issues the following notice: "Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct". And no graphs are drawn

In the same way, MAPLE 15 can plot the graph. Does anyone know what the problem is with MAPLE 2017 and how to fix it?

Sincerely thank

I noticed Maple changes 1/(4*(x-2)) to 1/(4*x-8)  when copying the output and pasting it back to the worksheet.  

But in the variable it does remain 1/(4*(x-2))  which is what I wanted, as I am going to look for this form later on and I do not want it expanded.

But I am curious why Maple does this. It is a little annoying when copying something from one worksheet to another and have it change.  Here is a small movie and the code to see it. 

Is there a way to prevent this automatic expansion?  Maple 2021.2 on windows 10




I use worksheet interface and have my display settings as the following

Unless I'm missing something where are the closed integral symbols? 

I couldn't find them but you can create them if you wanted to.  Then I tried adding them to a palette with DocumentTools and couldn't get it to work. 

`∮` `∯` `∰`

Creating the double and triple symbols can be done  easily by concatenating as many ∫  as you desire. But can't seem to add them to a palette.

DocumentTools:-AddPalette("Integral Symbols");

AddIcon("Closed Integral",`∮`)

I suppose you need to create a png file for that to work.  However the above would be simpler.

Maybe add functional closed integral operator symbols for the palette? 

how can i factor from 0.03 in this equation and show like the yellow equation



f := 0.03 - 0.03*cos(5.885*t)


l_1 := diff(x1(t), t) = 1/2*x1(t) - x2(t) + cos(4*t);
                 d          1                         
         l_1 := --- x1(t) = - x1(t) - x2(t) + cos(4 t)
                 dt         2                         

l_2 := diff(x2(t), t) = 3/2*x1(t) - 2*x2(t) - 1;
                    d          3                    
            l_2 := --- x2(t) = - x1(t) - 2 x2(t) - 1
                    dt         2                    

dsolve({l_1, l_2}, {x1(t), x2(t)});
    /          4               296              28          
   { x1(t) = - - exp(-t) _C1 + ---- sin(4 t) - ---- cos(4 t)
    \          3               1105            1105         

           /  1  \                                    
      + exp|- - t| _C2 + 2, x2(t) = -2 exp(-t) _C1 + 1
           \  2  /                                    

         93              36                /  1  \    \ 
      - ---- cos(4 t) + ---- sin(4 t) + exp|- - t| _C2 }
        1105            1105               \  2  /    / 

A := <<1/2 | -1>, <3/2 | -2>>;
                           [1            ]
                           [-  &uminus0;1]
                           [2            ]
                      A := [             ]
                           [3            ]
                           [-  &uminus0;2]
                           [2            ]

Eigenvectors(A, output = list);
              [                  [        /[2]\ ]]
              [[-1      /[1]\ ]  [        |[-]| ]]
              [[--, 1, { [ ] }], [-1, 1, < [3] >]]
              [[2       \[1]/ ]  [        |[ ]| ]]
              [                  [        \[1]/ ]]


I assume, I am getting a result in the format of:

Which implies, that the eigenvalues from dsolve are -1/2 and -1, which are correct according to the eigenvalues gotten from the matrix A.

But the eigenvector for the eigenvalue -1/2 appears to be incorretly given by dsolve, I should get 2/3 and no -4/3 in the top axis of the V_lambda(-1). For the second row I am getting -2 and not 1. It appears that Maple has multiplied the eigenvector for the eigenvalue -1 with -2? Or am I missing something?

If I am not missing anything; what appears to cause this issue?

First 53 54 55 56 57 58 59 Last Page 55 of 2141