MaplePrimes Questions

how to solve this case which like sylvester

NullSpace(A*X+X*A) = B

find X?

 any things like 

A*X + X*A = NullSpace^-1(B)

 

if i assume Matrix([[x1,x2,x3],[x4,x5,x6],[x7,x8,x9]]);

find NullSpace(Matrix([[x1,x2,x3],[x4,x5,x6],[x7,x8,x9]]))

then from the equation of it, put all equations = 0, and then solve them

and then substitue value of B into the result, can it be said NullSpace^-1(B)?

Can we define/set a range in Maple. e.g 


I have the following equation:

y = 1.048 + 1.02*x + 6.118*(z-4.041*x^2) + 16.22*(z^2) +6.241* (x*z)

The value of z is within 0.001 - 0.543, y is from 1 - 12 and x is from 0.001 - 0.7

How should I define it in Maple, so while solving equations it read the values within the given range? 


x := (1+v/2*cos(u/2))*cos(u);
y := (1+v/2*cos(u/2))*sin(u);
z := v/2*sin(u/2);

I converted this comment into a question since it was unrelated to the question in which it was put.
That way also other MaplePrimes users can have a shot at this new question.
Preben Alsholm
 

Dear Preben Alsholm 

Hi;

Hope you to be healthy and have nice times,

I have another problem and if it doesn't consuming your valuable times, please guide me.

