Tagged Questions

Questions on dealing with vector calculus functions of _Mathematica_ such as `Grad`, `Div`, `Curl`, `Laplacian` and their rperesentations in various coordinate systems.

learn more… | top users | synonyms

0
votes
0answers
56 views

Using NMinimize to solve the Multi Commodity Flow problem

I have the following optimization problem and I was hoping to solve it with NMinimize. Unfortunately, I don't get the desired result. I would like to get an optimal vector for $d_{e}$. This is how ...
2
votes
1answer
100 views

How to simplify an expression using a differential equation as an assumption?

I have a fluid dynamics flow described by a 3D velocity field such as u = {ux[x, y, z], uy[x, y, z], uz[x, y, z]} I want to calculate ...
2
votes
3answers
165 views

Hypergeometric function with a matrix argument

I am looking for the evaluation of a Hypergeometric function with a matrix argument as for example in Ploev and Edelman or as showcased in this Wikipedia article. From what I understand from ...
3
votes
1answer
110 views

Bad numerical approximations

I'm trying to do some calculations here, but for some reason Mathematica starts using numerical approximations that are no good for my work. Specifically: ...
1
vote
2answers
36 views

how to extract a matrix and a vector from a vector [duplicate]

I have a vector like this: V = {a*exp1+b*exp2+c*exp3, a*exp4+b*exp5+c*exp6, a*exp7+b*exp8+c*exp9} and I want to get : ...
-1
votes
1answer
57 views

Plotting vector field [closed]

I am trying to plot the gradient vector field of the function: $$sin\sqrt{x^2+y^2}$$ This is my code: ...
-1
votes
1answer
120 views

Finding the maximum of a gradient vector

Im trying to find the maximum of my gradient vector G[x,y], I've tried several options including FindMaximumValue, FindMaximum etc. but i couldn't find it. The full function is shown below, any help ...
-2
votes
1answer
164 views

How to compute the curl of a vector field?

Can anyone help me out with this? I am trying to get Mathematica to evaluate curl of a vector field. The vector field is ...
3
votes
1answer
98 views

Working with abstract vectors

I often need to compute derivatives or integrals involving N-dimensional vectors (where the dimension could be equal to 2 or 3 but is not particularly relevant for the sake of the derivation). The ...
3
votes
1answer
256 views

Plot a function based on Derivative - Gradient field

I'm trying to plot the gradient field of a function in such a way that is possible to change it easily, only editing the function. Consider the code: ...
0
votes
0answers
30 views

Is it possible to obtain vector differentiation results in terms of vectors rather than components? [duplicate]

When doing certain operations on vectors in Mathematica it is quite desirable to convert the results back into coordinate-free notation. Unfortunately, after going through the Documentation Center and ...
0
votes
0answers
70 views

How can I do symbolic vector differentiation in mathematica? [duplicate]

I want to do this calculation using Mathematica and I do not know how. I want to finf gradient of a vector function I Have this $$ f:R^n\to R $$ $$ \nabla_xf(x)=\nabla_x(x^tBx) $$ I want to get ...
0
votes
0answers
31 views

How to keep the vector product form?

When we deal with the vector product problem, we usually should define some fixed form for the vector product. For instance, I define it in this way, ...
8
votes
1answer
667 views

How do I plot a proper streamline plot, including spacings and line endings?

Mathematica includes two nice built-in tools to visualize vector fields, VectorPlot and StreamPlot. The latter is a useful tool, ...
0
votes
1answer
212 views

Symbolic matrix simplification and solving

I am trying to expand or solve an equation that contains matrices for a certain matrix, but it is not working: ...
1
vote
1answer
128 views

Rotating a set of 3D points against an arbitrary axis

I want to rotate a list of 3D points around unit vector (1,0,1) by the angle of Pi. How can I accomplish it?
1
vote
1answer
86 views

Simplify ignores DiracDelta contribution from Coulomb sources

Executing the following line, Simplify[Div[Grad[-1/Sqrt[{x, y, z}.{x, y, z}], {x, y, z}], {x, y, z}]] yields 0. However, this is incorrect, since ...
0
votes
1answer
112 views

