Items tagged with vector vector Tagged Items Feed

I have a vector, it is a zero vector apparently. I don't know what opperation lead to it being a zero vector. Any time I try to reduce this thing, I mess it up. I can't reduce it manually.

a=-((4*I)*sqrt(3)+4*sqrt(3)-5-7*I)/((4*I)*sqrt(3)+4*sqrt(3)-7-9*I)

b=-((390*I)*sqrt(3)+30*sqrt(3)-52-675*I)/(-3-4*I+2*sqrt(3)+(2*I)*sqrt(3))^3

c=-(1/2)*sqrt(2-2*I)*((6*I)*sqrt(3)+6*sqrt(3)-11-10*I)*(-2+sqrt(3))*sqrt(-2-2*I)*((2*I)*sqrt(3)+2*sqrt(3)-5-2*I)/(-3-4*I+2*sqrt(3)+(2*I)*sqrt(3))^3

 

<a/sqrt(a^2+b^2+c^2),b/sqrt(a^2+b^2+c^2),c/sqrt(a^2+b^2+c^2)>

 

The simplify command will reduce this bugger no problem. The moment I try rationalizing the denominator or anything like that I end up buggering the whole thing up. Maybe I should do a distance on two of the ordinates independently first. Anyway, I don't know what's up. 

Dear all;

 

Hello everybody, I need your help to dispaly some values obtained using my function f. When I run the code there is no results obtained. Many thanks.

restart:

# The vectors e(i) satify the folowing conditions
e(0)*e(1)=e(n-1) assuming  1<n;
e(0)*e(0)=e(2):
e(1)*e(1)=e(n-1) assuming  1<n: :
e(2)*e(1)=e(n) assuming  1<n:
#
for i from 1  to n-1 do
e(i)*e(0)=e(i+1);
end do:

# We define the function f
f:=e(0)->e(0)+(n-3)*e(1);
f:=e(1)->(n-2)*e(1);
for i from 2  to 3 do
f:=e(i)->(n+i-3)*e(i)+(i-1)*(n-3)*e(n-3+i);
end do:

for i from 4 to n do
f:=e(i)->(n+i-3)*e(i)
end do:

# We define the two vectors
x:=sum(alpha(k)*e(k),k=0..n);
y:=sum(beta(k)*e(k),k=0..n);

#Question : I would like to compute the following  but there is no display of the solution. 
(x*y);
f(x*y);
f(x);
f(y);
x*f(y);
f(x)*y;
f(x*y)- f(x)*y-x*f(y);

Hello all,

I am trying to create a matrix throgh a procedure but unfortunately I could not find how to define an empty matrix in Maple. This is how my procedure looks like,

Initial_Matrix := proc (BC::list, n, BP::list)
local M::Matrix, i;
M:=[][];
for i from 1 to numelems(BC) do
M:= <M,function_coeffs(BC[i], n, BP)>;
end do;

where the procedure function_coeffs returns a row vector. Logically, it should keep adding rows function_coeffs(BC[i], n, BP) to the emty matrix, but its not happening. Please help me out.

Thank you for your time.

Hi,

what i have to calculate in Maple for q=Mises-Stress and the third invariant of deviatoric stress as r???

It is for:

Link:
http://www.egr.msu.edu/software/abaqus/Documentation/docs/v6.7/books/stm/default.htm?startat=ch04s04ath116.html

 

Suppose we are given a matrix M[n*2n] of n linearly independent row vectors. Then I am trying to create a procedure to add n more linearly independent row vectors to this matix resulting in to a matrix M[2n*2n].

Consider this easy example, if the given matrix $M_{2*4}$ is

Matrix(2, 4, {(1, 1) = 4, (1, 2) = 1, (1, 3) = 0, (1, 4) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 4, (2, 4) = 1}) 
then we can add $[0,0,0,1]$ and $[0,0,1,0]$ to obtain the matrix $M_{4*4}$

Matrix(4, 4, {(1, 1) = 4, (1, 2) = 1, (1, 3) = 0, (1, 4) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 4, (2, 4) = 1, (3, 1) = 1, (3, 2) = 0, (3, 3) = 0, (3, 4) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 1, (4, 4) = 0})

 

Can I use rref form?Thank you for your help.

 

I am tryng to write and solve within Maple the equations of movement for a single body subject to a gravitatinal field, F_=-GMm/r2

But I get an error message when trying to define the Angular Momentum which doesn't make sense for me.

Thank you for any help on this topíc.


restart; with(Physics[Vectors]); conventions, Setup(mathematicalnotation = true);
conventions, [mathematicalnotation = true]
r_ := rho*_rho;
r_ := rho _rho
_rho(t);
_rho(t)
rho(t);
rho(t)
v_ := diff(rho(t)*_rho(t), t);
/ d \ / d \
v_ := |--- rho(t)| _rho(t) + rho(t) |--- phi(t)| _phi(t)
\ dt / \ dt /
a_ := diff(%, t);
/ 2 \
| d | / d \ / d \
a_ := |---- rho(t)| _rho(t) + 2 |--- rho(t)| |--- phi(t)| _phi(t)
| 2 | \ dt / \ dt /
\ dt /

