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

learn more… | top users | synonyms (2)

11
votes
2answers
639 views

Solving and Animating Three Body Problem

I am attempting to solve a three-body problem using the Lagrange formalism, with a 1/r potential. I started off by defining T and U (kinetic and potential energy) as follows: ...
1
vote
2answers
61 views

Minimize the sum-of-squares errors for a system of ODEs

I'm trying to minimizing the sum of squared errors on a 2-equation system of ODEs and empirical data. I keep getting an error and I think the root of the problem is my invocation of ...
2
votes
0answers
48 views

Problem solving a nonlinear partial differential diffusion equation

EDIT: actualy the nonlinear partial differential equations for interacting density distributions, including boundary conditions, should be given as $$ \frac{\partial\phi}{\partial t} = D ...
0
votes
1answer
152 views

A good choice to find solutions of ODE

In order to rephrase the question. I would ask you to take a look at the following problem: I have this system of ODE: $$ \begin{cases} \dot{x}(t)= x(t)[r_1(t) -b_1(t)x(t- \tau(t)) - c_1(t)y^m(t)] ...
0
votes
0answers
21 views

Singularity or stiffness error solving a system of differential equations

I'm trying to solve a system of differential equation relating to the elastic problem of a tapered beam. Here is my code: ...
0
votes
0answers
31 views

Why does DSolve not give the expected analytical solution for a system of DAE?

I am trying to solve the following system of DAE: X'[t] - y*Y'[t] + x*X[t] == 0 Z'[t] - y*Y[t] + x*Z[t] == 0 X[t] + Y[t] + Z[t] == 1 with initial condition: ...
1
vote
1answer
181 views

Locating Periodic Orbits

Here is the code for the numerical integration of an orbit. First the module for the definition of the equations of motion. ...
1
vote
1answer
157 views

Runge-Kutta Butcher tables

I like to have the Butcher's table for Explicit (or implicit as well) Runge-Kutta method of a fixed order. I do not understand reading ...
0
votes
0answers
16 views

How to compute DSolve or NDSolve for kind of piece wise function ?

I need to Compute the expression for two functions. q[t] and q1[t] over certain range of time. The code is given below. ...
0
votes
1answer
32 views

Solving a system of ode and plotting the nullclines

I have the following problem to solve it is a system of two ode's. I'm given that q=2/Sqrt[1 + 1600 y] and then the ode's are as follows. ...
1
vote
1answer
732 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
0answers
48 views

System of non linear differential equations with two unknowns

I am new to Mathematica and am trying to solve a system of 3 (non-linear) equations with two variables $x,\ y$ and one unknown function $f(x,y)$ of these two variables. The third equation defines the ...
0
votes
1answer
56 views

Problem in computing Dsolve

I want to compute my result using Dsolve. The code is given below. ...
0
votes
0answers
96 views

Using Runge-Kutta method to find arc length [duplicate]

I don't know much about using the Runge-Kutta method in Mathematica, so I have to ask a straightforward question: How can I use Runge-Kutta to numerically approximate the parameter t of this curve in ...
-1
votes
1answer
76 views

Laplace equation

I tried to run this code, but it doesn't work for NDSolve ...
1
vote
1answer
103 views

Why DSolve[] does not work?

