Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I received an unexpected error message when trying to minimize a function: evaluating

returns the error message

Error, (in @) too many levels of recursion

Why am I getting this message?  It's hard for me to see how minimizing a function involves recursion, unless Maple is trying to iteratively approximate a solution.

How to open batch maple output as mws in Maple Player?  I was used to doing this daily at work before retirement in my university office on a campus computer.  There was an "open" command in a file menu of some sort that I could click on, I seem to remember.   Am just fresh afte downloading Maple Player on my iMac and have it installed.  The file from my batch Maple output is form a run of batch Maple on the office machine that I can logon to remotely from home, as I did in this case.  I get rough looking graphic output (using the alphabet style graphics methods).  It would be good to have a more detailed and fine grained graphics to answer some basic questions that I am analyzing.  Thanks . . !!

Hello;
        i am wording on fluid dynamics, in which i can up a system of nonlinear partial differential equation with i am suppose to solve using implicit keller box method. i need an asistance on how to implement this in maple.

 

error module is maple.dll_unloaded

 

i have already called stopmaple(kv);

in

testfunction(string hello){....stopmaple(kv));

 

testfunction(1);

testfunction(2);

when call testfunction again , it got error

Trace  on computing the distinct-degree decomposition of the squarefree polynomial
f=x^17+2x^15+4x^13+x^12+2x^11+2x^10+3x^9+4x^4+3x^3+2x^2+4x belong to F_5[x].tell from the output only how many irreducible factors of degree i the polynomial f has, for all i.

 

Mapleprimes_Integral.mw

I have a question regarding following problem:

assume(a > 0, a < 1, t > 0, Z0 > 0, z > 0)

f1 := proc (z) options operator, arrow; 1/z end proc

proc (z) options operator, arrow; 1/z end proc

(1)

I_1 := int(f1(z)*ln((a*z+1)/(1+z/a)), z = 0 .. Z0); 1; MultiSeries:-asympt(%, Z0, 3)

2*ln(Z0)*ln(a)+(a^2-1)/(a*Z0)-(1/4)*(a^4-1)/(a^2*Z0^2)+O(1/Z0^3)

(2)

Using the representation which should hold for all a>0 and z>0

int(z*exp(t)*(a^2-1)/((exp(t)+a*z)*(exp(t)*a+z)), t = 0 .. infinity); 1; combine(%)

ln((a*z+1)*a/(z+a))

(3)

I'm calculating the result the other way around

int(z*exp(t)*(a^2-1)*f1(z)/((exp(t)+a*z)*(exp(t)*a+z)), z = 0 .. Z0); 1; I_2 := int(%, t = 0 .. infinity); 1; MultiSeries:-asympt(%, Z0, 3)

2*ln(Z0)*ln(a)+(a^2-1)/(a*Z0)-(1/4)*(a^4-1)/(a^2*Z0^2)+O(1/Z0^3)

(4)

plot(eval([I_1, I_2], a = 1/2), Z0 = 0 .. 10)

 

So the results are the same.

But if I calculate this with another function

f2 := proc (z) options operator, arrow; 1/(z*(z+a)) end proc

proc (z) options operator, arrow; 1/(z*(z+a)) end proc

(5)

I_3 := int(f2(z)*ln((a*z+1)/(1+z/a)), z = 0 .. Z0); 1; MultiSeries:-asympt(%, Z0, 3)

-dilog(a^2)/a-2*ln(a)/Z0+(1/2)*(2*ln(a)*a^2+a^2-1)/(a*Z0^2)+O(1/Z0^3)

(6)

int(z*exp(t)*(a^2-1)*f2(z)/((exp(t)+z*a)*(exp(t)*a+z)), z = 0 .. Z0); 1; I_4 := IntegrationTools:-Change(int(%, t = 0 .. infinity), t = ln(z)); 1; MultiSeries:-asympt(%, Z0, 3); 1; simplify(convert(convert(MultiSeries:-series(I_4, Z0, 1), polynom), polynom))

ln(a)*(ln(Z0)+ln(a))/a

(7)

I get another result :-/ The Integral doesn't even vanish in the limit Z0 -> 0

Though if I take the limit prior:

int(z*exp(t)*(a^2-1)*f2(z)/((exp(t)+z*a)*(exp(t)*a+z)), z = 0 .. infinity);

-dilog(a^2)/a

(8)

the result is correct. What is the problem here?



Download Mapleprimes_Integral.mw

 

Hi,

I want to compute a formula which is too complicated and it contains some variables. So I divide it into several parts. But it always turns out kernel connection has been lost. I looked maplesoft online help system and change the ConnectionType from 0 to 2. But it does not work. So how does this happen and how to solve?

I attach my maple file which appears error.
Thanks a lot!

Regards,

Yan

