mehdibgh

240 Reputation

6 Badges

7 years, 329 days

MaplePrimes Activity


These are questions asked by mehdibgh

As you know the matrix calculations of the Matlab is more powerfull than Maple, My question is that is it possible to call Matlab solver from Maple worksheet? I mean, in part of my Maple worksheet (for loop) I need to compute the eigenvalue and square root of a big matrix (order of 1000) that Maple couldnt while Matlab can do it. So it will be so nice if I can call Matlab directly from my ws rather than I export my data and import it to Matlab and calculate and return back to Maple (manualy It is tedious).

Is it possible???

How do you recommend to calculate the square root of big Matrices (e.g, 300*300) with Maple??

My machine couldnt calculate the square root of Matrices (9*9) as you see below:


 

``

restart

Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received shareman

 

with(LinearAlgebra):

``

A := Matrix([[1, 2, 3, 4, 5, 6, 7, 8, 9], [9, 8, 7, 6, 5, 4, 3, 2, 1], [1, 2, 3, 4, 5, 6, 7, 8, 9], [9, 8, 7, 6, 5, 4, 3, 2, 1], [1, 2, 3, 4, 5, 6, 7, 8, 9], [9, 8, 7, 6, 5, 4, 3, 2, 1], [1, 2, 3, 4, 5, 6, 7, 8, 9], [9, 8, 7, 6, 5, 4, 3, 2, 1], [1, 2, 3, 4, 5, 6, 7, 8, 9]])

A := Matrix(9, 9, {(1, 1) = 1, (1, 2) = 2, (1, 3) = 3, (1, 4) = 4, (1, 5) = 5, (1, 6) = 6, (1, 7) = 7, (1, 8) = 8, (1, 9) = 9, (2, 1) = 9, (2, 2) = 8, (2, 3) = 7, (2, 4) = 6, (2, 5) = 5, (2, 6) = 4, (2, 7) = 3, (2, 8) = 2, (2, 9) = 1, (3, 1) = 1, (3, 2) = 2, (3, 3) = 3, (3, 4) = 4, (3, 5) = 5, (3, 6) = 6, (3, 7) = 7, (3, 8) = 8, (3, 9) = 9, (4, 1) = 9, (4, 2) = 8, (4, 3) = 7, (4, 4) = 6, (4, 5) = 5, (4, 6) = 4, (4, 7) = 3, (4, 8) = 2, (4, 9) = 1, (5, 1) = 1, (5, 2) = 2, (5, 3) = 3, (5, 4) = 4, (5, 5) = 5, (5, 6) = 6, (5, 7) = 7, (5, 8) = 8, (5, 9) = 9, (6, 1) = 9, (6, 2) = 8, (6, 3) = 7, (6, 4) = 6, (6, 5) = 5, (6, 6) = 4, (6, 7) = 3, (6, 8) = 2, (6, 9) = 1, (7, 1) = 1, (7, 2) = 2, (7, 3) = 3, (7, 4) = 4, (7, 5) = 5, (7, 6) = 6, (7, 7) = 7, (7, 8) = 8, (7, 9) = 9, (8, 1) = 9, (8, 2) = 8, (8, 3) = 7, (8, 4) = 6, (8, 5) = 5, (8, 6) = 4, (8, 7) = 3, (8, 8) = 2, (8, 9) = 1, (9, 1) = 1, (9, 2) = 2, (9, 3) = 3, (9, 4) = 4, (9, 5) = 5, (9, 6) = 6, (9, 7) = 7, (9, 8) = 8, (9, 9) = 9})

(1)

MatrixFunction(A, sqrt(v), v)

Error, (in LinearAlgebra:-MatrixFunction) could not compute finite interpolating value by evaluation of (1/2)/v^(1/2) at eigenvalue 0 which has multiplicity greater than one in the minimal polynomial

 

``

``

``

``

``


 

Download askkk.mw

I have many linear equations as below(f,g,h,...,p, are linear of S,T,..,W):

y1=f(S[i,j],T[i,j],U[i,j],V[i,j],W[i,j]);

y2=g(S[i,j],T[i,j],U[i,j],V[i,j],W[i,j]);

y3=h(S[i,j],T[i,j],U[i,j],V[i,j],W[i,j]);

.

.

.

yn=p(S[i,j],T[i,j],U[i,j],V[i,j],W[i,j]);

Where (i,j)=(0,0),...,(I,J)

How ask Maple to write them in Matrix form as below:

AX=0

Where X is: X=Transpose{S[0,0],S[0,1],S[0,J],...,S[1,0],S[1,1],...,S[1,J],...,S[I,0],S[I,1],...,S[I,J],

                        T[0,0],T[0,1],T[0,J],...,T[1,0],T[1,1],...,T[1,J],...,T[I,0],T[I,1],...,T[I,J],...,W[I,J]}

    

I have an equation were build with many unfixed numbers of variables as below:

eq:=f(U[i,j],V[i,j],W[i,j],S[i,j],T[i,j]), i,j=0,1,2,...(changing during the program running)

How do you suggest to calculate the first derivative of such equation with respect to U[i,j],V[i,j],W[i,j],S[i,j],T[i,j]?

Suppose that i,j=22, in this way to calculate the metioned derivative with respect to U[0,0] only, I must write:

diff( eq(U[0,0],U[0,1],U[0,2],...U[22,22], V[0,0],V[0,1],...,T[22,22]), U[0,0])  ( More than 2420 terms)

As you see it needs cumbersome writings if i,j be constants, and because the i and j are changing during the program, I dont know how to cope with this derivatives in Maple??

I have a procedure, and I want to use it so many times whthin my equations, because all notation in my equations are in indicial format, it will be so nice if I can call the procedure in indicial form rather than paranthesis form.

Is there anyone knows how is it possible in Maple?

Rr is my procedure that accepts 2 inputs, So there is no problem to call it using Rr(x,y), I am trying to call it with Rr[x,y] (indicial form

)

Is it possible to do so in Maple? It will be so great to do so.

First 14 15 16 17 18 Page 16 of 18