I try to find the analytical solution to the ordinary differential equations as fellow: $$\frac{d^2 \theta}{d t^2} = \dfrac{-M\cos\theta\sin\theta(\frac{d\theta}{dt})^2 - ...
0
votes
2answers
81 views

differential equation of order2

I want to solve the following differential equation ode = x*y''[x] + y'[x] + 4*x^2*(x^2 - 1) y[x] == 0 But do not know how to actually solve it. Any suggestion?
0
votes
0answers
54 views

fourier frequency analysis for NDSolve of complex variable

i have a set of coupled equations for complex amplitudes c1 and c2, for which i'd like to find frequency and phases (and phase shift); i looked up for some nswers on the forum, but can't seem to find ...
0
votes
0answers
48 views
1
vote
1answer
88 views

Failed to create a special random-disturbed initial condition for partial differential equation

I try to use BSplineFunction to creat a initial condition for NDSolve which has periodic boundary condition, that is realized by ...
9
votes
1answer
611 views

NDSolve::ndcf: Repeated convergence test failure. How to solve?

I am trying to simulate a system of $n$ pendulums with some friction in Mathematica 9. This is the code I am using: ...
0
votes
0answers
59 views

Black Scholes PDE solved by method of lines

I've been trying to solve the Black Scholes partial differential equation (BSPDE) by transforming it to the heat equation (HE) form and using method of lines (discretizing spatial argument and leaving ...
0
votes
0answers
29 views
2
votes
4answers
226 views

Solving differential equation

I want to solve the following differential equation DSolve[y''[x] == (λ y[x])/Sqrt[-1 + x], y[x],x] But do not know how to actually solve it. Any suggestion?
1
vote
1answer
66 views

Solve system of equations with implicit function [on hold]

I am trying to solve a system of 3 (non-linear) equations with two variables $x,\ y$ and one unknown function $f(x,y)$ of these two variables. Also, the system includes the change of this function $f$ ...
0
votes
1answer
83 views

Solving differential equation with initial conditions

I want to solve the following differential equation with initial conditions: ...
18
votes
5answers
1k views

How to splice together several instances of InterpolatingFunction?

I have a set of InterpolatingFunction returned by NDSolve which are valid over different (but overall continuous) domains. How ...
0
votes
1answer
45 views

Solve system of differential equations and show the chaotic behavior of the solutions [on hold]

I want to solve a system of differential equations and show the chaotic behavior of the solutions. Here x[t], y[t] and ...
0
votes
1answer
63 views

Inverse Laplace Transform difficulty

I am asked to find the inverse Laplace transform of: $$F(s)=\frac{2s^2+s+13}{(s-1)((s+1)^2+4)}$$ I did the partial fraction decomposition by hand and got: $$F(s)=\frac{2}{s-1}-\frac{3}{(s+1)^2+4}$$ ...
1
vote
2answers
188 views

How to solve for “k”

So the question is given that the atmospheric pressure is modelled by the differential equation p'[h]=-kp[h], p[0]=1013, p[20]=50. find k and comstant of integration. My code was: ...
0
votes
0answers
16 views

ImplicitRegion Error in NDSolve?

I'd like to solve a differential equation with some initial conditions to get an eclipse curve for Earth's orbit using Newton's Law in a cylindrical coordinate. NDSolve gave me two more useless ...
0
votes
2answers
31 views

solving ode and plotting solution vs time

I'm currently trying to solve this system of two ode's and can't seem to be able to solve them. I am new to coding, so I don't know where I could have done a mistake. I'm attaching the code which I ...
5
votes
1answer
109 views

Difficulty in getting correct Gaussian curve for diffusion of point source

I want to solve diffusion of a point source numerically and check it against analytical solution. first I define initial profile, ...
2
votes
2answers
318 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
1answer
77 views

My simulation is not never finishing and it is taking so much time?

I am running simple code which involves calculation of two first order differential equations.It is taking so much time and i am unable to get the results as the simulation is never finishing. ...
0
votes
0answers
17 views

How to solve coupled first-order ODE with vectors?

m * v'[t] = q * (E + v x B) using the following: B = (0,0,B0), E=(Eperp,0,Ez). How would I find the general solution of ...
4
votes
1answer
133 views

Numerical solution of IVP for linear ODE with variable coefficient runs wild soon

Cross posted in scicomp.SE. A friend of mine showed me this initial value problem (IVP) for a linear ordinary differential equation (ODE) with variable coefficient: $$y''(x)=\left(x^2-1\right) ...
0
votes
1answer
47 views

DSolve not outputting any result

I'm trying to solve a system of differential equations using DSolve, but I am not getting any output beyond a repetition of my input: ...
-1
votes
2answers
82 views

How to make sure input is an ODE

I'm having trouble wrapping my head around why this isn't an ordinary differential equation. ...
0
votes
1answer
55 views

How To Use D[ ] to Define a System of PDEs

We have two variables $u_{1}$ and $u_{2}$ which are functions of $x,t$. That is, $u_{1}=u_{1}(x,t)$ and $u_{2}=u_{2}(x,t)$. Consider the the system of partial differential equations $$\frac{\partial ...
5
votes
1answer
133 views

How to plot the solution of a Partial Differential Equation?

My attempt. I need to solve numerically the Complex Ginzburg-Laudau Equation (CGLE): $$ \frac{\partial A}{\partial t}=\epsilon A-(1+i\beta)|A|^2A+(1+i\alpha)\nabla^2A $$ I'm using a uniform initial ...
1
vote
1answer
116 views

Solving coupled differential equations with an eigenvalue

I have 2 coupled differential equations with an eigenvalue Ei and want to solve them ...
1
vote
1answer
74 views

Constants from differential equations

I have a problem, I've solved 2 diff equations and I have: ...
0
votes
1answer
47 views

How to properly use conditionals inside NDSolve

Currently I am working with Mathematica to learn the program, and I'm confused what exactly is wrong with this code. ...
2
votes
3answers
548 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
0answers
31 views

Difference between DSolve and NDSolve [closed]

What's the difference between DSolve and NDSolve? I know both can be used but I'm not sure of the exact difference.
1
vote
3answers
96 views

DSolve with assumptions

I'm new to mathematica and was wondering if someone could help me. I have a differential equation which is: $\frac{dy}{dx} = 1/2 y(1 - y)(a - y + a y)$ It has initial condition $y(0) = y0$ and ...
4
votes
3answers
121 views

How to fit parameters in a system of coupled differential equations

Problem description: I'm studying a chemical reaction that can be modeled by a system of coupled differential equations, and I want to use measured data to determine the parameters that appear in the ...
1
vote
3answers
71 views

NDSolve parametric plot

Here is my code. s = NDSolve[{f'''[x] + f[x]*f''[x] == 0, f[0] == A, f'[0] == 0, f'[10] == 1}, f, x] I want to plot f''[0] ...
18
votes
1answer
436 views

Analogue for Maple's dchange - change of variables in differential expressions

Maple owns an interesting function called dchange which can change the variables of differential equations, but there seems to be no such function in Mathematica. ...