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

learn more… | top users | synonyms

0
votes
0answers
38 views

Evaluate function gives zero while plot shows non-zero function

I have obtained two functions, f and g, as two of the solutions to a certain set of PDE's (using NDSolve). I cannot provide the differential equations here because the system is very large, but will ...
1
vote
2answers
98 views

Mathematica won't plot the a function that is solved from NDSolve

I have used NDSolve to solve for a function (an angle of a triple pendulum), and now I wish to plot the derivative of that function (the angular velocity). ...
0
votes
1answer
64 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
4answers
151 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
45 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
2answers
120 views

NDSolve solution does not satisfy boundary conditions

I am trying to solve a 1st order non-linear ODE ...
1
vote
0answers
87 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 ...
1
vote
0answers
50 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
48 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
46 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
80 views

Why is my output equal to my input?

I tried to solve a system of differential equations with the code below: ...
1
vote
0answers
102 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))$$ ...
3
votes
2answers
127 views

Turning arrowheads off

My thinking was that Arrowheads[0] would turn off the arrowheads. ...
4
votes
1answer
108 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
41 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 ...
7
votes
1answer
231 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
41 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: ...
0
votes
0answers
67 views

Having difficulty with NDSolve [closed]

I've entered: NDSolve[{y'[x] == z[x], z'[x] == -y[x], y[0] == 0, z[0] == 1}, {y, z}, {x, 0, Pi}] The code above was copied and pasted from the Mathematic help ...
3
votes
1answer
66 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
56 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
0answers
90 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 ...
2
votes
2answers
92 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]: ...
0
votes
1answer
37 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 ...
3
votes
1answer
83 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
63 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
0answers
33 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
2answers
48 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
1answer
74 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 ...
0
votes
0answers
40 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 ...
0
votes
1answer
119 views

Help solving the following non-linear differential equation

After spending several hours exploring different options, I still couldn't not figure out how to solve the following differential equations. Any suggestions would be highly appreciated. ...
2
votes
0answers
78 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 ...
1
vote
1answer
88 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
2answers
163 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.. ...
2
votes
1answer
59 views

Evaluate function defined by DifferentialRoot

I have the following sequence of rationals that I want to find the generating function of: ...
0
votes
0answers
35 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
2answers
107 views

ParametricNDSolve and plotting problem

I have a problem with ParametricNDSolve: ...
2
votes
1answer
132 views

Combining Gravity Turn and Orbit Models

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

solving the heat equation

I'm trying to model the diffusion equation on mathematica, but I'm a complete beginner. So far I've managed to find the heat equation on the NDSovle examples but I'm not really entirely sure what all ...
0
votes
2answers
52 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!\)]\) ...
1
vote
0answers
58 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
0answers
25 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 ...
1
vote
1answer
34 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
44 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 ...
4
votes
2answers
87 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
1answer
111 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 ...
0
votes
0answers
34 views

NDSolve error message: Cannot find starting value for the variable y'

I'm attempting to solve a differential equation using NDSolve. This is the code: ...
2
votes
1answer
129 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
81 views

StreamPlot in a system x´ = f[x] + g[t] [closed]

StreamPlot can be used in the O.D.E system: x' = 2 x - 3 y + 3 t y' = 5 x + y - t How? ...
6
votes
1answer
216 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 ...
1
vote
1answer
50 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 ...