Items tagged with vector vector Tagged Items Feed

How to add values in vector? This is what I am trying in my code below. I know there are more than one values satifying my conditon.

What is the best way constructing such vectors. If there are other options, please let me know. Thanks.

Hi all,

I've been working on this problem for a while and only need to figure out why the function vector is returning zeros every time I call the procedure. The problem statement is in the file I attached to this.

If any of you could help me, I would be a happy guy.

 

File: 1.2.mw

The physics package is pretty awesome, but one thing that would be a big help is some functionality for entering and resolving vectors in angle/magnitude format (50 angle 45 degrees, or what have you.) Is there a way to do this?

 

Thanks

 

Joe

I'd like to plot the differences between terms in a sequence of vectors. Each difference term should start at the end of the last difference term, so that if I was to plot the actual term, the vector would meet at the end of the difference term. The sequence is limited in length and stored as a list.

Ex:

Suppose I have  [<1,1>,<2,2>,<3,3>]. The difference terms would thus be [<1,1>,<1,1>,<1,1>]. The first difference term would be plotted be plotted from <0,0>, the second starting at <1,1>, the third at <2,2>.

I can compute the difference terms, but I am not sure how to make the plot I desire. Is arrow(...) the answer, somehow?

Hello all,

Suppose I have a vector valued function f (dimensions of vectors is 2). I can use fieldplot to show me the function's behaviour over some window.

I'd like to do the same thing, but I want to compact the space down to the unit circle. Basically, given a vector <x,y>, I'd like to "fieldplot" f(<x,y>) but with each arrow f(<x,y>) appearing at <x,y>/(1 + Norm(<x,y>,2)) instead of <x,y>.

Is this possible?

I've got a vector x=[0.36,1.3279,1.6882] (1*3) obtained as x:=pseudoInverse(A)*b where A=[<2,4,1>|<1,-1,1>|<3,3,2>] where <2,4,1>,<1,-1,1>,<3,3,2> are the columns of A and b is [8,5,4](1*3). Now, when I find the rowspace of A using RowSpace(A) I get the row vectors <1,0,1> and <0,1,1>, neither of which are equivalent to x. How do I arrive at my result that x is in the rowspace?

In Maple15, in spherical coordinates, what statements plot covariant vector [1,r^2*sin(theta),sin^2(theta)] and its contravariant equivalent [1,sin(theta),1/r^2] within their respective bases vectors. Assume r = 2, theta = pi/6 and phi = pi/4. I presume each plot will display the same space vector.

Hi there,

I would like to compute and display the nullclines of a set of ordinary differential equations.

AFAIK, I can compute the nullclines in Maple by defining the equations and solving the system

e.g.:

# Define the equations
eq1 := u(t)*(1-u(t)/kappa)-u(t)*v(t) = 0;
eq2 := g*(u(t)-1)*v(t) = 0;

# Solve the system (i.e. compute the nullclines)
sol := solve({eq1, eq2}, {u(t), v(t)});

However, I am not quite able to imagine how to display them over a dfieldplot or a phaseportrait.

Attached is an example with some differential equations, and their vector field and trajectories: MaplePrimes_Predator_prey_model_nullclines.mw.

It can be use to illustrate how to (compute and) display the nullclines.

 

Thank you,

jon

I experienced strange operation of "union" for sets of vectors.

Mt1:=Matrix(2, 4, [[ 0,1, 0, 0], [ 0,  0,  1, 1]]); Ms := Vector[column](4, [8,4,2,1]); St1 := {}:

St1:= `union`(St1, {Mt1 . Ms});

I am surprised, because each execution of union adds new and the same vector <4 | 3> to set St1:

1

But after copying any set in the clipboard and pasting the set St1 has only one instance of vector <4 | 3>:

2

What does it mean?

Hi all,

 

I'm trying to initialize a 4-vector and promote it as a tensor in Maple 16 using the Physics package. My attempt so far has been unsuccessful :) 

This 4-vector needs to be a function of the coordinates X. What I'm trying to do is to allocate terms by terms my tensor as indicated in my attempt below. Everything looks fine (the query is telling me F2 is a tensor) but when I'm trying to compute contravariant quantities, I can't obtain an answer for a given coordinate (here F2[~n](X) is not equal du v_1(X))... Finally, the last line, computing the sum over all the indexes is not giving me an answer... 

 

