Mathematica Definition and 1000 Threads

Hello PF. I wanted to optimize time expenditure on some calculation that I'm doing. A natural way to do this can be by computing parts in parallel. Now when I say "parts" I don't mean that some computations are wholly independent of other computations. My situation is that I have a bunch of...

Hi everyone, Iam just wondering how to export the mathematica nb. file to MS word? Preferably word 2003 version. I have tried to use cut and paste method into word, however, the presentation in word does not look nice. Any advice is very much appreciated.

Hi, I've got an interpolating function which has been generated from using NDSolve and I'm trying to find all the values of x for which the y value is equal to 2. I've constructed a (much) easier example to show what I mean. Suppose I have a set of points which I have generated an...

Hi I have the following code: lorentz[A_, Ox_, Oy_, FWHM_, x_] := A (1/3.14) FWHM/((x - Ox)^2 + FWHM^2) + Oy; gauss[A_, Ox_, Oy_, x_, C_] := A Exp[-(x - Ox)^2 C] + Oy; Convolve[lorentz[1, 0, 0, 1, x], gauss[1, 0, 0, x, 1], x, y] It takes extremely long time for this to finish -- is it...

Dear all, I wonder if there is a function implemented in Mathematica to find the average value of a 2d function f(x,y). So averaged over a specified domain. Thanks in advance for any help! Sue

Hi all, I need help to edit the code so that the arrangment of the outputs match the one that Iam looking for. For instance: term[n_Integer, m_Integer] := Apply[Dot, Table[matg[[i]], {i, m + 1, n}]].hl[[m]] /; (n >= m); Subscript[lay, s_] := Sum[term[s, i], {i, s}]; In[1]...

I am trying to use Sow/Reap in Mathematica to replace Append in my code. The problem is that I need to define the tags within my code, but doing so doesn't reap the desired results. A simplified version of my code is Clear[tag]; Reap[tag = {a,b}; Sow[1, a];Sow[2,b];, tag] However, if I...

I have defined vectors which are to be column vectors of some matrix. How to tell Mathematica to construct such matrix whose column vectors are already defined vectors?

