Maple 14 Questions and Posts Maple 14 Questions and Posts Feed

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


I'm trying to solve the following non-linear ODE numerically:

by ececuting

but maple gives me this error-message:

"Error, (in dsolve/numeric/make_proc) Could not convert to an explicit first order system due to 'RootOf'"

I couldnt find any useful information in the manual. What does this error mean? Is there something wrong with my maple code or is there just no solution for this particulare differential equation?


Thanks in advance


I'm trying to fit some parameters to the data of multiple slightly different experiments.

I've written a function which returns the sums of the error for all experiments.


The function which has to be minimized does the following:

1) Set parameters for the model ODE based on the input parameters

2) calculate the difference (numeric solution of ODE <-> experimentel data)

3) repeat step 1+2 for all experiments

4) return the sum of all differences


The function works as aspected. But when I try to minimize it by calling:

I got an error:

"Error, (in dsolve/numeric/process_parameters) parameter values must evaluate to numeric, got A = A"

For some reason Maple isnt able to set the new parameters of the ODE.

Anyone got an idea how to fix this?


Thanks in advance.





Please find my Maple-worksheet attached below.


How do I solve this problem? I can not drawn it.

I would really appriciate if somebody could help me. 

Hi there,

I'm doing my first steps with Maple and am a little confused with a problem:
I want to solve something like

and expect the result to be z=-m*a. However I get empty brackets as result.
I tried to dig a little and I find that maple has problems isolating x1+x2+x3+x4

Error, (in isolate) x1*a+x2*a+x3*a+x4*a does not contain x1+x2+x3+x4


When I isolate it myself Maple even tells me they wouldn't be the same


I even reduced the problem to x1*a+x2*a and Maple still has problems with that. As this is a really simple problem I strongly assume the issue lies with the user. Can anyone help me there?




How to find the coefficients of the terms diff(eta1(xi1),xi1)*diff(eta2(xi2),xi2)eta1(xi1)^2,...

Dear all

I want your help in a problem.I want to substitute in a summarization ,the variables for Vi and ξi.Because the number of variables are more than 200 I want to know if there is a way to take them from the excel sheet and put the to my program.attached are the maple worksheet and the variables table

P.s. i have done it for 44 variables



I am using densityplot function to produce different plots of a function. The function depends on two variables q and yi .

densityplot(f, q = 0 .. 1, yi = 0 .. 1, grid = [100, 100], style = patchnogrid, color = "Red", axes = boxed, axesfont = [Times, Roman, 15], labelfont = [Times, Roman, 15])

I pass the command style = patchnogrid, but still I get these white lines which I believe are grid lines. 
In the example provided in maple tutorial, when they used this command, all gridlines were removed. Why are they not removed in my plots?

if the plot is not visible, click the link



I don't manage to animate a DensityPlot. Here is my program :

> restart;
> with(Statistics): with(plots):
> P := lambda->RandomVariable(Poisson(lambda)):
> animate(DensityPlot,[P(n),range=0..10],n=1..5,frames=5);

And here is the error message :

Plotting error, non-numeric vertex definition

Thanks for any help


I am new in Maple and I'm trying to get functions from system of equations.

Constants are defined in line 4 and equations are:

eq1 := E2 = fE2(1+(KaE2+Ca)/(1+KaE2*fE2+KaT*fT+KaDHT*fDHT)+KsE2*Cshbg/(1+KsE2*fE2+KsT*fT+KsDHT*fDHT))

eq2 := T = fT(1+KaT*Ca/(1+KaE2*fE2+KaT*fT+KaDHT*fDHT)+KsT*Cshbg/(1+KsE2*fE2+KsT*fT+KsDHT*fDHT))

eq3 := DHT = fDHT(1+KaDHT*Ca/(1+KaE2*fE2+KaT*fT+KaDHT*fDHT)+KsDHT*Cshbg/(1+KsE2*fE2+KsT*fT+KsDHT*fDHT))

