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

learn more… | top users | synonyms (3)

2
votes
1answer
149 views

NDSolve doesn't work for this simple differential equations

I try to use NDSolve to do this partial differential equations, but the result is wrong... The source term is a Gaussian pulse as time and an exponential decay with the depth z. The result should be ...
4
votes
0answers
91 views

Unknown (internal?) function with `DSolve`

I am solving With Mathematica 10 the following ODE system: DSolve[{B'[x] == -f[x]*Cos[l]*G[x], G'[x] == +f[x]*Sin[l]*B[x]}, {B,G}, x] The solution is almost ...
0
votes
1answer
75 views

Solving differential equation with Variable constans?

How do I solve the following differential equations using Mathematica under the assumptions: If Abs[r[t]]<50 ->A=0.5;B=2; AND ...
6
votes
1answer
210 views

3D FEM with holes

I'm trying to create a finite element mesh with some spherical holes. To my mind the behavior is erratic but generally fails. I'm using MMA 10.0.2. In the following 3 different sets of holes are ...
-4
votes
1answer
64 views

How do I solve and graph this coupled PDE?

How can I solve and graph v1 versus t and v2 versus t for a given x? Let A = 3, p = 2, B = 4. At t = 0, v1 = 2, v2 = 2. For the graph the range of t is [0,60].
1
vote
1answer
64 views

Assymptotically Solving ODE with free parameters

