After completing my sophomore year at college as a Physics major and I realized that the laptop I currently have is not able to keep up with the demand I put it under for school. The laptop I currently have has a <2 hour battery life and has been known to crap out (didn't know how else to...
While reproducing a research paper, I came across the following equation,
∂f/∂t−(H(f)(∂f/∂x)=0
where [H(f)] is hilbert transform of 'f.'
and f=f(x,t) and initial condition is f(x,0)=cos(x) and also has periodic boundary conditions given by
F{H{f(x′,t)}}=i⋅sgn(k)F{f(x,t)},
where F(f(x,t) is...
Wolfram made what I would call a gigantic move: they released the Wolfram Engine for Developers - for zero cost!
The idea of the Engine is primarily that you can call and run Wolfram Language commands from lots of different sources. I am not the kind of person to automate Mathematica, although...
Hello everybody!
I need a little help with FeynCalc. I think the problem is really simple but I can't find how to fix it.
I want to evaluate a trace coming from a Feynman diagram. Since the particles are all massless, I want to impose the condition on the momentum ##p^2=k^2=p'^2=0##.
I've...
1. Homework Statement
2. Homework Equations
3. The Attempt at a Solution
I used the NDSolve function from mathematic but its giving me problems. What is the correct way to enter the equation?
soln = NDSolve[{y''[t] = (-9.8/5)*sin (t), y[0] = 20, y, {x, 0, 12}}]
Dear friends:
It's strange that Mathematica can do the integral of ##\int_0^\infty dx~x~_2F_1(a,b,c,1-x^2)##, however, fails when it's changed to ##\int_0^\infty dx~x~_2F_1(a,b,c,1-x-x^2)##.
Are there any major differences between this two types? Is it possible to do the second kind of integral...
A problem that I have to solve for my Linear Algebra course is the following
We are supposed to use Mathematica.
What I have done is that I first checked that A is symmetric, i.e. that ##A = A^T##. Which is obvious.
Next I computed the eigenvalues for A. The characteristic polynomial is...
I calculated a simple complex function integration using Mathematica 11.
As the screenshot shows, the only difference between the two inputs is that I set a = 1 in the second one.
However, the two outputs differ by a factor of 2.
What's wrong? Is this a bug?
My objective is to make a list of functions and afterwards be able to make operations with those functions.
Hyfield[list_, bits_] := Module[{i, auxList, Hy},
auxList[a_] := List[];
For[i = 1, i <= bits*2, i++,
auxList[a] =
Append[auxList[a]...
Hi,
I am using Mathematica to calculate density of states and current of the Green's function times self energy in most simple form. I am not sure if I am getting current integral over energy implemented correctly. Shouldnt first current plot be a line with a slope? Below is my code...
I have some data points that I want to plot and suppose I have the function ##f(x) = x##, with the domain having the range ##0\leq x \leq 10##. Assume that the experimental values lies in the range ##4\leq x \leq 7##, how can I put a rectangular region to cover this range behind my plot so as to...
Hi there guys,
Currently writing and comparing two separate Mathematica scripts which can be found here and also here. The first one I've slightly modified to suit my needs and the second one is meant to reproduce the same results.
Both scripts are attempting to simulate the trajectory of a...
Hi, I need some helps. I try to run the code below but its keep give me the value of intial value.
Q = 1415; \[Mu] = 1000; \[Sigma] = 500;
double dx;
double x;
double fx;
double fxplusdx;
double I1;
double I2;
I1 = 0;
I2 = 0;
For[x = 0, x <= Q - dx, dx += 0.01,
fx = PDF[NormalDistribution[0...
Hi all! I'm new to Mathematica.
I have written a code for performing a convolution integral (as follows) but it seems to be giving out error messages:
My code is:
a[x_?NumericQ] := PDF[NormalDistribution[40, 2], x]
b[k_?NumericQ, x_?NumericQ] := 0.0026*Sin[1.27*k/x]^2
c[k_?NumericQ...
Suppose I have a differential equation
$$\ddot \phi + 3H (1+Q) \dot \phi + V_{,\phi} = 0$$
where ##\phi## is the inflaton field. ##H## is the Hubble parameter, ##Q## is just a number, ##V_{,\phi}## is the derivative with respect to ##\phi##, and initial conditions given by ##\phi[0] = 2...
I want to plot a function f(x) such that it's domain is from 1 to 10000 but I want the graph to read 1, 10, 100, 1000, 10000 in the x-axis, while the y-axis is normal.
##f(x) = 2\sin^{-1} ( e^{-\frac{60}{x}})##
From cosmology,
##H^2 = \frac{ρ}{3M_p^2} = \frac{1}{3M_p^2}(½\dot φ^2 + ½m^2φ^2)##
Suppose ##V(φ) = ½m^2φ^2##
where
##ρ## = density
##M_p## = planck mass
I want to graph ##H## vs. ##φ## but there is a ##\dot φ## and I know this is a differential equation, can somebody help me what to do here?
I'm very new to Mathematica/programming and I want to do a theoretical calculation using Mathematica,
suppose I have,
##Y=CX,~~~C=constant##
Now, I want to plot Y vs. X but X should run at every point since every point is a solution for Y, how should I do this? Before, I was thinking maybe I...
1. Homework Statement
I am a rookie to the QFT extension in Mathematica called FeynCalc, and tried to use that into solving some quiz. Soon I met a problem upon some condition presented in a problem which declares an relation of two same tensor with different indices results in some value when...
Hello to all. I've been trying to plot some graphs on Mathematica but I have faced some troubles so far.
I am working on this paper by Lattimer on binary systems:
https://arxiv.org/abs/astro-ph/0002203
I am trying to reproduce Figure 11 but something seems to go wrong.
Here's what I have...
1. Homework Statement
Use Mathematica to calculate the Gaussian curvature of the plane Ax+By+Cz=D, in which A, B, C, and D are constants and C≠0.
Use the following data:
2. Homework Equations
3. The Attempt at a Solution
First I found the line curvature. As here:
That code gets...
Mathematica is fully functional for daily astronomy data analyzing works, and it is much more convenient than python, so why can't them use Mathematica instead of python? It's like when you have warp drive you still travel through space using chemically powered traditional rockets, which is...
I need to define a function that integrates a function in some interval and returns its numerical value. I am not allowed to use built in functions for integrating in Mathematica.
This is my code that won't work:
http://pokit.org/get/?93856fe8f2070ba781028f8634b9ac3a.jpg
This is a code from a...
As I found on many websites, they suggest to use Chop[]. I tried that already but it doesn't work.
This is my output.
{{2 (Conjugate[
I.0] I.0 + (Conjugate[I.0] + Conjugate[Subscript[v, 1]]/Sqrt[
2]) (I.0 + Subscript[v, 1]/Sqrt[2])) Subscript[\[Lambda],
1] (Conjugate[I.0]...
I have an integration that doesn't have a solution in the table of integrals. So, I evaluated it using Mathematica using the command NIntegrate. However, when I compare the result with Monte-Carlo simulations, there is a very small constant gap between the two curves. Is it because of the...
So, I'm trying to plot a 3D "dipole" (an arrow with a small torus around it basically) in mathematica, and want to rotate it according to Euler angles...
I use this code for the rotation matrix:
rot[a, b, g] := RotationMatrix[g, {1, 0, 0}].RotationMatrix[b, {0, 1, 0}].RotationMatrix[a, {0, 0...
I recently plotted a piecewise function:
Plot[Piecewise[{{1 - Exp[-.002*t],
0 <= t < 120}, {-Exp[-.002*t] + Exp[-.002*(t - 120)],
120 <= t}}], {t, 0, 5000}, PlotRange -> {0, 0.25}]
I then defined the function which I am calling q[t_] as follows:
q[t_] := Piecewise[{{1 -...