restart; with(Physics):
Setup(mathematicalnotation=true):
Setup(coordinatesystems = cartesian):
Define();
g_[]:

F2[1]:=(x,y,z,t)->v_1(X);F2[2]:=(x,y,z,t)->v_2(X);F2[3]:=(x,y,z,t)->v_3(X);F2[4]:=(x,y,z,t)->1;

Define(F2[mu](X));

Define(F2,query);

F2[1](X);F2[2](X);

F2[~1](X);F2[~2](X);

F2[~nu].F2[nu]

 

Surely, something is wrong in my way of defining a spacetime tensor. What would be the solution to obtain a decent contravariant F2?

 

Another more or less related question would be to know how to compute a taylor expansion of a function itself i.e. 1/(1+f(X))~ 1-f(X) since f<<1? Would there be a built-in fonction to do such a thing?

 

Many thanks!

 

Vincent

Dear all,

I woul line integrate e gradient moltiply to a vector. I will try to explain better.

I have define this function:

phi(xi,eta)=(2*xi-1)*(xi-1)*(2*eta-1)*(eta-1)

and I have a vector a

v=(0,1)

I would like to applu a dot product between the gradient of phi an the vector and integrate the results. 

I have already try in many way but without succeed. Someone could please help me?

Thanks 

I have a nonlinear system with 4 equations and 4 unknowns. I am using fsolve. I know that there are multiple solutions for each variable but am only getting one. I need the others. what do I do??

This is my code:

R__1 := Matrix([[1, 0] , [0, 1] ]);

R__2 := Matrix([[1/2, sqrt(3)/2] , [-sqrt(3)/2, 1/2] ]);

R__3 := Matrix([[-1/2, sqrt(3)/2] , [-sqrt(3)/2, -1/2] ]);

R__4 := Matrix([[-1, 0] , [0, -1] ]);

R__5 := Matrix([[-1/2, -sqrt(3)/2] , [sqrt(3)/2, -1/2] ]);

 

d__1 := Vector( [ 0, 5.4] );

d__2 := Vector( [ 6.4, 4.539] );

d__3 := Vector( [ 11, 4.078] );

d__4 := Vector( [ 15.5, 2.079] );

d__5 := Vector( [ 19, 1.039] );

 

a := Vector( [ a__x, a__y] );

 

A__1:=R__1.a+d__1;

A__2:=R__2.a+d__2;

A__3:=R__3.a+d__3;

A__4:=R__4.a+d__4;

A__5:=R__5.a+d__5;

 

OO:=Vector([O__x,O__y]);

 

DA1:=A__2.A__2-A__1.A__1-2*(A__2-A__1).OO;

DA2:=A__3.A__3-A__1.A__1-2*(A__3-A__1).OO;

DA3:=A__4.A__4-A__1.A__1-2*(A__4-A__1).OO;

DA4:=A__5.A__5-A__1.A__1-2*(A__5-A__1).OO;

 

fsolve({DA1,DA2,DA3,DA4},{a__x,a__y,O__x,O__y});

Thanks for any tips you may be able to offer

 

Hi,

I have coded up a vector that is of my interest. The code runs witout any problem and gives me exactly what I want.

newtest.mw

All I want to know, is that if there are more efficient way to do so?

Any tricks, or better use of a particular function that I wasn't aware?

The only tiny bit of unsatisfactory is that, the (1-w) term is at the first term of the addition, is that possible to move it to the last term? Which is more conventional to read.

 

Thanks,

 

casper

Hi all,

Assume that we have a vector, namely v:=[1,-2,3,-4] and we want to construct special matrix namely Z, from Vector v as follow:

first row is 1, secnond row is -2,..., the end row is -4 namely
Matrix([[1,1,1,1],[-2,-2,-2,-2],[3,3,3,3],[-4,-4,-4,-4]]);

in the other word every row of matrix is 4 times corresponding componet in vector.(for example v is (1*4))

how can we do this?

best wishes

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

Hi,

How are you. Nice to meet you once again.

I need your help  to give an answer to  this question:

I have two vector v1:=[1,2,3,4,5,6,7,8]; and v2=[9,10,11,12,13,14,15,16];

I would like to construct a matrix whone first colunm is v1 and last colunm is v2 and all the rest is ZERO.

Many thinks.

 

 

 

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