/ 2 \ 2
| d | / d \
+ rho(t) |---- phi(t)| _phi(t) - rho(t) |--- phi(t)| _rho(t)
| 2 | \ dt /
\ dt /

eq[1] := -G*M*_rho(t)/r^2-a_ = 0;
/ / 2 \ 2\
| G M | d | / d \ |
eq[1] := _rho(t) |- --- - |---- rho(t)| + rho(t) |--- phi(t)| |
| 2 | 2 | \ dt / |
\ r \ dt / /

/ / 2 \
| / d \ / d \ | d |
+ _phi(t) |-2 |--- rho(t)| |--- phi(t)| - rho(t) |---- phi(t)|
| \ dt / \ dt / | 2 |
\ \ dt /

\
|
| = 0
|
/
Eq[1, 2] := seq(Component(lhs(eq[1]), n) = 0, n = 1 .. 2);
/ 2 \ 2
G M | d | / d \
Eq[1, 2] := - --- - |---- rho(t)| + rho(t) |--- phi(t)| = 0,
2 | 2 | \ dt /
r \ dt /

/ 2 \
/ d \ / d \ | d |
-2 |--- rho(t)| |--- phi(t)| - rho(t) |---- phi(t)| = 0
\ dt / \ dt / | 2 |
\ dt /
NULL;

L_ := `&x`(r_, m*v_);
Error, (in Physics:-Vectors:-&x) found the unit vector _rho present also as a function _rho(t); either one form or the other - not both - can be present in an algebraic expression

 

If I have an N-dimensional vector V and a polytope defined as the set of solutions to the equation A*x = b, where A is a d X N matrix, and b is a d X 1 vector, how can I project V onto the surface defined as above?

Thanks!

Hello,

I'm having an issue with this and I can't seem to fix it. Any help is greatly appreciated! Thank you in advance!!! For some reason mapleprimes won't let me upload the worksheet so I have pasted it below this message.

 

Kind regards.

Gambia Man

with(LinearAlgebra); UseHardwareFloats; with(plots); interface(rtablesize = infinity); with(Statistics);
L := 4; U := 1;
V := proc (x, y) options operator, arrow; piecewise((1/4)*L <= x and x <= (1/2)*L and (1/4)*L <= y and y <= (1/2)*L, U) end proc;
plot3d(V(x, y), x = 0 .. L, y = 0 .. L);

Vij := proc (ni, mi, nj, mj) local Xi, Xj; option remember; global U, L; Xi := 2*sin(ni*evalf(Pi)*x/L)*sin(mi*evalf(Pi)*y/L)/L; Xj := 2*sin(nj*evalf(Pi)*x/L)*sin(mj*evalf(Pi)*y/L)/L; return U*(int(int(Xi*Xj, x = (1/4)*L .. (1/2)*L), y = (1/4)*L .. (1/2)*L)) end proc;
HamilMat := proc (K::integer) local ni, mi, nj, mj, N, Hamil, Eigenvec, i, j, res; option remember; global Vij, U, L; N := K^2; ni := Vector(N); mi := Vector(N); nj := Vector[row](N); mj := Vector[row](N); for i to N do for j to K do res := (i+K-j)/K; if type(res, integer) = true then ni[i] := j; nj[i] := j; mi[i] := res; mj[i] := res end if end do end do; Hamil := Matrix(N, shape = symmetric); for i to N do for j from i to N do if i <> j then Hamil(i, j) := Vij(ni[i], mi[i], nj[j], mj[j]) elif i = j then Hamil(i, j) := Vij(ni[i], mi[i], nj[j], mj[j])+(1/2)*(ni[i]^2+mi[i]^2)*Pi^2/L^2 end if end do end do; Eigenvec := Eigenvectors(Hamil, output = ['values', 'vectors']), Hamil end proc;
SigFigEi := proc (Location::integer, SigFig::integer, VecSize::integer) local values, Eig, i; global HamilMat, OptK; Eig := Vector(VecSize); for i from 2 to VecSize do Eig[i] := HamilMat(i)[1][Location]; if evalf(Eig(i), SigFig+1) = evalf(Eig(i-1), SigFig+1) then OptK := i; break end if end do; values := evalf(Eig[i], SigFig); return values, OptK end proc;
BasisFunc := proc (location) local ni, mi, func, i, j, p, N, res, BasisSol; global HamilMat, L, OptK; p := evalf(Pi); N := OptK^2; ni := Vector(N); mi := Vector(N); for i to N do for j to OptK do res := (i+OptK-j)/OptK; if type(res, integer) = true then ni[i] := j; mi[i] := res end if end do end do; func := Vector(N); for i to N do func[i] := HamilMat(OptK)[2][i][location]*sin(ni[i]*p*x/L)*sin(mi[i]*p*y/L) end do; BasisSol := unapply(add(func[i], i = 1 .. N), x, y); return plot3d(BasisSol(x, y), x = 0 .. L, y = 0 .. L), func end proc;
BasisFunc(1);
Error, (in Vector) dimension parameter is required for this form of initializer

