I'm attempting to plot several solutions of this differential equation (I have uploaded my worksheet). I have used this series of commands before without issue, but for some reason I keep getting the error message: "Error, (in plot) incorrect first argument" ect.. Does anyone have any insight into what might be going wrong? Thank you.

Download ass_1_#9.mw

ass_1_#9.mw

f := x1^2/((x2^3)*(x3^2));

would like to extract 1/x2 from 1/(x2^3) and 1/x3 from 1/(x3^2)

and then assign 1/x2 to x2b and assign 1/x3 to x3b

and then replace 1/(x2^3) with x2b^3 and 1/(x3^2) with x3b^2 in f

only in this example has x2 and x3, in real case, it is unknown how many variables all division part

how to automatically to do above for equations with division

Respected sir,

     Maple package is very useful to all the Mathematics students to work with. however, i came to know how to solve

system coupled non-linear ODEs with boundary condtions in Maple 12.0. But my doubt is, the tool 'Boundary value problem solver' which methods and submethods are used to solve those BVPs. Sir please, clear my doubt. your response is highly helpfull for me. Please kindly do the need full.

Thanking you.

of 1/(x^(2*sin(x)^2)+x^(2*cos(x)^2)) over the interval 1..infinity? Is it possible to find that with Maple?

 119^(30^(8^33)) mod 630 with Maple? The WolframAlpha gives the answer 91, but Mathematica 9.0.1.0 fails with it.

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 trying to fit my data (in the form of two vectors "temp" and "resistance") on  resitivity as a function of the temperature "T" using a Dabye model.

The fit, intended to give as a result the two parameters "alpha" and "theta", doesn' seem to give any results.

Here my sample code:

>with(plots):
with(Statistics):

>rho[0]=1.33;

>T[0]=6;

> NonLinearFit(rho[0]+alpha*(T-T[0])*((T-T[0])/theta)^5*(int(x^5/(exp(x)-1)(1-exp(-x)), x = 0 .. theta/(T-T[0]))),temp, resistance, T,output=[parametervalues, residuals],parameterranges=[alpha=0.001..0.01, theta=10..100],initialvalues=[alpha=0.0027,theta=50]);

I am looking forward for your suggestions.

Thanks in advance

Hello, everyone. I have some problem with multithreaded calculation. I just need calculate eigenvalues of matrix m at various parameters (and then export to a file) using advantages of the parallelizing. The following code works but in serial way

restart: with(LinearAlgebra):

m:=ImportMatrix(cat(currentdir(),"m.txt")): # here is matrix m.txt

step:=1:

prc:=proc(k1,k2)

local u,i,j,nmc:

u:=Matrix(ceil(1+(k2-k1)/step),5):

u[1,1]:=k1:

for i from 2 to op([1,1],u) do

u[i,1]:=u[i-1,1]+step:

end do:

for i from 1 to op([1,1],u) do

nmc:=sort(Eigenvalues(m*u[i,1], output='list')):

for j from 2 to op([1,2],u) do

u[i,j]:=nmc[j-1]:

end do:

writedata[APPEND](cat(currentdir(),"u_",convert(k2,string),".txt"), [convert(Re(u[i]),list)]):

print(u[i,1]);

end do:

return finished:

end proc:

with(Threads[Task]):

Start(ArrayTools[Concatenate], 2, Task=[prc,1,20], Task=[prc,20+step,40]);

quit:

The Start(ArrayTools[Concatenate], 2, Task=[prc,1,20], Task=[prc,20+step,40]) function makes two tasks of calculation at the parameter ranges of 1-20 and 21-40. But in this case Start spends twice more time than simply prc(20+step,40). How to realize a multithreaded calculation?

By the way I don't need to use a Concatenate function in Start but without any procedure Start doesn't work.

The following code

restart:
A := Matrix(2,(i,j) -> a||i||j):
A + 1;
A + a;

produces the following output:

Why does Maple accept these additions, which mathematically are nonsensical? And why, if accepting these mathematically, does Maple not evaluate them analogously by putting the scalars either inside or outside the matrix in both cases?

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

Maple people:

The title is the question.  I would like to know if there is any Maple code available, either part of Maple itself or written by a user, to do computations with quaternions or dual quaternions.

I could Google this and probably find something, but I'll probably find a more helpful and less outdated answer here.

I am teaching a student about the subject and I'd like something to help me teach and help him learn.

I have a lot of experience with Maple but I am not a "computer person", so if the code involves fancy "libraries" or something beyond regular Maple worksheets, I may need a few tips how to use it.

GS

Hello everybody,

This one is an easy one but I cannot find the answers.  Suppose that I clik on the template for evaluating a function.  But mine have two (or more) variables (f=f(x,y)).  I have the x=a as one placeholder.  But I need to have an another inplut placeholder for y=b.  Can somenone tell me how to do it under Windows and under Linux.

Thank you in advance for your help.

--------------------------------------
Mario Lemelin

Maple 17 Ubuntu 13.10 - 64 bits
Maple 17 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

ifactor(op(1, convert(0.999987406876435, fraction)));
ifactor(op(2, convert(0.999987406876435, fraction)));

ifactor(op(1, convert(0.999919848203811, fraction)));
ifactor(op(2, convert(0.999919848203811, fraction)));

x1:=2;
x2:=3;
x3:=7;
x4:=17;
x5:=173;
x6:=709;
x7:=5347;
x8:=18713;

i think that there is a need to distinct all factor into a list of x1,x2,x3.... depending on all factors in these two decimal

0.999987406876435 = x4*x2^3*x5/(x3*x1^4*x6)
0.999919848203811 = x1*x8/(x3*x7);

then replace factors with x1,x2,x3....

with(Groebner):
with(LinearAlgebra):
T := lexdeg([x1,x2,x3],[e1,e2,e3]);
GB := Basis([e1+.999987406876435, e2-.999919848203811],T):

Error, (in LinearAlgebra:-Basis) invalid input: LinearAlgebra:-Basis expects its 1st argument, V, to be of type {Vector, {set(Vector), list(Vector)}} but received [e1+HFloat(0.9999874068764352), e2-HFloat(0.9999198482038109)]

with(Groebner):
with(LinearAlgebra):
T := lexdeg([x1,x2,x3,x4,x5,x6,x7,x8],[e1,e2]);
GB := Basis([e1-x4*x2^3*x5/(x3*x1^4*x6), e2-x1*x8/(x3*x7)],T):

Error, (in LinearAlgebra:-Basis) invalid input: LinearAlgebra:-Basis expects its 1st argument, V, to be of type {Vector, {set(Vector), list(Vector)}} but received [e1-x4*x2^3*x5/(x3*x1^4*x6), e2-x1*x8/(x3*x7)]