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

learn more… | top users | synonyms

0
votes
1answer
88 views

Misunderstanding of Function Notation

I have: sol = DSolve[{y'[t] == y[t]^2 - y[t], y[0] == c}, y[t], t]; y[t_, c_] := y[t] /. First[sol]; Then: y[t,c] Produces ...
5
votes
3answers
253 views

Numerically Simulating a Simple Differential Equation: RC Circuit

I need to learn to simulate a differential equation in Mathematica. The differential equation that I need to simulate is complicated, without an analytical solution. However, I am brand new to ...
2
votes
1answer
70 views

Breaking out of NDSolve

I am solving a coupled set of differential equations with NDSolve for 6 unknown functions of time. At a certain point in time, the system hits a singular point where the potential governing the ...
2
votes
0answers
121 views

Volterra integral equation

I have to find an approximate numerical solution for the equation $$ F(x) - \lambda \int\limits_1^{x} \text{d}s \;s^2 F(s) Z(x-s) = G(x) $$ $$Z(s) = (\psi''(1-2\ h\ i\ s)- 0.5 \psi''(1-2\ h\ i\ s))$$ ...
1
vote
0answers
111 views

Plotting phase portrait from NDSolve (i.e., must plot x'[x], where $x=x(t)$)

I'm doing a double pendulum system and I wish to plot the phase portraits for each angle. I've successfully used NDSolve to solve for each of the angles. Let's just ...
4
votes
1answer
121 views

Can anyone help me to get the solution of this first order ODE?

$\ y'(x)=-\sqrt(y(x)^3)$ with initial condition $\ y(0)=-2$. In fact, this is an example here in page 3-4 ...
1
vote
0answers
69 views

How to solve Fuzzy Differential Equations

How to solve Fuzzy Differential Equations using DSolve or NDSolve, Or Differential Equations with fuzzy number as initial conditions
1
vote
1answer
64 views

Plotting ODE solutions with two parameters

I have a system of two first-order ODEs with two parameters, a and b and I'm trying to plot the solution curve in the X-Y plane at a given initial condition for various parameter values. Can someone ...
0
votes
0answers
71 views

Solving PDE with conditions

I am new to Mathematica, so it maybe a silly question. So, I'm trying to solve the following PDE with Mathematica: $u=u(x,y)$ $\partial_x ^2 u-6\partial_x \partial_yu+9\partial_y^2u=x^2+y^2$ ...
1
vote
0answers
99 views

Why is my output equal to my input?

I tried to solve a system of differential equations with the code below: ...
3
votes
2answers
147 views

Turning arrowheads off

My thinking was that Arrowheads[0] would turn off the arrowheads. ...
8
votes
1answer
276 views

How can I improve speed of code that generates an orbit plot?

I am attempting to generate an orbit for the "time 2 Pi map" of a forced damped pendulum. The example comes from page 57 in the Chaos text by Alligood,Sauer and Yorke, see Figure 2.7. The code ...
1
vote
0answers
66 views

Direction Field for an Ordinary Differential Equation [duplicate]

I have found links for plot the direction field of a differential equation. For example: How can I plot the direction field for a differential equation? And: Plotting the direction field of a ...
5
votes
2answers
3k views

How to handle NDSolve::ndsz problem (singularity problem)

I have 2 second order differential equations (non-linear). The physics behind them is correct. I verified the equations many times. It is a solid pendulum with a mass-spring at the end of it. Now, ...
4
votes
2answers
180 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
0answers
44 views

The use of Evaluate in NDSolve and Plotting the answer [duplicate]

I am having a bit of hard time understanding the use of Evaluate when plotting the output of NDSolve. From the documentation, I enter: ...
15
votes
1answer
1k 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 ...
1
vote
1answer
404 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 ...
2
votes
3answers
461 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
161 views

NDSolve and strange “nonlinear coefficients problem”

I'm stuck solving the following problem. I defined two functions as follows: $$ \varphi(\lambda) = \frac{\left((\lambda-2)^2-1 \right)^2}{4}$$ $$ \gamma(\lambda) = \varepsilon^2 ...
3
votes
1answer
79 views

System of differential equations on a directed graph

I want to use Mathematica to model blood flow through a capillary network. The approach we are using is going to be to create a directed graph, and have a system of differential equations at each ...
1
vote
0answers
73 views

Maxwell propagator problem 2D [closed]

Hi folks I`m trying to turn certain matlab code to mathematica but for some reason it doesnt really want to do it could you please assist me a bit. Thanks The matlab code is as follows: ...
2
votes
2answers
116 views

Can't use my own discontinuous coefficient within NDSolve

I want to evaluate the pendulum's motion with changing equilibrium state given by eqilValue[t]: ...
3
votes
1answer
113 views

Troubles with FindFit for ODE

I need to use FindFit for ODE with arbitary number of parameters. This code (from documentation) works: ...
0
votes
1answer
38 views

Coefficients of a DiffEq [closed]

I'm attempting a Do[] recursion on Coefficient[] with the end goal to use it to find the coefficients of an nth order diffeq. The following is a basic example of what I'm trying to do, subbing a ...
0
votes
1answer
75 views

Plotting and Analysing derivatives of solution of differential equation [closed]

Using NDSolve I found the solution for the differential equation that governs the dynamics of a simple pendulum. How do I plot the derivatives of the solution? Also, I would like to find the ...
0
votes
1answer
785 views

how to solve second order nonlinear coupled differential equations using NDSolve with hyperbolic function

i have to solve some solitons scattering through this coupled equations. i need to get two different graph, but still the graph did not come out. and also the equations quite complicated containing ...
0
votes
0answers
34 views

Finding a specific number of eigenvalues in Mathematica

I'm trying to learn how to determine a specific number for eigenvalues for a transcendental equation. For example, ...
1
vote
1answer
107 views

using Module in NDSolve

I am trying to solve a set of two wave equations where solution of first equation defines initial conditions for the second one. And I would like to use Module to define those equations as a function ...
1
vote
2answers
53 views

Strange behavior of the Evaluate function

I am solving a system of differential equation whose unknown are saved under variable vec={C00,C01,C10,C11} and the solutions are saved under variable ...
1
vote
2answers
172 views
6
votes
1answer
269 views

Poisson PDE using FDM in Mathematica

The problem modeled by the Poisson PDE is related to the torsion of prismatic beams and I use the Finite Differences Method. When I use Mathematica to solve it using a quadratic mesh it is working ...
0
votes
0answers
72 views

Help deciphering Neumann Value

Version 10 has an option for NDSolve called "NeumannValue" which is supposed to allow you to specify boundary conditions of the form $\vec{n} . (c \nabla u + \alpha u -\gamma) = g-q u$ but there is no ...
2
votes
0answers
87 views

Hermite method in Mathematica

When I solve the example below, in the results it is mentioned that MMA 10 has used Hermite method. I cannot find anything about this in MMA documentation. Is there nothing about Hermite method ...
2
votes
2answers
216 views

ParametricPlot3D line not smooth

Multivariate parameters appear to present a jagged appearance of integrands (using default Runge-Kutta ODE integration intervals?) in ParametricPlot3D plotting on a single argument. Higher Mesh.. ...
1
vote
1answer
97 views

Problem solving a differential equation

I want to solve the following differential equation: y''[x]+a*y[x]=-b*u[x], y[0]==0, y'[0]==0 for y(x). u[x] is a vector, ...
2
votes
1answer
87 views

Evaluate function defined by DifferentialRoot

I have the following sequence of rationals that I want to find the generating function of: ...
1
vote
1answer
68 views

Series expansion of InterpolatingFunction obtained from NDSolve

I am trying to obtain a series expansion of the numerical solution of a differential equation. I encounter difficulties going beyond first-order expansions which I believe might be due to my inability ...
3
votes
1answer
251 views

Combining Gravity Turn and Orbit Models

I have a mathematical model for the motion of an orbiting spacecraft about Earth: ...
0
votes
0answers
72 views

DSolve echoes input

I'm attempting to solve a set of coupled, non-linear ordinary differential equations using DSolve, but Mathematica simply echoes the code that I put in: ...
1
vote
0answers
66 views

Updating procedural fitting algorithm to more efficient style?

For some time now, I have used this procedural programming approach to fit data in my research that is too cumbersome for the built in Mathematica functions. (Although this simple example works fine ...
0
votes
2answers
66 views

Derivative in function form

I have a function in mathematica defined as b[n_, x_] := x^n/n!/\!\(\*UnderoverscriptBox[\(\[Sum]\), \(k = n\), \(\[Infinity]\)]\*FractionBox[\(x^k\), \(k!\)]\) ...
0
votes
0answers
261 views

Problem with Plot and ParametricNDSolve

This is my first question on this site, I hope it is written in an understandable way. Let's start. My aim is to plot a certain parametric region, with $x$ and $y$ coordinates given by the following ...
5
votes
2answers
101 views

How to add (energy) constraint when using NDSolve to Equation of Motion

To simplify my problem, I will try and solve the Equation of Motion for a particle in a 1D Harmonic Potential. energy[x_, p_, m_, ω_] := p^2/(2 m) + (m ω^2)/2 x^2 ...
0
votes
0answers
26 views

Identification of a parameter in a first order differential equation [duplicate]

I wan't to identify a parameter in a first oder differential equation. My goal is proof of the calculation principle i use. First I calculate a response where the paramater is set to 0.3. Later this ...
0
votes
0answers
68 views

Solving an ODE, where the coefficients are implicit functions of time, not in closed form [duplicate]

I have an ODE, say of the first order for simplicity, and of the form, $A.x'[t] =B$. The coefficients are functions of time and x[t], not in closed form. I can define the coefficients, $A$ and $B$, ...
1
vote
1answer
47 views

NDSolve with WhenEvent, resetting system for a prolonged period

I am currently working on a complex system where I would like to (for the lack of a better description) reset (part) of the system for a certain period. Say I five ODE's: ...
1
vote
0answers
54 views

A Stupid Question Maybe - EigenNDSolve

I would like to find out how accurate and how useful the chebyshev (collocation?) method is in finding many eigenvalues to a second-order ODE in one go. Specifically, I used the Mathematica package ...
2
votes
1answer
132 views

Analytic approximation of NDSOLVE

I have to solve the following differential equation L'[t] == L[t]^(-3/2) - 0.1 I tried DSolve: ...
0
votes
1answer
112 views

Tough Calculation, novice mathematica user

I have an equation, that I've been calling $b_N(x)$ that satisfies the following identity: $$-Nb_N(x)^2=(x-N)b_N(x)+xb_N'(x)$$ where $b_N'(x)$ is the first derivative. I take the derivative then of ...