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:

prcLeastSquareLA:=proc(data,degree)

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 
 setShortestDistNode
 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 

can...

I adopted the approach suggested in an earlier post (http://www.mapleprimes.com

/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:

A[mu,1]:=A[mu]

A[mu,2]:=B[mu]

.

.

.

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:

http://hopf.chem.brandeis.edu/yanglingfa/pattern/bruss/index.html

Thank you

Kyle

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 

http://www.maplesoft.com/support/help/Maple/view...

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

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

restart:with(plots):

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

Kyle

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;

return

this function is saved as maple_fun.m.

in Maple,

restart:
with(Matlab):
openlink():
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?

Hi,

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

please help me.it 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

theta[i](x):=subs(m=i,dtheta[m](x)):
theta[i](x):=value(%):
eval(-theta[i](x),x=0);