Hello all, first time posting here.  Anyways, I am trying to make a velocity vector plot from known u and v velocity components at known x and y locations i.e. at (0.1, 0.1), u=-0.5 m/s and v=0.43 m/s.  My calculations have been performed with a FORTRAN code, and the results output in matrix format.  What I am left with are 4 vectors: x location, y location, u velocity component, and v velocity component.  I know that MATLAB has a specific function, called "quiver," for creating a velocity vector plot from known data points (through the use of vectors containing those points).  I have trouble believing that there is not an equivalent function, however convoluted it may be, in Maple to create these plots.  I have tried to manipulate the fieldplot, vectorplot, and arrow functions to no avail.

Any ideas on how to create a velocity vector plot from known data points?  I am using Maple 18.

Here is a link to an example of what I am looking for:

https://people.sc.fsu.edu/~jburkardt/py_src/spiral_data/spiral_matplotlib.png

Thanks, Brandon.

Hello, I am working on an elektromagnetism assignment, and have hit a slight problem with the Physics[Vectors] package. The cartesian unit vectors are set to i, j and k. but at my school we have always used x, y and z for both the cartesian unit vectors and the cartesian coordinates. Is there a way to change this so that i am able to use _x, _y and _z istead?

 

P.S. This is my first post so sorry if it is badly done.

Hi everyone,

I wanted to know if it was possible to store symbolic variables in a vector because I have a function of 10 variables which I want to write with Σ, Π and A(l), A being my vector of symbolic variables...

Thank you for your time


Benjamin

Hello!

Vectors are defined by a lenth and an angle.

The arrow command wants me to give it a starting point and an endpoint in coordinates, for example:

with(LinearAlgebra);

with(plots);

gu := arrow([5, 10], [15, 20], .2, .4, .1, color = green)

I'm dealing with vectors from electrical circuits, so the information I have would be the lenth of the vector, and an angle given in degrees (which I could convert to radians if needed for the plot).

The goal is to plot all the vectors in the circuit on the same display so that it can be seen how they are in relation to each other, and to visualise the solution. This means I would need to find out the endpoint of the arrows to start the next arrow at that location with its angle in reference to the same axis. Does it make sence? I'm not sure about the wording of this.

Here's a picture of a calculated vector lenth and angle: http://imgur.com/hJtNeGg

So, the -45 degrees would need to be in reference to the y-axis for example, then the next vector can be placed at lets say perpendicular in reference to the vector sol_1, starting at its endpoint.

The only other way I have to draw this, excluding doing it by hand, is using the object tools from microsoft word, and to be honset I dislike using word for math stuff. I would much rather learn the syntax for Maple to do this.

 

I hope someone knows what I mean, Please ask if somethig is unclear,

Thanks!

 

 

 

 

 

 

I want to make a dot product like

Eq6 := (sum(delta[i].v[i], i = 1 .. 3)).(sum(delta[j].w[j], j = 1 .. 3))

and get the complete result like

Eq7 := v[1].w[1]+v[2].w[2]+v[3].w[3]

but instead I get

Eq6 := ((Vector(3, {(1) = 1, (2) = 0, (3) = 0})).v[1]+(Vector(3, {(1) = 0, (2) = 1, (3) = 0})).v[2]+(Vector(3, {(1) = 0, (2) = 0, (3) = 1})).v[3]).((Vector(3, {(1) = 1, (2) = 0, (3) = 0})).w[1]+(Vector(3, {(1) = 0, (2) = 1, (3) = 0})).w[2]+(Vector(3, {(1) = 0, (2) = 0, (3) = 1})).w[3])

Is there a way to force the dot product operation all the way?

 DotProductCommutativity.mw

I define unit vectors as

delta[1] := Vector(3, {(1) = 1, (2) = 0, (3) = 0})

delta[2] := Vector(3, {(1) = 0, (2) = 1, (3) = 0})

delta[3] := Vector(3, {(1) = 0, (2) = 0, (3) = 1})

Then I manually take the dot product as

Eq3 := delta[1].delta[1]+delta[1].delta[2]+delta[1].delta[3]+delta[2].delta[1]+delta[2].delta[2]+delta[2].delta[3]+delta[3].delta[1]+delta[3].delta[2]+delta[3].delta[3]

and get

Eq3 := 3

But when I do this

Eq4 := sum(sum(delta[i].delta[j], j = 1 .. 3), i = 1 .. 3)

I get

Eq4 := Vector(3, {(1) = Vector(3, {(1) = 1, (2) = 1, (3) = 1}), (2) = Vector(3, {(1) = 1, (2) = 1, (3) = 1}), (3) = Vector(3, {(1) = 1, (2) = 1, (3) = 1})})

How can I make Maple give me the first result using the sum command?

 

Here the potential of maple 2015 to the quantitative study of the decomposition of a vector table is shown in two dimensions. Application for the exclusive use of engineering students, which was implemented with embedded components.

Atte.

Lenin Araujo Castillo

Archivo Corregido:  Decomposición_Vectorial_Corregido.mw

1 2 3 4 5 6 7 Last Page 1 of 15