Questions related to the application of Mathematica to problems in physics.

learn more… | top users | synonyms

1
vote
1answer
27 views

Implicit region is not “valid” in ParamatricNDSolveValue function?

I am trying to solve a geometric problem with relation to my Schrodinger equation and its boundaries. Here is my code: ...
0
votes
0answers
79 views

Plotting: Magnetic field lines

I have data which I am trying to plot as 3D field lines. I need to plot points at coordinates i, j and ...
1
vote
0answers
30 views

How to evaluate Poisson Bracket? [closed]

I'am new in Mathematica so please don't yell at me. I want to define Poisson Bracket in Mathematica 9 and evaluate it for two functions,but it didn't worked. Could you help me please? ...
0
votes
1answer
63 views

How to plot and visualize a single linear vector in 3D? [closed]

I have a vector (in physic) designated as F1=250cos(60)i+250cos(60)j+250cos(45)k, and i would like to see it in a 3D graphic with the axis centered at the origin, after what i would include other ...
1
vote
1answer
55 views

Particle motion with discrete values for the potential

Sorry if the question seems trivial, but I am new to physical simulations in Mathematica. I have a numerical estimation of an electric potential $V(x)$ in 3D, to be more precise I have a matrix ...
0
votes
0answers
61 views

How to filter my signal?

As you can see, my signal is superimposed by an other sinus which has a slightly longer frequency. I would like to filter this out. I think I can do the filtering process, but there is one step that I ...
1
vote
1answer
185 views

NSolve fails to solve a nonlinear equation

Mathematica can not solve this: g = 9.82; ω = 0.5; h = 5; y0 = 1; v = 0; τ = 0; NSolve[h + v t - (g t^2)/2 == y0 Sin[ω t], t]; The error code is: ...
0
votes
0answers
32 views

Failure to evaluate functions of derivatives of Theta functions

I define a variant of the EllipticTheta function (multiplied by a constant factor) : ...
1
vote
1answer
41 views

Variables within functions within DSolve

Hi I am just beginning to learn Mathematica and this is my first time I have been exposed to any type of coding. I am encountering a problem in a basic physics problem. For example we are always ...
3
votes
0answers
133 views

Numerical integration: complicated 2D integral seems to be poorly estimated

In the course of some physics research I've been working on, a very annoying integral has appeared that I'm having difficulty evaluating numerically. Any help you could offer would be greatly ...
0
votes
0answers
26 views

Can't Load .rst file

