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

learn more… | top users | synonyms

1
vote
1answer
96 views

Mathematica doesn't know the answer to this differential equation system

I am having some trouble with solving a somewhat heavy differential equation system, which is consisted of 10 variables and other abstract parameters as follows: ...
0
votes
0answers
34 views

Finding derivatives in a differential equation system

I have a differential equation system with 10 variables whose functional forms are: $\lambda_R = \alpha_0 + \alpha_1 P_R$ $E_R = \theta_0 + \theta_1 P_R$ $\lambda_X = \beta_0 + \beta_1 P_X$ $L_D = ...
1
vote
1answer
56 views

Code for time derivative or time differential equation

I have three variables $X, Y, Z$ whose functional form is the following (I also have a function for $Z$ which is somewhat complicated but is not essential for my question here; so it is omitted.): ...
1
vote
1answer
111 views

Definite Integral over Bessel Function

Hello I am interested in evaluating the following integral. ...
2
votes
0answers
89 views

Puzzling NDSolve[] behavior for PDE (smooth solution, inconsistent with boundary conditions)

Consider the following: NDSolve[{D[z[x, y], x, x] + D[z[x, y], y, y] == 0, z[x, 0] == Sin[x], z[0, y] == Cos[y]}, z[x, y], x, y] {{z[x, y] -> ...
0
votes
0answers
61 views

NDSolve with mixed boundary conditions

I am trying to solve a spherical Laplace equation subjected to Neumann boundary conditions in the lower half ($θ < π/2$) and constant potential in the other half (say, $f(R,\, θ) = 0,\, {\rm if}\ θ ...
2
votes
1answer
223 views

ParametricNDSolve KKT Constraint Directly

I am trying to directly solve for an optimal control using Pontryagin Maximum, where H is the Hamiltonian, in the presence of a KKT constraint. The constraint is ignored. I am wondering if a direct ...
1
vote
1answer
142 views

Smoothing a single rectangular pulse in NDSolve

Hello I am interested in smoothing out the sharp edges of a single rectangular wave centered at the origin. I am looking for a smooth function that looks almost exactly the same as the square wave ...
3
votes
0answers
65 views

Increasing MaxPoints in NDSolve results in memory issue

I am interested in increasing "MaxPoints" in NDSolve's "MethodOfLines" in attempt to increase the resolution of the plot of the solution of linearly damped wave equation with transparent boundary ...
2
votes
0answers
167 views

How to solve a nonlinear coupled PDE with initial and some boundary values

I would like to solve the following nonlinear coupled PDE with a mix of initial conditions and boundary values: ...
0
votes
0answers
65 views

No solution with DSOLVE

I have been trying to solve the following system of ODEs: \begin{eqnarray} ...
0
votes
1answer
67 views

Substitution of a function in a differential equation

Suppose we have a differential equation as follows: g[t_] := t/2 m[t_] := t/3 A[t_] := {{g[t], 1}, {1, m[t]}} p'[t] == A[t]. p[t], p[0] == {1, 0} ...
3
votes
5answers
236 views

Output is Input for a Differential with Sign?

Here's the problem: For one of my classes, we're supposed to use Mathematica to solve the equation y''=y-b*y' for a variety of b values and specified boundary conditions and plot the outputs. There I ...
5
votes
1answer
158 views

NDSolve and Piecewise functions

I want to clarify upfront, This is not a please help me get the answer to my homework question. I know what I need to do I just do not know how to get mathematica to do it. I am working on an ...
1
vote
0answers
33 views

How to update Plot and some output inside NDSolve callback (StepMonitor)

I have the following problem. I am calling NDSolve to solve a system of first order differential equations (nothing fancy), except that the number of equations measures in thousands and the number ...
2
votes
1answer
68 views

Question about using NDSolve with a matrix equation

I tried to use NDSolve to study a matrix-valued differential backward equation. The following program seems to work: ...
1
vote
1answer
128 views

NDSolve giving the wrong solution?

I'm considering the non-linear second order ODE DE $=0$, with DE given by ...
1
vote
0answers
94 views

NDSolve Wave Equation - Triangular Wave Pulse Inital Condition

I am trying to solve a simple damped wave equation with transparent boundary conditions with triangular shaped piecewise function as my initial condition. I understand that there are issues with this ...
1
vote
1answer
121 views

Can I use the values of interpolating points from NDSolve's output? And How?‎

The output of NDSolve is a (list of)InterpolatingFunction(s). In standard output format the ‎elements of interpolation are not ...
2
votes
2answers
108 views

Can Mathematica solve this equation with NDSolve?

Here is my code: ...
0
votes
0answers
48 views

How to remove error at large time in NDSolve

I am kind of amazed by the fact that Mathematica V9 can't solve a linearly damped wave equation at large time t=200 while Mathematica V8 does without giving me a huge deviation. Intuitively, an ...
4
votes
2answers
144 views

Coordinate free differential forms package

I am looking for a package in Mathematica which can handle differential forms in a coordinate free manner. I am aware of several packages which do differential forms, but it seems that for all of them ...
1
vote
2answers
104 views

How to explain NDSolve that functions must be positive only?

I have a simple system of differential equations to solve but both functions should be positive: x[t] > 0 and y[t] > 0 for all t. I was not able to explain NDSolve this fact so it goes x[t] < 0 ...
2
votes
1answer
182 views

How to solve time dependent Optical Bloch Equations for a three level system?

I am trying to solve the optical block equations for three level system which are time dependent, which i generated using the "Atomic Density Matrix" package. When i use the command to solve this 9 ...
2
votes
2answers
147 views

How to get the series coefficient of a function defined by a differential equation [duplicate]

I need a Taylor series approximation for $x(t)$, which is defined by the following differential equation. $\frac{dx}{dt}=-k_1 x+(1-x)k_2 e^{-k_3 t}$, $x(0)=0$ ...
3
votes
1answer
144 views

How to propel the integration of time a little bit further? Numerical solution can not evolve to the max time

I try to solve a nonlinear partial differential equation. I obtain a numerical solution which can not continue to the max time I set, I always receive message NDSolve::ndcf: Repeated convergence ...
2
votes
1answer
157 views

Coupled Nonlinear Differential Equations Problem

I have a big problem solving a set of coupled nonlinear differential equations using NDSolve. Solving the equations by themselves works quite well but if I want ...
3
votes
1answer
211 views

Second order differential equation

Hey guys I need someone to give me a hand on this. I don't know if this is too complicated or it's just the lack of knowledge I have on Mathematica. I'm trying to solve the following equation but ...
9
votes
3answers
198 views

Wrong answer from DSolve?

I was trying to solve the initial value problem $$u'(t) = \sqrt{u(t)} + \frac{1}{n+1}, \, u(0) = 0$$ using DSolve: ...
1
vote
1answer
86 views

Solution of differential equation in terms of incomplete gamma function

I need help in solving equation 15 and 16 either manually or in Mathematica to get the solution in terms of the incomplete gamma function. This is what Mathematica tells me. I can't understand ...
0
votes
0answers
72 views

Solving a system of DAE on mathematica

I am having trouble solving a system of Differential Algebraic equations of mathematica, the solution I get is just zeros although it should give me an answer, here is my code: ...
2
votes
1answer
71 views

How to use initial fixed timestep, then decrease it according to dependent variable, while spatial stepsize is fixed

I am trying to solve an advection equation. I want to force constant spatial step size (x dimension) with the “MethodOfLines” option, whereas I want to use initially fixed time step size 0.01 then ...
4
votes
1answer
100 views

Why I can not get the plot when I use NDSolve`ProcessSolutions?

Why I can not get the plot when I use NDSolve`ProcessSolutions? Anyone can give me a clue. Thanks a lot! ...
-2
votes
1answer
80 views

Trouble to assign solution of differential equation to a function

I have written a function to Solve the following differential equation where the number n (called in my code Nphoton) is a variable. Nevertheless I have troubles to assign the solutions ...
9
votes
2answers
716 views

Solve Laplace equation using NDSolve

I am new to Mathematica, a friend recommended this software and started using it, in fact download the trial version to know. I recently did a program in C to calculate numerically the solution to ...
1
vote
1answer
102 views

Save output of NDSolve to a file

I am solving following differential equation using Mathematica, ...
7
votes
1answer
135 views

Simpler code evaluating Dickman's function?

I've read here that Mathematica 10 can Obtain symbolic solutions to delay differential equations. Would that help in numericaly evaluating Dickman's function $\rho(u)$ ? It is a delay differential ...
3
votes
1answer
164 views

reconstruct a 3D curve from discrete curvature and torsion

I tried to reconstruct a 3D curve with given curvature and torsion. I saw some threads talking about using runge kutta. However, as far as I see that they required curvature and torsion were ...
20
votes
1answer
465 views

Calculating a potential function using the finite element method

This is my first attempt to use the Finite Element method available in version 10. There are questions and I am very open to suggestions. My example is flow around a cylinder which is a well known ...
0
votes
0answers
107 views

How can one use differential boundary conditions with helmholzSolve?

I can not get the helmholzSolve function provided by Mark McClure and user21 to work for a case that I want to constrain the spacial derivative of a boundary. In particular, I am looking at the ...
-1
votes
1answer
106 views

Differential Equations: Solving a second order ODE with DSolve

I have the following equation h''[η] + (2 a'/a) h'[η] + (k)^2 h[η] == 0 2 a'/a = 1.3551 with boundary conditions ...
1
vote
3answers
129 views

When NSolve fails due to a differential situation?

c = 1.1111; y[x_] = x - c Sin[x] NSolve[y[x] == 0, x] The method has procured no result. Successive derivatives were plotted in an attempt to fix the problem. ...
0
votes
1answer
62 views

DSolve not returning “trivial” solutions [duplicate]

When I enter this DSolve[y'[x]^2 + y[x]^2 == 1, y[x], x] the answer I get is ...
1
vote
1answer
111 views

NDSolve fails for certain choices of parameters and solve range

I'm trying to solve a pair of coupled ODEs with NDSolve. I know roughly what the solution should look like (both should give periodic functions, pi/2 out of phase, the amplitude of which damp towards ...
2
votes
1answer
70 views

Numerically solving systems of (first order) linear delay differential equations

I was wondering, how to solve such a system with mathematica? I found the NDsolve function, but in the reference is not mentioned if it is possible to solve a system of equations.
4
votes
2answers
125 views

NDSolve not returning the expected solution

I'm trying to simulate a simple circuit with Mathematica. The equation of the circuit is $R \dfrac{dQ}{dt} + \dfrac{Q(t)}{C} = f_{sig}(t)$. This is the definition of $f_{sign}$, and the function ...
1
vote
1answer
80 views

Numerically integrating solution obtained from NDSolve method

In the following example, $u(x)$ is found numerically using NDSolve method. ...
8
votes
2answers
372 views

Driven simple harmonic oscillator — amplitude of steady state motion

I am (partly as an exercise to understand Mathematica) trying to model the response of a damped simple harmonic oscillator to a sinusoidal driving force. I can solve the differential equation with ...
4
votes
1answer
253 views

PDE with Stefan Conditions, a.k.a variable boundary

I want to solve the one-dimensional one-phase Stefan problem, but I don't know how to make Mathematica understand the conditions. If you are not familiar with what I'm asking please refer to this ...
1
vote
1answer
108 views

Real and Imaginary parts of solutions to a complex linear ODE system

Consider a complex linear ODE system $x'=Ax$, where $$A=\left( \begin{matrix} 0&1\\ -2&-i \end{matrix}\right). $$ One can first find the eigenvalues and eigenvectors using ...