Questions on the symbolic (DSolve, DifferentialRoot) and numerical (NDSolve) solutions of differential equations in Mathematica.

learn more… | top users | synonyms (3)

6
votes
1answer
236 views

Second order differential equation with boundary conditions solving repeatedly

I want to solve a second order differential equation in the interval[-1:1], which does not have a analytic solution, \begin{eqnarray} y''(x) &=& k \phi^2(x)y(x) \\ \phi(x) &=& ...
7
votes
1answer
446 views

Modeling neural excitation with a non-linear differential equation

I think I have a special problem and I am not sure how to search for an answer, so I thought I would try here. I am working with the so called FitzHugh-Nagumo model which describes very simple ...
8
votes
0answers
2k views

Integro-differential equation [closed]

I have to numerically solve a nonlinear partial integro-differential equation using Mathematica. This is my equation, $$\frac{\partial y(x,t)}{\partial t}=\int_{-\infty}^\infty K_0(|x-u|) ...
1
vote
1answer
109 views

non linear Least Squares with model function given by ODE's

I have to fit to experimental data a model function given by a linear combination of functions y1(t) and y2(t) a1*y1(t) + a2*y2(t) with a1 and a2 ...
2
votes
1answer
203 views

How can I invoke the solution of NDSolve to determine a parameter in my equation just inside NDSlove?

I am trying to solve a differential equation by NDSlove for $h(x,t)$. It reads $$h_t=h_{xx}-V_h-\lambda(t)$$ where $V_h$ is a given function of $h(x,t)$ denoted by ...
1
vote
0answers
34 views

NDSolve error with a given system of equation [closed]

I know that my question will be very similar to NDSolve giving NDSolve::deqn: message when given my system of equations . Sorry for that, I am so new on Mathematica. I have tried to use the ...
3
votes
2answers
111 views

FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions

I found the deflection of an orifice plate (circular plate with a hole) subject to uniform pressure using Mathematica's NDSolve functionality. The plate is fixed at its outer circumference and is ...
0
votes
0answers
55 views

NDSolve fixed step problem