replace is not happening for specific vector valued expression

Given x[t_] := {Subscript[x, 1][t], Subscript[x, 2][t]} following replace works fine x[t] /. x[t] -> 0 0 but following ...
3
votes
1answer
106 views

Integration over a convex combination of a region: $\int_{\Omega} (w_1 z_1 + w_2 z_2)^{1-\sigma} d (z_1, z_2)$ where $\Omega = \{ z_1 + z_2 = 1\}$

Take $w,z\in R^{n}$. I am interested in integrating (as generically if possible) $$\int_{\Omega}(w \cdot z)^{1-\sigma} d z$$ Where the domain of $\Omega$ is $1$ dimension, and includes the convex ...
1
vote
0answers
97 views

Simple Jacobians, Gradients, etc. with arbitrary length vectors/matrices?

Is there any way (or a package built for it) can do simple operations with vectors and matricies of arbitrary size, but conforming extents? For the simplest example to test, given an arbitrary vector ...
1
vote
1answer
161 views

Blank VectorPlot

I'm trying to draw a VectorPlot of a function that admittedly probably is not very well suited for such a plot. But at the moment I just get a white sheet and I ...
2
votes
0answers
50 views

Draw Vector Tangent to Point on Parameterized Curve [duplicate]

I define the vector function $\vec{r}(t)$ by r[t_] = {t, Cos[t], Sin[t]} The command r'[t] gives the vector tangent to the ...
3
votes
1answer
596 views

Normal and Tangent of Acceleration in 3D

I need to figure out how to find the Normal and Tangent of acceleration. I know the formula for the tangent of acceleration is $((Acceleration . Velocity)/(Velocity.Velocity))*Velocity$ and the normal ...
1
vote
1answer
129 views

Difference between vector's Dimensions

What's the difference between projection= Table[Conjugate[autoVett[[k]].someVector1], {k, 1, r}] Dimensions[projection] that gave me as results ...
7
votes
2answers
404 views

How to declare a 3D vector variable?

How can I do vector calculations without telling Mathematica the vector entries? I have very many arbitrary linear combinations in $\mathbb{R}^3$ which I want to perform some general calculations on ...
0
votes
1answer
217 views

Projection of a set of 4D points to the 3D space

Consider the following set of points in the 4D space: ...
1
vote
2answers
282 views

Doing vector manipulations on Mathematica

This is hopefully a simpler version of this previous unanswered question of mine. Let me just focus on the two expressions $F_2^{(s)}$ and $F_3^{(s)}$ given in A.3 and A.4 of page 19 of this paper. ...
4
votes
1answer
377 views

Checking if an expression is equal to zero

I have 6 points in $\mathbb R^3$ as follows: ...
2
votes
2answers
418 views

CoordinateTransform from Cartesian to Spherical

I have a set of coordinates on the sphere in terms of cartesian coordinates. I would like to convert them to spherical data points. However, when I try to convert the point {0, 0, 1} it tells me that ...
2
votes
0answers
227 views

Derivatives of list elements

Could someone explain the odd behavior of the Derivative function when drawing arguments from lists? We have, ...
6
votes
1answer
740 views

Smoothing/Averaging 2D Vector Fields

I have a list of 2D vectors defined by {{x,y},{u,v}} and would like to smooth or average the vectors. For example here are 2 vector fields, the second has noise ...
2
votes
2answers
346 views

Why are these flow lines cut short?

I want to plot some flow lines for a certain vector field, but Mathematica doesn't seem to plot the entire flow line, only a small piece of it. How can I fix this? ...
3
votes
1answer
150 views

Extending D for four-derivatives