I've been working with the FeynRules 2.1 package in Mathematica 10. If I download the MC4BSM_2012 FeynRules files and run the Mathematica Notebook I've got no problem in Ubuntu (A lot of warnings ...
2
votes
0answers
132 views

How to control pedestrians to enter the simulation region with random functions?

I am working on simulating crowds of a cross-typed region, given that pedestrians entering from 4 gates of north, east, west and south, ...
0
votes
1answer
74 views

How can I use manipulate on multiple lists within a function?

So I have my N-Body Problem set up like so, and I can plot it in ListPointPlot3D, but I cannot figure out how to stick this into a manipulate function. ...
4
votes
2answers
193 views

Embedded Graphics: General Relativity

Is there a way to make the following graph with two dimples? Like the central indent and another indent to the side to simulate a central mass like the sun and another smaller mass like the Earth in ...
2
votes
0answers
85 views

Position from 3 dimensional acceleration data? [closed]

Per a previous post on Mathematica StackExchange here I'm trying to utilize 3-axis accelerometer data to determine approximate 3D position. I wanted to start with a simple case. To this end I've ...
3
votes
1answer
145 views

Introduction to Vectors and NDSolve

I want to use NDSolve for a complicated equation that contains vectors. I am still new to Mathematica, and have really enjoyed what I have learned thus far. Concerning this question, I have already ...
5
votes
3answers
278 views

Numerically Simulating a Simple Differential Equation: RC Circuit

I need to learn to simulate a differential equation in Mathematica. The differential equation that I need to simulate is complicated, without an analytical solution. However, I am brand new to ...
5
votes
2answers
158 views

Multidimensional integral to compute the end-to-end distribution of a FENE ideal chain

How do I evaluate the mean end-to-end squared distance of a FENE ideal chain at fixed inverse temperature $\beta$ in the canonical ensemble? This quantity is defined as the mean value ...
1
vote
0answers
71 views

Lens profile in a grid

So I have a FDTD code which works and its all good the propagator of the code it self is: ...
0
votes
0answers
113 views

Solving for a Transfer Matrix

I am currently working on a program that will solve a transfer matrix and extract several elements. Particularly, the transfer matrix is characterized by four elements in a ...
1
vote
0answers
75 views

Maxwell propagator problem 2D [closed]

Hi folks I`m trying to turn certain matlab code to mathematica but for some reason it doesnt really want to do it could you please assist me a bit. Thanks The matlab code is as follows: ...
0
votes
1answer
78 views

Plotting and Analysing derivatives of solution of differential equation [closed]

Using NDSolve I found the solution for the differential equation that governs the dynamics of a simple pendulum. How do I plot the derivatives of the solution? Also, I would like to find the ...
1
vote
0answers
69 views

Creating an orrery

I am interested in creating an orrery, ideally with SystemModeler. I have been experimenting with the MultiBody library since that allows you to do animations. The most relevant example seems to be ...
2
votes
2answers
456 views

draw part of parabola using graphics

today I was doing my physics homework. one of the questions was draw velocity-time and position-time graphs from the acceleration-time graph. So for the first step I draw acceleration-time using ...
0
votes
0answers
49 views

SetDelayed::write is Protected error on Linux not PC

I have 2 computers, a PC and an Ubuntu computer both with Mathematica 10. I am running the exact file with the exact support files on both computers, but just on the Ubuntu I get the following error ...
1
vote
0answers
117 views

Finding the point of best focus for a ball lens

This is my project (simulation): A set of parallel rays goes trough a ball lens (circular lens). Because of spherical aberration, refracted rays will not meet in a single point. Find the distance ...
5
votes
2answers
165 views

Pull out scalars from NonCommutativeMultiply in commutator of SU2 spin algebra

I am trying to work on spin chains in Mathematica. That is, N spin-1/2 on a one-dimensional lattice with periodic boundary condition (N+1=1). I need therefore the SU2 spin algebra $S_z, S_+, S_-$ with ...
10
votes
2answers
459 views

How to calculate position from 3 dimensional acceleration data?

I discussed this problem yesterday with Pickett here. We concluded that the best solution may be just Interpolate and then ...
2
votes
1answer
138 views

2D Fourier transform of a few (4) disjoint discs on a plane

I'd really appreciate some advice. Short Version I'm trying to calculate the following $$ ...
8
votes
1answer
291 views

Calculate inertia tensors

Moments of inertia and inertia tensors are well-known characteristics of rigid bodies in physics and applied mathematics. How to calculate them in Mathematica? Is possible to use regions which was ...
5
votes
1answer
276 views

NDSolve and Piecewise functions

I want to clarify upfront, This is not a please help me get the answer to my homework question. I know what I need to do I just do not know how to get mathematica to do it. I am working on an ...
7
votes
1answer
285 views

How to draw quantum circuit with mathematica?

How to draw quantum circuit as the following one with mathematica? I only find the quantum package which is designed for MMA 7.0 as the author said. But the version of my Mathematica is 9.0 and when ...
-1
votes
1answer
326 views

Navier–Stokes equations [closed]

how write Stokes equation in mathematica? the Stokes equation is in in this link http://en.wikipedia.org/wiki/Navier-Stokes_Equations
7
votes
3answers
163 views

Substituting values in 3j symbols

I believe these three expressions should give the same answer of -1/Sqrt[3], but the first gives +1/Sqrt[3] whereas the second ...
3
votes
1answer
123 views

Interpolation function rejected by ParametricPlot3D

I was playing around with the bouncing ball example, and decided to make a 3d version. I decided to use ParametricPlot3D for a simple independent 3D ...
1
vote
0answers
114 views

Solving PDEs with complicated boundary conditions

The system I'm trying to solve is $$\nabla^2 C_{(r,\theta)} =0$$ $$C_{(\infty,\theta)}=C_0$$ $$ [ \frac{\partial C_{(r,\theta)}}{\partial r} \cos(\theta)+\frac{1}{r}\frac{\partial ...
1
vote
1answer
57 views

How to make a discretized NMinimize more precise

I am using Mathematica for physics research and I want to minimize a Hamiltonian equation with respect to two variables (I have also discretized the problem). I have a single constraint. When I plot ...
0
votes
1answer
102 views

Help with Computation-heavy triple integration

I'm Trying to use Mathematica to symbolically integrate a function with respect to x, y, and z. For some reason, this computation either takes about, 60s and returns the input, or takes many hours ...
0
votes
2answers
78 views

Numerical simulation - getting slower the more data I have [closed]

I'm simulating the trajectory of a charged particle. I have the coordinates saved in the list $coor$ ({{0,0,0},{0.1,0.1,0.003}, etc}) as the time goes on this array is getting bigger and the ...
1
vote
0answers
88 views

Help wanted in evaluating an integral

Trying to evaluate an integral. To be specific, I'm trying to find the expectation value of $x^2$ for the ground state of hydrogen. Here is the code I used: ...
0
votes
1answer
211 views

Solving for the time-evolution operator in a periodically driven system

I am looking at the Hamiltonian $$H(t)=\begin{pmatrix} 0 & e^{i\Omega t}\\ e^{-i\Omega t}& 0\end{pmatrix}$$ I am trying to solve for the unitary operator $U(t,0)=\mathcal{T}\exp(-i\int_0^t ...
15
votes
1answer
1k views

Animation of double pendulum

Sadly, I am a completely newbie. I am studying Physics and in our theoretical physics class we got the task to solve the double pendulum using Mathematica. We just got the program, but no introduction ...
3
votes
1answer
124 views

How to simulate the 1-D motion of a particle under constant acceleration? [duplicate]

I am trying to represent the position of an object under freefall, but am completely lost. Here is my attempt so far: ...
1
vote
1answer
306 views

Traces of products of Pauli matrices

Given a matrix $M$ of shape $2^L*2^L$, I would like to compute all the traces $\text{Tr}( M.(\sigma^{n_1}\otimes\sigma^{n_2}\otimes\ldots\otimes\sigma^{n_L})) $ for $n_1=1...4$, $n_2=1...4$, ..., ...
1
vote
0answers
72 views

Solving for a charge distribution from a dipole moment [closed]

I've been looking around but I haven't been able to find any way of solving this conceptually simple problem. In short I am trying to find the charge density in a dielectric in an external electric ...
1
vote
2answers
146 views

Making mathematica do regulated integrals

Consider the integration, $\int _0 ^\infty dx\ x \tanh( \pi x) \sqrt {x^2 + a^2 } $ where $a$ is a real number. This integral is divergent. We note that an asymptotic expansion of the integrand ...
1
vote
1answer
111 views

Operator operation

There are two basic operators r={x,y,z} and p=-I*h*{D[#,x] , D[#,y], D[#,z]} &. The angular momentum operator is defined as ...
8
votes
3answers
475 views

1D quasicrystal: points on a line nearest to points on a lattice

I have a simple lattice / line manipulation: ...
1
vote
2answers
152 views

How to fit exponential function which is not converging

I performed an experiment where I measured diffusion in a lipid bilayer. I performed the measurement by fluorescently labeling some molecules in the bilayer, and the irradiating them with light ...
0
votes
0answers
198 views

Boundary Value Problem- using NDSolve or another method

I am trying to solve a set of coupled partial differential equations, with defined boundary conditions using mathematica. Here are the equations and the boundary conditions. ...