Working example here: Drive folder (have both files in the same directory! Notice: the line <<variables' in the file seems to throw an error for me, but ...
3
votes
1answer
47 views

DSolve: Why does Mathematica does a separate integration from 0 instead of 1 and then from 1 to upper limit?

I am trying DSolve in Mathematica for an ODE. Here is the command I am trying and the result I get DSolve[{x'[t] + y[t] x[t] == 0, x[0] == x0}, x[t], t] $$ ...
5
votes
1answer
125 views

List of unknowns in a differential equation

I often find myself wanting to see what variables/constants are undefined in a large set of differential equations, is there a function that will do this for me? As an example, for the following ODE ...
2
votes
1answer
123 views

Solving SDE: $\frac{dy(t)}{dt}=(c+\sigma_wW(t))y(t)+\epsilon(t) $ in Mathematica

I want to solve this differential equation $\frac{dy(t)}{dt}=(c+\sigma_w W(t))y(t)+\epsilon(t) $. For details see ...
-1
votes
0answers
64 views

Solving a system of first order differential equation with DSolve

How to solve this system of differential equations analytically using Mathematica (using DSolve) $$x'(t)=\frac{a+s}{(1-t)d}x(t)-\frac{b}{(1-t)d}y(t)$$ and ...
2
votes
1answer
79 views

Fails to apply NonlinearModelFit on a numerically evaluated model defined as a system of ODEs

I call for your help with because I'm getting troubles to fit a model defined through a system of ODEs. The system of ODEs is as follows: ...
6
votes
1answer
94 views

Why should the spatial derivative order of the ODE *not* exceed two?

Following this question I came across this strange behaviour. Let me define a 1 D interval implicitely ...
11
votes
2answers
122 views

Integration of a PDE solution over a custom-defined subdomain

Preamble: In my previous question I raised a problem of a custom-chosen integration of the solution of a PDE. In that case it was an integration along a line lying within the PDE domain. The solution ...
1
vote
1answer
44 views

solving coupled differential equations with parameters subject to optimization problem

I have an equation, but I don’t know how to program it in Mathematica: x[0]==1; y[0]==1; x'[t]==(a+b) x[t] y[t]; y'[t]==a b y[t]/x[t]; where ...
0
votes
1answer
144 views

Numerical solutions of Gross-Pitaevskii equation

I am trying to solve the Gross-Pitaevskii equation numerically and I am trying to use NDSolve. Nevertheless it is become a dead end. The problem is boundary and ...
2
votes
2answers
119 views

How can I check that the general solution obeys the original differential equation?

For example: I solved the following differential equation: y''[x] == (λ x^(3/4) y[x])/Sqrt[1-x] that 0<x<1 with the ...
0
votes
2answers
187 views

Solving differential equation with initial conditions

I want to solve the following differential equation with initial conditions: ...
1
vote
1answer
30 views

Reverse order in evaluated list

I solve the following set of ODE and plot solution of S1 vs k2 (Figure 1). ...
2
votes
1answer
68 views

Increasing NDSolve computation range for 3D stationary PDE

I am trying to solve 3D Schrodinger equation (with mass equals to 1) and with no potential (I will add a potential later). My goal is to set boundary conditions on the outer boundary of a sphere and ...
1
vote
1answer
697 views

Solving a complex-valued differential equation with NDSolve

I am trying to solve $dx/dt=\sqrt{1+(ix)^{1.8}}$ for initial condition $x[0] =-0.9877 + i 0.1563$, where $x$ is a complex variable. I would like to plot the imaginary part of the solution versus the ...
0
votes
0answers
224 views

Differential equations with a complex variable

Is Mathematica able to handle ordinary differential equations where the variable itself is complex? I am looking for solutions of ODE systems of the form $$\left\{ \begin{align} i\frac{da_1}{dt} ...
2
votes
0answers
40 views

NDSolve with NIntegrate [duplicate]

I'm new to Mathematica and I'm trying to use it to solve some equation like: NDSolve[x'[t]+NIntegrate[x[tau],{tau,0,t}]==20,x,{t,0,10}] But it keeps giving me ...
1
vote
1answer
47 views

Matrix issue with NDSolve

I am solving a differential matrix equation with NDSolve. ...
8
votes
1answer
415 views

Optimizing Monte Carlo simulation of a Pred-Prey model

My assignment and code As part of an assignment for one of my classes, I'm trying to run a "massive" Monte Carlo simulation in Parallel on the follow model: ...
10
votes
3answers
2k views

help to plot Poincaré section for double pendulum

I am reading a book about classical mechanics. In the chapter about chaos, it gives the simplified and scaled equations for double pendulum as $$ \frac{d}{dt}\left[ \begin{matrix} \alpha \\[3mm] ...
2
votes
1answer
89 views

Problem in understanding PDE solution

I solved the PDE shown below and I got the solution. I am new to this tool. In the solution, I have a constant multiplied by a vector like function (please refer second term of the solution). What ...
1
vote
1answer
46 views

Trouble solving a system of partial differential equations

I am trying to solve a system of linear PDEs analytically and I dont get the answer. I am new to this tool, appreciate your help ...
6
votes
1answer
304 views

Dipolar magnetic field lines inside a cylinder

I'm drawing the magnetic field lines of a rotating dipole (magnetic field distorded by emission of radiation). Currently, my Mathematica 7.0 code is fully working, but I'm having a constraint which I ...
9
votes
1answer
127 views

Kernel Crash while using “ValuesOnGrid” method of InterpolatingFunction

Bug introduced in 10.0.0 and fixed in 10.2.0 I was experiencing some kernel crash while using "ValuesOnGrid" method of an ...
18
votes
2answers
283 views

Getting rid of spikes in the PDE solution

Bug introduced in 10.0 and persisting through 10.2 or later Note: In 10.0, Rationalize[fd, 0] was needed or mesh generation would fail. Preamble: I am solving a ...
9
votes
1answer
101 views

MMA breaks down when `DSolve` -ing a third-order linear ODE

My code is: DSolve[y'''[x] + (x^3 + a^3) y[x] == 0, y[x], x] It should be noted that similar code not only meets the similar problem but also returns more ...
4
votes
1answer
153 views

OutputResponse in V9 vs V10

Bug introduced in 10.0.1 and persisting through 10.2 or later I have same code that works in version 9.0.1, but not in version 10.0.1.0 (according to @Sektor, it works in 10.0.0.0, too). ...
5
votes
1answer
77 views

In a PDE solution let us integrate out one of two coordinates

Here is a trapezoid within which I am solving a simple Laplace equation with the Dirichlet boundary conditions set at its top and bottom: ...
1
vote
0answers
44 views

Transform a vector DE into a set of scalar DEs [closed]

I have the vector equality: {{a,b},{c,d}}.D[{x1[t],x2[t]},t]=={{e,f},{g,h}}.{u1[t],u2[t]} I want to solve this using DSolve. ...
37
votes
2answers
893 views

Numerically solving Helmholtz equation in 3D for arbitrary shapes

Context While studying manifold Learning I got interested in finding the eigenvectors of the Laplacian. (also in connection to this problem of solving the heat equation) Following this and that ...
0
votes
1answer
52 views

Problem with solving coupled differential equations

The following code seem to produce an output that appears to be independent of the parameter "fg". It is possible that the way the code is written might be causing it, but how can one fix this ...
12
votes
1answer
705 views

Numerically solve 2nd order differential equation with singularity

Consider a second order differential equation with a potential that diverges at some generic value in the variable. For example: $$-y^{\prime\prime}(s)+\frac1{\mathrm{cn}{(s\mid k^2)}}y(s)=0$$ where ...
1
vote
1answer
780 views

Integro-differential eqn with double integral

I am looking at the following variation of a integro-differential, with y[0]=1. The output is not great, any solutions to this? ...
0
votes
1answer
63 views

NDSolve There are fewer dependent variables, {V[S,t]}, than equation

I am trying to solve the BlackScholes PDE for Barrier option. It works fine for european barrier, but errors out on american boundary condition ...
11
votes
0answers
186 views

Possible bug with FourierTransform linearity

Each component is easily transformed but the sum is not: FourierTransform[f''[x], x, k] FourierTransform[f'[x], x, k] ...
7
votes
1answer
256 views

Problem with Neumann condition in quarter disc

Bug introduced in 10.0 and fixed in 10.2 So I'm following the available examples in version 10 for FEM, The plane stress operator is shown as this ...
6
votes
1answer
98 views

code is ignoring WhenEvent

I would like the value of x[t] to be equal to 0 when it has a negative value so I used WhenEvent but it has no effect on my code ...
0
votes
0answers
42 views

Is there a way to pre-calculate a slow animation and then run easily with no lag? [duplicate]

I have an animation of a triple pendulum on a moving base, however even after trying to change the whole coding to make it faster by using module and all it is still pretty bad as I need it to be ...
-2
votes
2answers
81 views

Error in the solving of the simultaneous Ordinary Differential Equation [closed]

As I am trying to get the solution of the simultaneous ordinary differential equations ...