I have a first order ODE, $F(y[x], y'[x], A, B) = 0$. I want to solve this numerically (with the boundary value $y[1] = i$ ) but also with the following requirements: A and B are free parameters which ...
0
votes
0answers
38 views

Extract x-values from InterpolatingFunction object? [duplicate]

Ok, So I have a numerical solution to a differential equation, and I am trying to get the data points out from it. So far, I have only been able to get the y-values with the following code: ...
1
vote
1answer
132 views

Generating a boundary mesh from a list of points

I have a list of points from which I am trying to generate a mesh: ...
37
votes
2answers
895 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 ...
2
votes
1answer
131 views

NDSolve with shooting method

Here are my equations: ...
15
votes
1answer
223 views

Eigenfunctions of the Laplacian on an arbitrary mesh

So, I've constructed a mesh over which I'd like to find eigenfunctions of Laplace's equation with a free boundary (a zero Neumann boundary condition along the edge). Mostly because I figured an ...
0
votes
2answers
67 views

Plotting direction fields along with trajectories

I am reading "Dynamics and Bifurcations" by Hale and Kocac, and in chapter one the following diagram is shown: Question is, how can I achieve the same plot using Mathematica? I've tried ...
8
votes
1answer
155 views

Using a Mathematica index as a DiscreteVariable in NDSolve when solving a coupled set of ordinary differential equations

Context Since the explanation below of the problem to be solved is lengthy, let me preamble this by saying that I have code that works to solve the problem, but I don't know whether (1) it's ...
1
vote
0answers
37 views

Using an arbitrary area bounded by lines for PDEs?

I have a set of lines I pulled from an image that make a closed loop, pretend it's a unit square for reference. Is there a way to use such an arbitrary region for solving PDEs, and is there a way to ...
0
votes
2answers
58 views

Plotting a Forced Vibration Differential Equation [closed]

I'm attempting to plot a differential equation solution but mathematica is giving me trouble. The plot is showing up blank. Here is what I have. ...
0
votes
0answers
44 views
4
votes
1answer
145 views

Solving a chemical system using NDSolve

I have to solve a set of coupled differential equations and am not sure how to proceed. It is a model of a chemical process where I'm interested to find the species concentration as $r$ varies. For a ...
1
vote
0answers
84 views

Solve PDE with DiracDelta function

This question is related to Solve wave equation using DSolve?. I think it is more appropriate to open another question on this topic. In the question, the compatibility conditions can be introduce ...
1
vote
0answers
102 views

How can I improve the solution of my PDE?

I want to improve the last code to find a better solution for umi. For now I'm using only MaxStepFraction. Also, when I put no ...
0
votes
2answers
82 views

Setting up two ordinary differential equations [closed]

Just having a total blank at the moment, I realise that maybe this question would best be suited for a math forum but i assume Mathematica is needed in order to solve these equations. So the question ...
1
vote
2answers
134 views

How to write absolute value and square root in one equation?(runge kutta) [closed]

I got some problem when I run my program. It hangs when I run it. This is my coding for runge kutta Geometric mean(RKGM). ...
1
vote
1answer
69 views

Trying to find solution for system of ordinary differential Equations [closed]

I've got 2 functions: A'[t] = (0.5/(1 + (B/k1)^n)) - 0.2*A B'[t] = (1/(1 + (A/k2)^n)) - 2*B with conditions of ...
0
votes
1answer
136 views

nonlinear differential equation

For such a equation $$-y''(x) + x^2\cdot y(x) - |y(x)|^2\cdot y(x) = 0$$ with initial conditions: $y(-500) = 0$ and $y(500) = 0$. I write in mathematica ...
1
vote
1answer
82 views

Diff eqs can be solved?

I'm trying to solve a pair of coupled differential equations, but am having trouble making DSolve evaluate my equations. ...
3
votes
3answers
127 views

Using NDSolve within Manipulate

I'm trying to use NDSolve inside manipulate. Specifically, I have a series of differential equations with coefficients, k1, k2, ...
5
votes
2answers
298 views

Solve wave equation using DSolve?

The wave equation with pined-pined boundary has well-known solutions, but directly input the equation into Mathematica does not return them. $$EI\cfrac{\partial^4 w}{\partial x^4} + ...
5
votes
1answer
195 views

How to program efficient undershoot/overshoot

I would like to solve the following boundary value problem for $y(x)$ for a fixed value of $k$ between $0 < k <1$: $$y'' + \frac{3}{x} y' - y + \frac{3}{2}y^2 - \frac{k}{2}y^3=0 \\ y'(0) = ...
1
vote
0answers
83 views

How to incorporate the boundary conditions into the differentiation scheme in MMA?

Let that we want to numerically solve the following PDE \begin{equation}\label{sde} -r V(S,t)+r S \frac{\partial V(S,t)}{\partial S}+0.5 S^2 \text{sigma}^2 \frac{\partial ^2V(S,t)}{\partial ...
5
votes
1answer
224 views

ODE fitting to dataset

So, I have a ODE system, it is a complex biochemical kinetic mechanism with six species changing over time. ...
1
vote
1answer
120 views

Power::infy: “Infinite expression 1/0. encountered” (ODE Solution using NDSolve)

While solving the following ode, Eq1 = f'''[x] - f'[x]*f'[x] + f[x]*f''[x] - K1*(2 f'[x]* f'''[x] - f[x]* f''''[x] - f''[x]*f''[x]) - f'[x] == 0 with the ...
1
vote
1answer
144 views

NDSolve and assuming

I am trying to solve the system of PDE with NDSolve and I got an error message like this NDSolve::deqn: "Equation or list of equations expected instead of ...
1
vote
1answer
49 views

Unexpected behavior of WhenEvent in NDSolve?

In Mathematica 9.0.1, here is a simple code to numerically solve an initial-value problem: ...
0
votes
1answer
76 views

Problem using NDSolveValue

Good afternoon, fellas. Could anyone here help me solve this problem? ...
0
votes
0answers
97 views

Plotting a numerical integration

I have the following code to calculate a numerical integral for any given a: ...
1
vote
2answers
92 views

Plot ParametricNDSolve

I am trying to plot the probability of a molecule being in a certain state via the numerical solution of the Schrödinger's equation of molecule in an oscillating electric field (in a rotating system). ...
4
votes
2answers
183 views

Evaluation problems with NDSolve using the fixed step ExplicitEuler-Method

Motivation I am using Mathematica to implement System Dynamics (SD) models, e.g. differential equation models as applied to management/economic problems following the "methodology" of Jay W. ...
1
vote
1answer
49 views

Solving ODE having function which changes conditionally

I want to solve the following set of ODE with a piecewise function A[t] which changes if a certain condition is satisfied. ...
1
vote
1answer
81 views

ParametricNDSolve without initial condition?

For example, ParametricNDSolve[{y'[t] == a y[t]}, y, t, {a}] returns a ParametricFunction which is actually inexplicable: ...
2
votes
0answers
157 views

Choosing FEM element type and mesh refinement

I have only recently been introduced to Mathematica's(v10.0) FEM capabilities. I understand that for solving PDEs on non-uniform shapes via NDSolve, Mathematica ...
3
votes
1answer
174 views

Adding noise to nonlinear ODE system

I have a problem with adding a noise to 3d order ODE nonlinear system. I solved the system numericaly and got the periodic solution and now I need to add random pertubations. Due to the lack of ...
0
votes
0answers
39 views

Differential equations: Solving a PDE

I need to find an analytic solution to this PDE and for some reason nothing works. This is the equation: ...
1
vote
0answers
63 views

Solving non-linear coupled ODE with constraint using NDSolve

I am trying to use NDSolve in Mathematica to solve the following coupled ODEs from a paper I am trying to understand. I have the following 5 differential equations: ...
0
votes
2answers
125 views

NDSolve issue with initial and boundary conditions

While solving the heat equation in one spatial variable $u_t = u_{xx} $ (x goes from 0 to L) with the initial temperature distribution $T_0 \frac{x(L-x)}{L^2}$ , and with neumann boundary conditions ...
0
votes
0answers
37 views

Series Expansion of Differential Operators [duplicate]

We know about the possibility of series expansion of a function in Mathematica. Same as below, where a function f[a x] is expanded around zero up to 5th order. ...
2
votes
1answer
58 views

Plot directional field with a initial vlue

I do not have enough reputation to comment on other questions, hence I have to ask my own one I want to solve the following nonlinear first order ODE, $$(1 + x) (0.5 y[x]^{-0.5} (x - y[x])^{0.5} - ...
1
vote
0answers
65 views

Orthogonality relations of Hermite polynomials

The Hermite polynomials are orthogonal. $$ \int_{-\infty}^\infty H_m(x) H_n(x) e^{-x^2}\, \mathrm{d}x = \sqrt{ \pi} 2^n n! \delta_{nm} $$ Does Mathematica not use this relationship? Because running ...
0
votes
0answers
53 views

Error when using NDSolve in SIR model

I have the following code for a SIR model: ...
0
votes
1answer
55 views

NDSolve::bcsol?

I'm writing a code to numerically solve a differential equation. It has to repeat the process several times since there is a variable involved which is unknown, and the equation is fit to some ...
0
votes
1answer
65 views

Symbolically solve system of equations involving matrix

I need to solve the following system of equations symbolically for a_i, b_i, and c_i: ...
2
votes
0answers
65 views

The idea behind Stiffness switching method with NDsolve [closed]

Does the Stiffness switching method with NDsolve switch just between multiple variants of 4th order Runge Kutta method or it uses also other methods?