hi friends!
I have a 2*2 matrix(namely Mcell), whose elements vary as a complex function of "nu".
for each cycle of a "for" loop, h[j] := abs(Re(Mcell[2, 2])), whose diagram is periodic and has many crosses with y=1.
i just wanted the first two points(their difference) and used

L := Roots(h[j]-1, nu = 1 .. something, numeric); and
w[j] := L[2]-L[1];

in a "for" loop.

i check the answer, they are all true(w gets 20 values in 20 iterations).

Hello,

After applying "solve" to an equation system with lots of sin, cos etc., the result looked very unreadible for my non-math-eyes...basically, a lot of RootOf() and _Z4 expressions...The equation system has 6 variables rot_x1, rot_y1, rot_z1, rot_x2, rot_y2, rot_z2 and I have solved for rot_x1, rot_y1, rot_z1

...I then applied convert(answer, radical) and it looks like this, which is something I can understand  - a few fractions and root-signs etc...(see pic attached)

Question 1: Is it ok to do this conversion or do I lose some solutions without knowing it??

Question 2:

I then tried to eval the converted answer by assigning rot_x2, rot_y2 and rot_z2 with values....Actually, I expected that the result would be easy-to-read values for rot_x1, rot_y1 and rot_z1, but instead I get  something with roots and "signum" which I dont really understand (esp. the signum)....Therefore, I applied convert(piecewise) to it....Does that make sense what I just did???....Can solutions go lost or change when i apply the "conversion" function?????....

THANKS for any help....and sorry, I know I should better improve my Maths, but I guess there is no time for this atm...;)...

         
 

 

Hello,

I need to solve a system of ode in matrix. I've been trying with the proc instruction for three days, but it doesn't work! (I'm quite new into Maple).

My system as the form:

M*y'' + C*y' + K*y = F.

M,C,K are matrix of (n*n) (n is between 2 and 200 depending of the problem). y is the unknown, a vector of n degrees of freedom depending of time. F is a Vector of n components.

I've try with:

Hi, would anyone mind explaining the use of k:=op(J)[1] in this code below. 'Int((3*polylog(5/2, -exp(b*(u-i)))/b^(5/2)+i*polylog(3/2, -exp(b*(u-i)))/b^(3/2))*r^2, r = 0. .. 10.)'; J:= convert(%, rational); K:='op(J)[1]'; with(codegen); makeproc(K,r): KK:=optimize(%, 'tryhard'); 'Int( KK(r) , r=0.6 .. 10, method = _d01akc, digits=8)'; evalf(%); Thanks MS

Hi, I need to put a 2-dimensional quadric (given as an equation including x's & y's) into normal form. How can I do this? I skimmed through the Geometry/2D-Euclidean directory of the Maple help but didn't find a solution.

Thanks in advance.

Hi,

I am an absolute beginner with Maple and am appreciating any help.

Basically, I want to use maple because it can deal with algebraic expressions. For example with Matrixmultiplication I don't just want to see the final result, but also I want to see the algebra behind each matrix element. 

So I have defined the following 2 matrices :

rot_x := matrix([[1,0,0], [0, cos(alpha), -sin(alpha)], [sin(alpha), cos(alpha), 0]]); 

rot_y := matrix([[cos(alpha), 0, sin(alpha)],[0,1,0],[-sin(alpha), 0, cos(alpha)]]);

Hi everybody, I have a series and I want maple to truncate small variations of a variable of the order greater than n. For example, series1:=1+x+x^2+O(x^3) = 1+x+x^2. Thank you,

hi, i want to solve these 4 equations and i have two problems

1-when i use the function SOLVE , the output is empty brackets, why?

2-how can i add constrains or inequalities to make maple consider them in the solution

the solve line is

Is there a faster way to do the following code? This takes 6.5 hours to run for n:=10...

> with(Statistics):
n := 1:
A := Array(1..10000000):

> for k to 10000000 do
t := RandomTools:-Generate(list(float(range=0..1,digits=20,method=uniform), n)):  # chooses 6 random variables from the interval [0,1]

t := sort(t):   # sort them in ascending order

Hi all,

I am using the Physics package. I have an expression with many terms and need to contract indices with the Simplify command. It returns the following

> KLMN := Simplify(subs(substtot, G4Z2s));
Error, (in Check) unable to handle tensor expressions as 4*I
 

Has anybody come across this problem?

Thanks in advance.

hello, im just trying to do a solve function in  maple 13 for this equation.....1.5=2.4713(1+e^-t+1)*(1-1.23934*e^-0.2145(t+1))^3

I want to solve for T but it keeps telling me that the solutions are lost when i try to solve, please helppppp.  im freaking out here.  and the "e's" in the equation are the exponent function, not another variable.

hello,

1.how can solve the riccati equation for nonlinear time invariant systems, please give example with 4X4 matrices.

2.for the same problem how one can get  the plot for input and states by new method called (theta- D). suboptimal nonlinear control technique.

 

 

Hi all, Is there any way I can check what my maple is doing? It's running this code for almost two days now.

 

Hi!

I need help in solving following 4th order system of 02 ODEs which are coupled

where, '  indicates derivative w.r.t  time t , and A,B are constants.

Hey Everyone!

I cant figure out why Maple can give eigenvectors for a matrix A, through Eigenvectors(A); but I cant use LinearSolve() to do this same calculation. (I get <0.,0.> as the answer).

Here is what I am doing (note, Eigenvectors(S) gives correct answer).

with(LinearAlgebra);

S:=<<7476.45,303.62>|<303.62,26.19>>;

I2:=IdentityMatrix(2);
charPoly:=Determinant(x * I2 - S);
EV:=solve(charPoly=0,x);

         produces: EV := 7488.802926, 13.83707411