I have some Basis function,say W1(t),W2(t),W3(t),W4(t) which are orthonormal and i want to write a program that can approximate the integral of W(t)=[W1(t),W2(t),W3(t),W4(t)] again by Wi's, in other word
int(W(t')dt',t'=0..t)≈PW(t), where P is knows as integral operational matrix. the following is my attempt and unfortunately has no real solution!!!!

 

restart:
>              # Definition of 3th B-Spline
>
> piecewise(x>=0 and x<=1,1,0):
> h[1]:=unapply(%,x):
>
>             # Definition of 3th B-Spline
>
> h[2]:=simplify(int(h[1](x-t),t=0..1)):
> hh:=unapply(%,x):
>
>              #Definition of 3th B-Spline

int(hh(x-t),t=0..1): 
> simplify(%):
> N:=unapply(%,x):
>
> J:=1:     # Number of base function is 2^J+2
> phi:=linalg[matrix](2^J+2,1):
> for i from -2 to 2^J-1
> do
> N(2^J*x-i)*h[1](x): #for deleting out side of[0,1]
> simplify(%):
> phi[i+3,1]:=unapply(%,x):
> od:
>
> w[1]:=phi[1,1](x):
> w[1]/sqrt(int(w[1]^2,x=0..1)):
> W[1]:=unapply(%,x):     
>
> for i from 2 to 2^J+2
> do
> kk:=0:
> for j1 from 1 to i-1
> do
> aa[j1]:=int(phi[i,1](x)*w[j1],x=0..1):
> bb[j1]:=int(w[j1]^2,x=0..1):
> kk:=kk+aa[j1]/bb[j1]*w[j1]:
> od:
> w[i]:=simplify(phi[i,1](x)-kk):
> w[i]/sqrt(int(w[i]^2,x=0..1)):
> W[i]:=unapply(%,x):                 #Orthonormality process
>
>
> for j from 1 to 2^J+2 do
> for io from 1 to 2^J+2 do
> f[j]:=int(W[j](s),s=0..x);
> c[io][j]:=int(f[j](x)*W[io](x),x=0..1);
> od;
> od;

 

 

Best Wishes

I am new to Maple and I have difficulties simplifying expressions like the following: (exp(a+b)+exp(a+c))*exp(-a).

I would expect to see exp(b)+exp(c) but nothing happens when I use the simplify() function.

I googled but didn't find a solution.

Hi everyone,

 

At some point in my calculations, I want to do some simple manipulations wiht the use of the «clickable» features in M17.  For the moment, I can see directy how to do it by hand and I do it by copy-paste and use the mous to change it a little bit.  But I just cannot find a simple way to do it by simple click.

 

You will find in this attach question, the maipulations I want to do.  It is starting from equaiton (7).  Am I at the edge where it is simple by hand than with the features?

 Space_manipulation.mw

 

Thank you in advance.

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

Hello,

I was wondering how (or what is the best way) to write a worksheet in which a change of formula is used when a certain value on the y-axis is reached.

So for example: if there is a mass-spring system with damping in it, I would like to change the value of the damping when the displacement/velocity/acceleration has reached a certain value.
So when I apply a force to the mass-spring system, and the acceleration for example is LOWER than 0.2 m/s^2 I use a value of X % damping, but when the value of the acceleration is HIGER than 0.2 m/s^2 I want to apply Y % damping. So in time the curve will increase (when low damping is used) and the curve will decrease (because high damping is used, because the y-value is higher than 0.2 m/s^2), and so on...

I hope somebody has a 'simple' idea. I know what I want to do, but I don't know how to put this down in a formula which I can write in Maple.

Greetings,
Frank

Hi, i need help. I'm currently working with Taylor and Maclaurin series in Maple. I can easily compute the sum by typing in fx : 

taylor( (ex,x=0, 5) , and then I get the first 5 numbers of the series. But I would like Maple to write the series as a sum from n =0 to infinity fx.  I can't figure out how to do it. Can it be done? 

Thanks for helping.

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

Dear guys. I want to solve this equation analytically:

diff(Q(t),t)^2 = ln(t)^(b) + b*ln(t)^(b-1)

I think it is impossible. So I assumed that diff(Q(t),t)^2 = diff(P(t),t) and now I can solve it for P(t) easily. I obtained:

P(t) = t*ln(t)^(b)

Now, I want to know that there is anyway to obtain Q(t) using the latter relation and  diff(Q(t),t)^2 = diff(P(t),t) ?

Also I have a function as V(P)=P^(-2). What can I say about the behaviour of V(Q)?

Thanks a lot.

A common spreadsheet computation is to take values in column A in a spreadsheet and multiply by values in column B, e.g. A1*B1=C1; A2*B2=C2.  How do I do this in a Maple matrix?

I am trying to reduce a tensor expression: ωiωjUi,j 

For which I have tried the following

restart; with(Physics)

Setup(dimension = [3, `+`], coordinatesystems = X, spacetimeindices = lowercaselatin):

`The dimension and signature of the tensor space are set to: [3, +] `

 

`Default differentiation variables for d_, D_ and dAlembertian are: `*{X = (x1, x2, x3)}

 

`Systems of spacetime Coordinates are: `*{X = (x1, x2, x3)}

 

`Defined objects with tensor properties`

(1)

omega[i] := d_[j](U[k](X))*ep_[i, j, k];

Physics:-LeviCivita[i, j, k]*Physics:-d_[`~j`](U[`~k`](X), [X])

 

Physics:-LeviCivita[l, m, n]*Physics:-d_[`~m`](U[`~n`](X), [X])

 

Physics:-LeviCivita[i, j, k]*Physics:-LeviCivita[l, m, n]*Physics:-d_[`~j`](U[`~k`](X), [X])*Physics:-d_[`~m`](U[`~n`](X), [X])*Physics:-d_[`~l`](U[`~i`](X), [X])

 

-(-Physics:-d_[n](U[i](X), [X])*Physics:-d_[k](U[`~n`](X), [X])+Physics:-d_[m](U[i](X), [X])*Physics:-d_[`~m`](U[k](X), [X])-Physics:-d_[i](U[k](X), [X])*Physics:-d_[l](U[`~l`](X), [X])+Physics:-d_[k](U[i](X), [X])*Physics:-d_[l](U[`~l`](X), [X]))*Physics:-d_[`~i`](U[`~k`](X), [X])

(2)

continuity := [D_[l](U[l](X)) = 0]

[Physics:-d_[l](U[`~l`](X), [X]) = 0]

(3)

red_eq := subs(continuity, expr)

-(-Physics:-d_[n](U[i](X), [X])*Physics:-d_[k](U[`~n`](X), [X])+Physics:-d_[m](U[i](X), [X])*Physics:-d_[`~m`](U[k](X), [X]))*Physics:-d_[`~i`](U[`~k`](X), [X])

(4)

Question 1. Am I using the continuity condition correctly? How do I use this condition correctly? If I change the index for the expression, substitution does not work correctly. 

Question 2. How do I expand the red_eq term in terms of the basis to give out the full expression? 

Question 3. I would like to eventually replace U by (A-Amean) in the current expression. How do I implement this? 

Download term8.mwterm8.mw

I am running into what looks like an interference between the Physics package and linear algebra, specifically when using the . (dot operator) to stand for matrix-vector multiplication:

I have a column vector and a matrix and need to evaluate the product matrix.vector. Usually this works fine and gives me the column vector with the result, just as expected. In this particular case (having with(Physics[Vectors]) at the top of the sheet) I get an error message:

Error, (in Typesetting:-delayDotProduct) invalid input: Physics:-Vectors:-`.` expects its 1st argument, a, to be of type Or({algebraic, procedure}, `=`), but received Matrix(3, 3, ...)

The Matrix it complains about is in fact the one I am sending it. It was created from a Physics:-Vectors object using the Component function (so it should be a regular Maple Matrix). In fact, the Matrix is the result of a LinearAlgebra:-MatrixInverse operation.

The weird thing is: This works when I am using 1-d input in Worksheet mode (which I am usually doing). In this case I am creating a sheet in Document mode using 2-d input as it is a sheet I will use in some teaching activity later this year, so I want it to look more polished. I can replace the dot with the relevant MatrixVectorMultiply function and it will work, but that is a kludge I don't think  should be necessary.

Anybody ever seen this? BTW, I am doing this on Maple 15 using the version of Physics that came with it.

Thanks,

Mac Dude

 

First 1473 1474 1475 1476 1477 1478 1479 Last Page 1475 of 2426