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

learn more… | top users | synonyms

1
vote
1answer
133 views

Runge-Kutta 2nd Order ODE Solver

Suppose I have a 2nd order ODE of the form y''(t) = 1/y with y(0) = 0 and y'(0) = 10, and ...
0
votes
1answer
97 views
0
votes
2answers
110 views

Plot slope field in Mathematica for a difficult differential equation

I am trying to plot slope field for $\ddot{a}(t)$ for this differential equation but I don't know how to $\frac{d}{dt}\ddot{a}(t)=2\frac{\dot{a}(t)}{a(t)}\ddot{a}(t)$ I want to show that already ...
2
votes
1answer
100 views

Finding the roots of a function with a parameter

I am using a differential equation method for finding the roots of a function. My function is much more complicated but I can illustrate my problem with a trivial example. The problem occurs when ...
4
votes
3answers
340 views

Creating a 3D List Line Plot From Discrete Points

Given the following Runge-Kutta ODE solver and the graphical output below, how do I get a 3D line plot instead of a 3D point plot? I see that there is no ListLinePlot3D function, so I thought it might ...
1
vote
2answers
92 views

How to set initial guess in NDSolve

I'm trying to solve a BVP using NDSolve and I want to impose the starting initial guess. I googled and looked at the help, but I couldn't find this option. ...
0
votes
0answers
63 views

System of differential equations - solving for a constant

I will reformulate my question in the hopes it will be more intelligible. First, I defined the functions: ...
1
vote
1answer
91 views

Solving differential equation - initial condition

