I have a Matrix 'A' in maple with symbolic entries and I need to use that Matrix in matlab M file. The symbols with in 'A' are already initialized in M file.

How is this possible?

For Matrices 'A' of small dimension, I can use the codegeneration package to convert 'A' to Matlab Matrix and then copy paste.

But for matrices of 1000x1000, lets say, does not seem as a option.

Hello. 
I have a workcheet, that is incompleete, when using

  > with(XMLTools);
  > ParseFile("FILE");

i get the wollowing error, but i cant seem to fix it by my selv in an xml editor. 
Here you have the file:

9.mw

Hope you can help

Hello,

I have two sets of data representing two function that depends on x with a parameter A.

I need to do a fit on both data series at the same time so to fit with the best parameter A.

Here is how I do a fit on one function

> f(x):=A*cos(x-B)^(2);               

> g(x):=A*(cos(x-C)^(2)+ sin(x-C)^(2))^2; 

> fit1 := Fit(f(x), r, x, parameternames = [A, B, C...

hey friends..
Am trying to solve 2d poission equation with finite difference method and using different scheme is got

u_i-1,j -8u_i,j + u_i+1,j + 3u_i,j-1 +3u_i,j+1 = -h *16

and i have made code for this which is

> restart;
> with(LinearAlgebra);
> printf(`The lefthand endpoint on the X-axis is 0.\n`);
> OK := FALSE;
> while OK = FALSE do printf(`Input...

r-x^4
u-(x^3)*y
v-x*y^3
w-y^4

Step 1
with(Groebner):
K := {r-x^4,u-(x^3)*y,v-x*y^3,w-y^4};
G := Basis(K, 'tord', degrevlex(r,u,v,w));

Step 2
G intersect k[r,u,v,w]

G2 := Basis({w*G[1],w*G[2],w*G[3],w*G[4],w*G[5],w*G[6],w*G[7],w*G[8],w*G[9],w*G[10],w*G[11],w*G[12],w*G[13],w*G[14], (1-w)*K[1], (1-w)*K[2], (1-w)*K[3], (1-w)*K[4]}, 'tord', degrevlex(r,u,v,w));

G2 has too much though include the correct solution

correct solution...

3*rho1 - 2*rho2 + rho3 - rho4 = -1

4*rho1 +   rho2 - rho3        = 5

original without cost function:

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,(y1^1000)*(y2^1)*(y3^1)*(y4^100)- x1*x2*w + 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));

after have cost function 1000*rho1 + rho2...

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,x1*x2*w - 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));
(1,1,0,2)

subs({rho1=1,rho2=1,rho3=0,rho4=2},3*rho1 - 2*rho2 + rho3 - rho4);
subs({rho1=1,rho2=1,rho3=0,rho4=2},4*rho1 +   rho2 - rho3       );

can also reduce with G8,G7

R1 := Reduce((x2^(5+1))*(w^1), {G[7...

Here is my sheet. non_linear_P2_last_q.mw

In this problem I have achieved every thing what I aimed, thnaks to the wonderful people here!!!

I am looking for a way to do it in a more structured way meaning the calculation of Jacobian in the attached sheet is somewhat manuel, though I wonder does maple has some other way?

Hi there i have function as:

f(B,F)=

where: F,B>0 and 

,

such that: 

h=B^2;

I need to plot the surface (contour plot) f(B,F) at 0 and -ve points i.e f(B,F)<0.

Thanks

Hi all I have the following ODE

ODE_1:=diff(w(r),r)+R*(diff(w(r), r))^3=K/r;

ODE_T_1:=collect(algsubs(w(r)=u(r)+(r-1)*(U-0)/(delta-1), ODE_1),diff(u(r), r)) ;

 and when I try to do the following integration,

eq1:=int(phi[i](r)*ODE_T_1,r=1..delta) assuming delta > 1;

it gives the following error

"

Error, (in assuming) when calling 'int'. Received: 'wrong number (or type) of arguments:

i am doing a Q Laplace and finally summation it

however, i find an AccurateQSummation which is not like Summation

Should i use it, if not, use SumTools' summation how to simplify it?

restart;
with(QDifferenceEquations):
qexp := QPochhammer(-(-p*t), q, infinity);
f:= qexp*t;
b := subs(t=q^n,f);
QSimplify((1-q)*AccurateQSummation(b*q^n,n=-infinity..infinity));

with(SumTools):
QSimplify((1-q)*Summation(b*q^n,n=-infinity..infinity));

Is it possible for someone to provide the roller coaster maplesim model?

Hi

i know to make an asymptotic expansion one uses asympt(1+x,x); how do you make it go to -infity rather than plus infinity?

Hi,

I'm trying to modify the size of the window where an animation is shown. Tipically, this window has axis going from xmin to xmax and ymin to ymax but this causes some scaling problems when xmax-xmin>>ymax-ymin. I would like to focus in a certain region of my animation and change this viewpoint depending on the position of the elements of the animation. Can anybody help me?

This is my code:

> tf := 5;
> nfr := 201;

I have a huge PDE system included reaction-diffusion system with 4 independent and 12 independent variables. I read on the maple help for pdsolve/numeric/method. I found it is only for single partial differential equations in two independent variables. Is there anyone knows the ability of Maple for such a system?