KsT = 0.10e11; KaT = 4.6*0.10e6; KsE2 = 3.14*0.10e10; KaE2 = 4.21*0.10e6; KsDHT = 3*0.10e6; KaDHT = 3.5*0.10e6;

fT, fE2 and fDHT are variables, not functions (i.e. fT is not f(T) ) and I am trying to get fT=f(E2,T,DHT,Ca,Cshbg), fE2=f(E2,T,DHT,Ca,Cshbg) and fDHT=f(E2,T,DHT,Ca,Cshbg).

When I type:

eliminate({eq1, eq2, eq3}, {fE2, fT, fDHT})

Maple gives me a blank field. No error, no other comment.

I have no idea where I'm making mistakes.

Any suggestion is appreciated.


Thanks in advance.


Is it possible to solve linear matrix inequality in Maple?
For example, using Matlab and yalmip we can easily
solve Lyapunov equation
    A'P + P * A <0, P> 0

Is it possible to solve this equation in the Maple?



Hi i have two questions about paths in graph package: 1. suppose we have: path1=[1,3,5,7,9] path2=[1,2,3,6,7,8,9] if we want to create a graph that it's edges are edges of path1 and edges of path2,how can we do this? edge set of our graph should be {[1,3],[3,5],[5,7],[7,9],[1,2],[2,3],[3,6],[6,7],[7,8],[8,9]} 2. suppose p is an arbitrary path on a given weighted graph G how can i calculate weight of p (that is sum of it's edges weights) in maple ? thanks very much for your help


RTABLE(151446624,MATRIX([[tf[1,1] = 24110714285.0000/(1000000000000*s+1071534859*218^(1/2))], [tf[2,1] = .843883770817266e20*218^(1/2)/(2500000000000000000000000*s^2+4870187206000000000000*218^(1/2)*s+487482532797170793867)], [tf[1,2] = .841458646091250e19*218^(1/2)/(125000000000000000000000*s^2+350952468875000000000*218^(1/2)*s+50692506828944068713)], [tf[2,2] = 87390625000.0000/(5000000000000*s+3705223967*218^(1/2))]]),Vector[column])


How can I read the individual elements tf[1,1], tf[2,1], tf[2,1], tf[2,2] form this RTABLE?


i have 2 questions:

Question 1. suppose we have generated a random weighted graph with 100 vertices and 2300 edges and found shortest path frome node 1 to node 100 using function ShortestPath(G,1,100). this function returns a path like [1 3 8 2 9 100] but don't get us value of shortest path!

in other words if weights are costs of travelling on edges,we want to find minimum cost of travelling from node 1 to node 100.

how can i find value of shortest path (minimum weight/cost) ?

Question 2. i want to create a sub graph of G by removing the edge that has maximum weight on shortest path found in Question 1. how can i find such edge and how to create such sub graph?


I want to create a NxN circulant matrix of a system with fixed boundaries. This is an idea of what I'm trying to do (but it doesn't work):

> N := 8;

> A := Matrix(N, shape = Circulant[[2, VectorCalculus[`-`](1), ZeroVector[row](VectorCalculus[`+`](N, VectorCalculus[`-`](3))), VectorCalculus[`-`](1)]]);



In other words, I want to make a circulant matrix of the elements 2, -1, 0, 0, ..., -1 (with N-3 zeros). Then I want to set the top right and bottom left elements to zero (that's the fixed boundary portion). The first command almost works for N=4 -- except with extra square brackets in the matrix -- but doesn't work for N=8 or 16 (N is always a power of 2). The second and third commands seem to be bad syntax (I haven't really tried to figure them out yet).

hi dears. i cant solve this question by Analytical method such variational itteration method (VIM) or homotopy perturbation method (HPM), homotopy analysis method(HAM)

1 2 3 4 5 6 7 Last Page 1 of 34