I want to solve the differential equation: DSolve[{R''[ρ] + 2 /ρ R'[ρ] + (1 - (l (l + 1))/ρ^2) R[ρ] == 0}, R[ρ], ρ] This works fine and outputs: ...
0
votes
0answers
66 views

How can I get the time steps chosen by NDSolve when the Method -> “ExplicitRungeKutta” option is given?

I am trying to solve an ODE with NDSolve using the "ExplicitRungeKutta" method. I need to know exactly which time steps ...
0
votes
1answer
83 views

How do I fit NDSolve result to experimental data? [duplicate]

I have two coupled ordinary differential equations which I solve numerically in Mathematica, but now I want to fit the solution with experimental data. I do following, (the experimental data ...
0
votes
0answers
43 views

How to use results of Dsolve and get derivative of that? [duplicate]

How can I use the result of DSolve and differentiate it' result? Here for example I want to, first find urf from solving: ...
0
votes
0answers
46 views

Calculate inverse laplace transform by integration with four singular points?

I have a function U[s,x] to do inverse Laplace transform for analytical solution, as following ...
2
votes
1answer
77 views

NDSolve::ndsz problem

I'm using NDSolve to solve a set of coupled differential equation which depend on a variable x. I noticed that when I set the range of x from a small value to a large value, I obtain solutions. But ...
0
votes
1answer
47 views

Using variable with delay in WhenEvents in NDSolve

I have a system of differential equations which is solved correctly. However, I want to trigger two events, one of which depends on the "time elapsed" after of the first event trigger. Is it possible ...
0
votes
0answers
58 views

Solving a system of three differential equations

I'm trying to create this coloured graph which I'm going to animate, but I'm stuck with solving a system of three differential equations: ...
0
votes
1answer
108 views

Better way to check conservation in NDSolve

The partial differential equation in the NDSolve below is conservative (by construction) of $\int_{-1}^{+1} p(t,z)\ dz$ for all $t$. ...
0
votes
0answers
52 views

Can't Solve a system of differential equations, Solve takes too long

My question is probably not really interesting but I really need help on this because I won't have access to a more powerful computer. So my model contains 7 equations : ...
1
vote
0answers
38 views

RSolve works slow and does not give an analytical solution for solving high order linear difference equations

RSolve@@{{a[n]-a[n-10]==101}~Join~Table[a[i+1]==101+10i,{i,0,9}],a[n],n} I want to get the exact form of general formulas, but ...
6
votes
1answer
561 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 ...
0
votes
1answer
136 views

Zooming into Parametric Animation Plot that Uses InterpolatingFunctions

I have an interplanetary trajectory solver that calculates a spacecraft route from Earth to Mars, and I am trying to zoom into the spacecraft when it reaches its Mars intercept orbit but am having ...
1
vote
0answers
47 views

Combined Elliptic-Parabolic PDE in 1D [closed]

Here is an example of solving a system of parabolic and elliptic PDE: http://reference.wolfram.com/mathematica/tutorial/NDSolveDAE.html (see the subsection "Combined Elliptic-Parabolic PDE in 1D"). ...
0
votes
0answers
49 views

PDE combined with ODE 1D

I try to solve the following system of PDE coupled with ODE: $$\theta_t - a\theta_{xx} + b\kappa_a(\theta^4-\varphi)=0,$$ $$-\alpha\varphi_{xx} + \kappa_a(\varphi - \theta^4) = 0,$$ $$-a\theta_x + ...
11
votes
3answers
416 views

Determine frequency of oscillations

I am wondering how I could determine the frequency of oscillations of a differential model equation? How could I find the frequency from this example given in Mathematica Documentation: ...
0
votes
0answers
44 views

Ideas for NDSolve?

I'm currently trying to find a numerical solution to a differential equation of the form: D[W[X], {X, 4}] ==(-(1/(delta + G - (G X)/L)^2) + 1/(delta + (G X)/L)^2) ...
1
vote
1answer
101 views

Bilateral Laplace Transform in Mathematica (to be used with Moment Generating Functions)?

I am trying to solve a linear ODE for a Kolmogorov forward equation to get a stationary distribution of a random variable. The easiest approach may be to transform the ODE with a two-sided Laplace ...
0
votes
0answers
69 views

Problem with using the results of ParametricNDSolve in a ParametricPlot3D

I new to Mathematica and I have a problem with certain equations. I want to plot a trajectory in spherical coordinates given by the equations: ...
0
votes
0answers
77 views

Why do identifiers of the form ${\tt c}k\$nnn$ appear in the rules returned by my function?

I want to define a function, PolynomialInterpolation to find the coefficients of a polynomial interpolation function. My trial is as shown below: ...
2
votes
2answers
122 views

Find minimum value of function which calls NDsolve [duplicate]

I have defined a function which takes a coefficient for a differential equation as an argument, solves the differential equation with NDSolve, and returns the value of the solution at time t=10. I can ...
1
vote
1answer
114 views

How to choose a solution of a differential equation?

I am trying to solve the following differential equation: ...
2
votes
0answers
66 views

How to specify initial condition including integral equations in the case of the dirac equations?

Hi I am trying to solve a system of differential equation with NDSolve. The problem is I have like $3$ (maybe $5$) boundary conditions, but only $2$ differential ...
1
vote
0answers
61 views

Stiffness or Singularity

So trying to run some simulations based on a model I found in literature. I am by no means a Mathematica expert, or a mathematician for that matter. Initial I was getting a ...
4
votes
0answers
66 views

Why do NDSolve and OutputResponse not evaluate non-analytic functions numerically?

test = OutputResponse[TransferFunctionModel[1/(1 + s), s], Exp[-(1/t)], {t, 0, 10}] does not evaluate. If Exp[-(1/t)] is ...
0
votes
0answers
73 views

Nested For loop with NDsolve

I am solving two simple differential equation using NDsolve. ...
0
votes
0answers
34 views

How can I draw a bifurcation for the following system [duplicate]

How can I draw a bifurcation for the following system between b and T
3
votes
2answers
219 views

Symbolically solve a system of coupled second order differential equations

I have posted a similar question in another forum where the general consensus seems to suggest that it is not possible to symbolic solve a system of coupled second order differential equations with ...
0
votes
0answers
42 views

Problem with NDSolve and NDSolve::ivone error

I was trying to solve an diffusion-reaction rate model for the oxidation of Silicon (based on the nice paper by Almeida et.al.) with the help of NDSolve. In this model the boundary conditions are a ...
0
votes
1answer
180 views

Parallel computation of NDSolve and a problem

I used this code to solve a system of differential equations: ...
1
vote
1answer
76 views

How to study the sign of an expression including a function, with conditions on the function?

I'm trying to study the sign of an expression, having information on some of my parameters. I don't know how to code these conditions. For example, I want to know when this expression: ...
3
votes
1answer
113 views

How to change a parameter dynamically inside NDSolve

I have no idea how to dynamically update a parameter inside NDSolve. For example consider following first-order DE: ...
1
vote
1answer
139 views

To solve Cylinderical PDE

I want to solve this PDE. I have tried to solve it with NDSolve but found error 'Boundary values may only be specified for one independent variable. Initial ...
0
votes
1answer
138 views

Parallel computation of a system of differential equation

I need to solve a system of differential equations which contains, at minimum, ten thousands of differential equations. Actually, I am trying to solve density matrix equation. It's like something ...
4
votes
0answers
177 views

Modelling Hysteresis with a Differential Equation

I want to implement the bulk ferromagnetic hysteresis model (mostly the Jiles-Atherton Model), see http://drum.lib.umd.edu/bitstream/1903/6043/1/PhD_99-1.pdf page 44 equation (30). The needed ...
5
votes
2answers
481 views

Solving the Helmholtz equation in polar coordinates

I'm trying to (as a simple starter) find the solution to Helmholtz equation in polar coordinates -- I already know what it is and can derive it by hand, but just want to start here before asking ...
4
votes
1answer
114 views

NDSolve - sampling for result during the computation

I am using NDSolve for a Langevin dynamics problem. I want to the know long term behaviour of my system ($t>1$) but it has to be simulated with very small time steps ($dt\sim 10^{-9}$). An example ...
0
votes
0answers
61 views

Plotting the solutions of a system of differential equations on a simplex

I want to solve and plot the solutions of the following system on the triangle simplex: x'[t] == x[t]*(2*x[t] - 2) y'[t] == y[t]*(2*x[t]) z'[t] == z[t]*(2*x[t]) ...
0
votes
1answer
178 views

DSolve will not apply assumption `m ∈ Integers`

I am trying to solve a linear second order ODE using DSolve which involves an arbitrary integer m. ...
2
votes
1answer
140 views

Understanding of method for NDSolve

I used automatic method for NDSolve. Then I asked myself - which method Mathemathica prefered? I got an answer for this question on this forum, that I need to use: ...
2
votes
1answer
215 views

Different Poincaré maps for same Equation

I've seen the Wolfram Poincaré example. I'd like to modify this sample in order to have a list of Poincaré maps. In the example the event is captured only when ...
0
votes
0answers
47 views

Automatic method for NDSolve [duplicate]

How can I understand which method Mathematica automatically chooses to use for numerically solving the following two equations? 1: ...
0
votes
0answers
125 views

Calculating the Lyapunov Exponent/Spectrum for the Duffing Oscillator

Some fellow students and I are trying to calculate and plot the Lyapunov spectrum for the Duffing oscillator (currently a test approach). Each of us has Mathematica knowledge based purely on past ...