I'm currently trying to consistently define rules for extending D[] to four-derivatives. As 'backend' I'm using the package TRACER (http://library.wolfram.com/infocenter/MathSource/2987/), which can ...
6
votes
2answers
1k views

Is it possible to do vector calculus in Mathematica?

I am trying to rearrange and manipulate some vector differential equations in Mathematica. As far as I understand you have to tell Mathematica that a variable is a vector by specifying the components ...
15
votes
4answers
6k views

Finding unit tangent, normal, and binormal vectors for a given r(t)

For my Calc III class, I need to find $T(t), N(t)$, and $B(t)$ for $t=1, 2$, and $-1$, given $r(t)=\{t,t^2,t^3\}$. I've got Mathematica, but I've never used it before and I'm not sure how to coerce ...
1
vote
1answer
318 views

How to find this convolution?

How can I use Mathematica to find the convolution $f*f$ for ...
2
votes
3answers
479 views

How do you show a cone inside of a sphere?

I have a sphere of radius 5 and a cone (i.e. z=√(x^2+y^2)) inside of it. I can show them separately, but I'd like to show them together. I've tried a few tricks with opacity, implicit specifications, ...
3
votes
1answer
617 views

Plotting a complicated vector field

If we consider the vector $\left ( A \cdot \nabla \right) \: B$, we have in Cartesian coordinates $$\left ( A \cdot \nabla \right) \: B = \left ( A \cdot \nabla B_x \right ) e_x + \left ( A \cdot ...
2
votes
1answer
872 views

Lie-Bracket of two vector fields

I'm new to Mathematica (installed couple of hours ago) and I need to compute a few Lie brackets between two vector fields $f$ and $g$. $$ f\left(\mathbf{x}\right) = \left( \begin{array}{c} ...
4
votes
2answers
657 views

How to change coordinates of a differential operator?

I'm doing a basic quantum mechanics problem and am trying to learn how to do it in Mathematica. Any help would be much appreciated. $\vec{L} = \vec{x} \times \vec{p}$ where $\vec{x}$ has ...
12
votes
4answers
2k views

How do I plot the unit normal field for a surface?

The question is pretty much in the title; I'm about to teach my multivariable calculus students about orientations on surfaces, and I would like to be able to show them pictures. Any ideas?
6
votes
2answers
983 views

Singularities using VectorPlot

I am trying to plot a vector function of a fluid flow given by $\vec{V} = (\frac{-\cos(\theta)}{r^2},-\frac{\sin(\theta)}{r^2})$ I am trying to plot it in Mathematica using the command below, I ...
2
votes
1answer
348 views

Vector Analysis: Why are the Grad, Laplacian and Div being evaluated to zero

For the non linear partial differential equation below, why are the Gradient, Laplacian and Divergence being evaluated to zero despite using the VectorAnalysis ...
0
votes
1answer
8k views

Why I get the “Set::write: ”Tag Times in is Protected." error?

Why I get the error below with this code: ...
1
vote
1answer
810 views

Calculation of Jacobians for a finite element method simulation

I'm attempting to calculate symbolic expressions for the (Newton's method) Jacobian of a Galerkin finite element method using Mathematica. It would seem that this is a perfect application for symbolic ...
2
votes
1answer
345 views

Laplacian or Grad of an InterpolatingFunction

For obtaining a Schlieren image from an equation for density, I need to calculate the first derivative of density and make a contour plot. The below code snippet calculates the first derivative w.r.t ...
3
votes
2answers
1k views

How to substitute numeric values in a symbolic Jacobian matrix?

I have a multi-variate function from $\mathbb{R}^n\to\mathbb{R}^n$. Choosing any desired initial vector, we can produce the corresponding function value, which is a vector as follows. The main problem ...
16
votes
2answers
7k views

How to make Jacobian automatically in Mathematica

If we have two vectors, a and b, how can I make Jacobian matrix automatically in Mathematica? $$ a=\left( \begin{array}{c} \text{x1}^3+2\text{x2}^2 \\ 3\text{x1}^4+7\text{x2} \end{array} ...
2
votes
2answers
366 views

Finding a vector field from a set of points

I have a data set of the type $\{x,y,z\}$ where $(x,y)$ is a point and $z$ is the "value" or magnitude at that point. This gives me a triangular sort of shape since $(x,y,z)$ is not defined along the ...
10
votes
2answers
547 views

Creating randomly oriented planes

I would like to create randomly oriented planes. This is how I'm attempting to do that: I create a 2 random unit vectors, $\mathbf{v}_1$, and $\mathbf{v}_2$, in the $x$-$y$ plane I assume that if I ...