with(Groebner):
with(LinearAlgebra):
T := lexdeg([x1,x2,x3],[e1,e2,e3]);
hello1 := proc(xx,yy)
return MatrixMatrixMultiply(xx,yy);
end proc:
hello2 := proc(xx,yy)
return xx+yy- MatrixMatrixMultiply(xx,yy);
end proc:
m1 := Matrix(3, 3, {(1, 1) = -.737663975994461+0.*I, (1, 2) = -.588973463383001+0.*I, (1, 3) = .330094104689369+0.*I, (2, 1) = -.588012653178741+0.*I, (2, 2) = .320157823261769+0.*I, (2, 3) = -.742792089286083+0.*I, (3, 1) = -.331802619371428+0.*I, (3, 2) = .742030476217061+0.*I, (3, 3) = .582492741708719+0.*I});
m2 := Matrix(3, 3, {(1, 1) = -.742269137704830+0.*I, (1, 2) = -.590598631673326+0.*I, (1, 3) = .316590877121441+0.*I, (2, 1) = -.593533033362923+0.*I, (2, 2) = .360143915024171+0.*I, (2, 3) = -.719732518911068+0.*I, (3, 1) = -.311054762892221+0.*I, (3, 2) = .722142379823161+0.*I, (3, 3) = .617863510611693+0.*I});
m3 := Matrix(3, 3, {(1, 1) = -.751491355856820+0.*I, (1, 2) = -.574908634018322+0.*I, (1, 3) = .323636840615627+0.*I, (2, 1) = -.575794245520782+0.*I, (2, 2) = .332066412772496+0.*I, (2, 3) = -.747123071744916+0.*I, (3, 1) = -.322058579916187+0.*I, (3, 2) = .747804760642505+0.*I, (3, 3) = .580574121936877+0.*I});
AA := hello1(m1, m2);
BB := hello2(m1, m2);
GB := Basis([e1- AA,e2- BB],T):
NormalForm(m3, GB, T);

i have got alot of mixed and high degree derivatives. For example:

u[x]*u[x,t]*eta[x,t]+u[]^2*u[x]*eta[x]+kis(x,y)u[x,t]^2*u[]+eta(x, y)*u[]*u[x]^2+ksi[x,t]*u[x]^2*u[x,t]+......

like this alot of terms

my question is how can i solve divided by the derivative of the u(x,t) partial differential equations system and so  how can i find eta(x,t,u) and ksi(x,t,u) 

Hi MaplePrimers!

I have a simulation in MapleSIM, exported as a compiled procedure in maple using -LinkModel(), and -GetCompiledProc.

I'm trying to do parameter estimation on my MapleSIM model.  Within a optimization scheme, I call the MapleSIM model, and it will output a curve.  Using a least squares method, I compare this measurements to synthetic experimental data (I know the actual values), and generate an objective function.  The optimization algorithm will try different parameter values, and try to minimze the objective function.  When the curves are exactly the same, the objective function will be zero.

The problem I am having is certain parameter sets will cause the model to require very small steps.  I wish to put a timeout on these experiments, because speed is important.  However, I would also like to see the results up to the point of requiring very small steps.  For timeout, I was using code along the lines of:

out:= timelimit(30,cProc(params = PData)); #simulate with 30s limit

where PData are the parameter guessses, and cProc is the compiled MapleSim model.

I would like 'out' to be assigned whatever the results were after 30 seconds, even if the model had not finished integrating.

 

Thanks in advance for any help!

I want to define a band matrix using n vectors. Detial given maple sheet. please help me

Help.mw

I want to define a matrix of M by M which have some condition. detail given in maple file. please help me.

Help_Maple-1.mw

I wan to define a square matrix for M and r detial given in Maple sheet. Please Help....Help_Maple.mwHelp_Maple.mw

FirstEigenVector := Matrix(3, 1, {(1, 1) = -.736895432967255+0.*I, (2, 1) = -.588906969844997+0.*I, (3, 1) = -.331924240964690+0.*I});
SecondEigenVector := Matrix(3, 1, {(1, 1) = -.589856901397123+0.*I, (2, 1) = .320280857681335+0.*I, (3, 1) = .741275257969058+0.*I});
ThirdEigenVector := Matrix(3, 1, {(1, 1) = .330233185410229+0.*I, (2, 1) = -.742030156443046+0.*I, (3, 1) = .583384341736151+0.*I});
LHS := ProjOfEigenVector;
LHS := Matrix(3, 3, {(1, 1) = -.736895432967255+0.*I, (1, 2) = -.589856901397123+0.*I, (1, 3) = .330233185410229+0.*I, (2, 1) = -.588906969844997+0.*I, (2, 2) = .320280857681335+0.*I, (2, 3) = -.742030156443046+0.*I, (3, 1) = -.331924240964690+0.*I, (3, 2) = .741275257969058+0.*I, (3, 3) = .583384341736151+0.*I});
RHS := c1*FirstEigenVector+c2*SecondEigenVector+c3*ThirdEigenVector;
RHS := Matrix(3, 1, {(1, 1) = (-.736895432967255+0.*I)*c1+(-.589856901397123+0.*I)*c2+(.330233185410229+0.*I)*c3, (2, 1) = (-.588906969844997+0.*I)*c1+(.320280857681335+0.*I)*c2+(-.742030156443046+0.*I)*c3, (3, 1) = (-.331924240964690+0.*I)*c1+(.741275257969058+0.*I)*c2+(.583384341736151+0.*I)*c3});
solve([LHS[1][1] = RHS[1][1], LHS[2][2] = RHS[2][1], c1^2+c2^2+c3^2 = 1], [c1, c2, c3]);

 

after calculated the projection matrix, 

it is a 3*3 matrix on left hand side

however, combination of eigenvectors on right hand side is 3*1 matrix

when calculated c1,c2,c3 under the condition c1^2+c2^2+c3^2 = 1

how to know whether LHS[1][1] = RHS[1][1], or LHS[1][2] = RHS[1][1] or

LHS[1][3] = RHS[1][1]

First 228 229 230 231 232 233 234 Last Page 230 of 2230