Hey guys, I'm using Mathematica to plot some graphs and I'm having a bit of a hard time. First I have to solve the following equation, 2/3 a^2 b''[a] + (1 - w[a]) a b'[a] - (1 + w[a]) (1 - 3 c w[a]) b[a] Boundary condition b'[0.0001]=0 Where; w[a_] := 2*a^(3*(1 + c))/(1 +...

Hey guys, I'm using Mathematica to plot some graphs and I'm having a bit of a hard time. First I have to solve the following equation, 2/3 a^2 b''[a] + (1 - w[a]) a b'[a] - (1 + w[a]) (1 - 3 c w[a]) b[a] Boundary condition b'[0.0001]=0 Where; w[a_] := 2*a^(3*(1 + c))/(1 +...

Hi guys, i have expressions of the type (X*a)(X*b). I want replace this by X^2(a*b). So i tried building a block which does nothing but %//.(X*a_)(X*b_)->X^2(a*b). However, this works only if a is distinct from b. If a and b are equal if HAVE to use the replacement command...

Hi I am a user of Mathematica as a tool for calculation in Structural Mechanics. I had a problem recently in the integration of expressions a bit long, and honestly I do not understand what might be causing the slowdown. I want to just evaluate this: This is for the Raileygh-Ritz...

Hi, I wrote some code on a computer at school. dim = 50; Ediff = ConstantArray[0, {dim, dim}]; \[Lambda] = .5; steps = 10; Ediff = ConstantArray[0, {dim, dim}]; L = 1; For[xs = 1, xs <= dim, xs++, For[ys = 1, ys <= dim, ys++, Ediff[[xs, ys]] = Ediff[[xs, ys]] +...

here is the pseudo java (processing) code I'm trying to run float x0,y0,x,y,t; x0 = 2; //arbitrary value y0 = 5; //arbitrary value x = x0; y = y0; for(int t = 1; t <= 100; t++){ //one hundred iterations t =atan2(y,x); x += cos(t); y += sin(t); print( x...

Is there a way to bring up all of the notebook you're currently working with in Mathematica? Sometimes, a notebook I'm using disappears and I can't call it back up.

Hi, I have some Mathematica code I want to translate in Matlab, does anyone know how to translate the Table [expr,{i,a,b}]to Matlab? I made one with a loop, but loops are very slow and I am going to analyze large datasets with the functions.

Hello, how do you define a function or make a list with n elements, where n is any? More precisely there is a function in mathematica, BellY BellY[n,k,\{x_{1},...x_{n-k+1}\}] which gives the partial Bell polynomial. I would like to define in mathematica the complete Bell polynomial defined...

Hello, I am currently working on a problem, but at the moment I am stuck. I just don't know how to solve the problem so I hope someone can help me with it. This is the question 1a : Consider the standard mass-damper-spring system: m y''+γ y'+k y=u where the constants have the following...

pl help Here I have written a mathematica code. It is a small example of my actual problem. In this code, 't' is an arbitrary constant which is being fixed at the top of the programme by assigning a value. All the other variables 'a','b','c' depend upon 't' . and 'h1','h2' are the roots of...

Hi, I really need some help here. Right now I am plotting points on a 3D Scatter plot chart in Mathematica. I want to assign each of these points with a value which will be the label. Basically each point has 4 variable in the parameter. Its x,y,z cartesian coordinate position and the last...

Say I have a function: f(x,y) I have values for f and y such that I can make a plot of f vs y... however, from the values of f and y, I am interested in finding what x is. The x value should be a constant that is optimized to fit the f vs y graph. Is there a way I can do this in...

The following creates a button window to select a notebook to run. When a button is pressed Mathematica finds the notebook but cannot evaluate it. The following error occurs Could not process unknown packet "1" CreateDialog[ DialogNotebook[{Column[{TextCell[ Row[{Style["I would like...

Good day, I am relatively new to mathematica, and I am using it to fit a formula to a dataset that I have. To begin with, I am having issues with a very simple part, namely with a ListPlot. I have an xls file with 2 columns and 26 rows, corresponding to an x and a y value. Using import, I...

This question is not directly related to QM although my reason for asking it is I'm trying to compute wavefunctions on my PC -- please excuse me if this is the wrong place to ask. I am having trouble using the DFT in Mathematica although I don't think my problem is directly Mathematica related...

Homework Statement Plot the position vs time graph of an object attached to a spring. m = 0.5 kg k = 60 N/m x[0] = 3 m v[0] = 4 m/s Homework Equations d²x/dt² + (k/m)*x = 0 The Attempt at a Solution I solved the differential equation using DSolve and plotted the resulting...

Hello, I have the following Laplace function: \frac{\left[2\sqrt{a\,s}K_1\left(2\sqrt{a\,s}\right)\right]^N}{s} where N and a are constants and K1 is the modified Bessel function of the second kind and order 1. I need to find the inverse Laplace transform at a specific value. Can I do...

I'm trying to make a little manipulate/interactive box that shows the vector divergence of the E-field coming from a sphere. I have no idea how to start as I'm really new to Mathematica. Does anyone have any pointers? I can't find anything particularly helpful on the Wolfram reference or...

Mathematica says: "NDSolve::alliv: "The function x[t] was specified without dependence on all the independent variables. Each function must depend on all the independent variables."" But function x depends only on t, there is no other variable specified in the code. What could be the...

So I want to be able to count the number of neighbour pairs of opposite sign in a list of 1's and -1's (eg. 1 for {1,1,1,-1,-1} and 3 for {1,-1,1,-1,-1}), but I can't work out the pattern matching syntax in Mathematica even after reading the documentation. I was thinking to use the Count...

How you define vector function in Mathematica? For example, f is a vector function and f=(xy,yz,zx). How to define this in Mathematica and then how to calculate the value of the components of f for any number x, y, and z? For scalar functions it goes as this: f[x_]:=x^2 f[4] Any...

All, I have a button panel (code below) where a pressed button initiates a particular Mathematica script. When the button is pressed the script runs (as it should) but the panel closes. I wish for the panel to remain open so that I can click another button. The panel should close only when...

I'm looking for a function to plot things in 3D using cylindrical coordinates but apparently Mathematica doesn't have one?

I'm using a predefinied program... but I don't know because give me this error: Clear["Global`*"] Remove["Global`*"] << "ErrorBarPlots`" xdata = {180, 200, 220, 240, 260, 280, 300, 320, 340, 360} {180, 200, 220, 240, 260, 280, 300, 320, 340, 360} ax = First[xdata] 180...

Hi I have a set of data points in units of (time, voltage), and they have the form of a Gaussian when I plot it. I would like to normalize my data set, i.e. find a factor C that I multiply on to the voltage-data such that the area is 1. However, is there a way to numerically integrate data...

I have problem with for loop in Mathematica given in the notebook attached. It evaluates only first step and then shows problems and the kernel could not stop running. What could be the problem? What is the correct syntax?

What I want to do with Mathematica is this: 1. Given are 4×1 matrices e and u, and 4×4 matrix V. Matrices e and V are constant. 2. Solve the system u=V.f for f, where f is 4×1 matrix. 3. Define 4×1 vector g such as g1=f1*Exp(e1*t), g2=f2*Exp(e2*t) etc, where g1 is the first component of g, f1...

I have the next function: z^3-2xz+y=0 and I want to find taylor expansion of z(x,y) at the point (1,1,1), obviously I need to define F(x,y,z) as above and use the implicit function theorem to calculate the derivatives of z(x,y), but I want mathematica to compute this to me. I tried the Series...

hey guys i have a problem regarding plotting several graphs in one plot via mathematica. I want this plot's x-axis to be in a log scale so i used loglinear command. here the plot is b[a] vs. a. w[a] is part of my differential equation and it depends upon a and c. i have written the program...

Hey everyone, I am struggling to properly insert my simple model syntax in mathematica for a population growth. The formula is as following: Rate of growth over time = birth rate*N - death rate*N + net migration rate*N I have been trying to add saturation to the model for the population N...

Worth it or not?: Bertrand Russell's "Principia Mathematica" Evening, As a physics undergrad, I feel that my understanding of mathematics is lacking – not in terms of how to do something, but in terms of why you'd do something. For example, why take the integral of Schrödinger's Equation...

hello guys, I have a CSV file with hyuge set of data devided into 3 cols. Y position X position and corresponding intensity values. I want to plot the image , my thesis is laser- microscopy. so a 2D plot ( density plot ) based on these values. The following is what I have tried so...

I have a function of two variables f[x_][y_]=x^2+y; d=Take[NestList[f[-2],0,10],-2]; now i want to repeat this process by changing -2 to -1.99, then -1.98 all the way up to 0 and storing values in a table, so i can later make a plot. I can't figure out how to repeat this process.

say I have a function f(x)=x^7+k I want to iterate x(n+1)=f(xn), 50 times starting with x0=1 and keep the last 10 values i want to assigh these values to some function''y'' and be able to use y to calculate this iteration for any input constant K. ie say K=0, K=0.5 etc, and then i want to store...

Hello, From some matrix operations I get "AlgebraicNumber[-4 I, {0, -(1/16)}]" as output. I need to get i/4 and I don't know how. If I write "N[AlgebraicNumber[-4 I, {0, -(1/16)}]]", I get 0. + 0.25 i but it is not what I want neither. I would like to get just i/4 as output. Can...

Hi I am trying to solve Newtons equation for a particle in (x, y, z) using NDsolve. Here I what I have so far: sol = NDSolve[{ x''[t] == acceleration[x'[t], y[t], z[t]], y''[t] == acceleration[y'[t], x[t], z[t]], z''[t] == 0, x[0] == 2, x'[0] == 0, y[0] == 0, y'[0] ==...

I need to NDsolve four differential equations. But I find my code does not work. Who can give me some suggestions? Thanks!

Hey All, I'm trying to create a 3-D rectangle in Mathematica with the following measurements: Mass M=1.5 kg, and sides of length a=10 cm (parallel to the x-axis), 2a (parallel to the y-axis), and 3a (parallel to the z-axis). Let one corner be at the origin, and let the three adjacent edges...

Can MMA take limits as 2 variables approach a value? I want to take the limit as h and k both approach 0 for f(x+h,y+k)-... Is there some trick to doing this?

Long story short I have an integral which is something like this: (in Mathematica code) Integrate[1/{(1 - b^2)*{((1 - y)*z + y)^2 + (1 - y)^2*(1 - z)^2} + 2*(b^2 + 1)*{(1 - y)*(1 - z)*((1 - y)*z + y)}}, {y, 0, 1} ] Written without assumptions. b is actually only in [0,1] and I think I...

EDIT: I solved it. Sorry for waste of time! :( Here is the function ψ[x]: 0.0000551056 a0.25 (30240. a0.5 x - 80640. a1.5 x3 + 48384. a2.5 x5 - 9216. a3.5 x7 + 512. a4.5 x9)(-((a x2)/2)) I am trying to integrate ψ[x]*ψ[x] over all space (-infinity to +infinity). I don't get any error...