Maple 16 Questions and Posts

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

Hi, I need to update a procedure prcLeastSquareLA so that it includes error estimate information as output in it. So far for the procedure I have this:


local vars,y,A,V,k,e,i,j,v,c,lstVars;  

vars := seq(c_i, i = 0 .. degree);

y := unapply(`+`(seq((c_i)*t^i, i = 0 .. degree)), t);

 for k to nops(data) do 

e_k := y(data[k][1]) = data[k][2];

v_k := [coeffs(lhs(e_k))];

How to extract affine terms from a differential equation after linearization?

Vanderpol Eqn:




To be precise, I linearized...

Modify the shortest path problem by rewriting a procedure Path, non-Recursion
that print out the shortest path from a to e base on the information from
input: destination "e" and setShortestDistNode
output: path : a -> b -> d -> e the total cost: 9
hint: copy everything that we did in shortest path problem, and rewriting a procedure Path, with while loop

I adopted the approach suggested in an earlier post (

/questions/137320-Plotting-A-System-Of-Ode) for the known results,

but some part of curve is missing.

The output I got is

Can I create an tensor in the Physics package in Maple 16 with a free index like i=1..N:







and so on? The index is not one of the type gauge, spacetime or spinor.

Can you expalin the difference between the gaugeindices, spinorindices,and spacetimeindices in the Physics package in Maple 16? 

Hi, there, maple built-in help just has an example of Brusselator model in 1D case. But how can I intepreter the diffusion terms that Brusselator has in 2D case? Could you please give me an example of the code to run Brusselator simulation in 2D case?

Brusselator model:


Thank you


Hello every one,
I got an error while trying to solve numerically a nonlinear system of odes using 
the dsolve command. This error is very common,
"Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging".
I tried the different technique explained in

This question was asked before but because of the curiosity, I 

bring it in the light again. We have a system of odes


a := 1; b := .5; d := 1; omega := .4; h1 := 1+a*cos(x); h2 := -d-b*cos(x+omega);

F := Q-1-d;

de:={alpha*(diff(f(y), y, y, y, y))+G*(diff(theta(y), y, y))+B*(diff(phi(y), y, y))

+6*beta*(diff(f(y), y, y))*(diff(f(y), y, y, y))^2+3*beta*(diff(f(y), y, y, y, y))*

(diff(f(y), y, y))^2 = 0,

In matlab, a function is returned a structure H, say H.a, H.b, ...

In maple, I want to call this matlab function, but what returned to maple is a Record:


Hs:= getvar("H");

Hs:= Record(a=..., b=..., ...);


So how can I convert this record to a normal maple array? so I can access to each element by Hs[1]. 


Thank you


Hi there, I'm trying to pass some calculation to Matlab then retrive the result by getvar(), but seems something is getting wrong.

my matlab function is very simple:


function out =maple_fun(a, b)

out = a + b;



this function is saved as maple_fun.m.


in Maple,

DirectoryLocation := " cd  '/Users/Kyle/Documents/MATLAB/'  ":

I am playing around with the PolynomialFit command in maple.
If I want to get a return of 1% then I can simply plug in the return at time t
for v in the expression below. That works for all different returns.

When I use ImportMatrix to import a MATLAB mat file, I get a list of names and values. Should I write a "for" loop to assign them in Maple or there is an easier way to do it?


I want to parallel the following is good that part1 solved by one thread and part2 solved by other thread.

please help must be noted that part1 and part 2 are independent "at each" iteration.


dtheta[m](x):= theta[m-1](x)+phi[m-1](x);
dphi[m](x):= theta[m-1](x)-2*phi[m-1](x);

for i from 2 by 1 to 10 do;

# part1


3 4 5 6 7 8 9 Last Page 5 of 40