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

learn more… | top users | synonyms (3)

2
votes
1answer
125 views

solving Stokes equations for Force=0 boundary condition

Say you want to solve the movement of a particle in a viscous fluid such that the force on the particle is zero. In general, the solution is trivial since one has no-slip boundary conditions, ie the ...
0
votes
1answer
54 views

Using DSolve for r(t) y'(t) = y(t) r'(t) [closed]

I am fairly new to Mathematica and I was trying to solve this equation r(t) y'(t) = y(t) r'(t) using DSolve. How can I adapt it to the DSolve function formalism?
0
votes
1answer
79 views

DEigensystem and NDEigensystem simply not working in Mathematica 10.0 [closed]

I am trying to use DEigensystem and NDEigensystem to check some theoretical derivations I have recently done. I am running ...
1
vote
2answers
120 views

How to plot coupled Differential equations with changing parameters?

I'm a new user of mathematica. What I'm trying to do is following: I have two coupled non-linear differential equations. Say $x'(t)$ = $A− (B +1)x + x^2 y$ $y'(t) = Bx − x^2 y$ What I want is to ...
0
votes
0answers
46 views

Overdetermined Differential Equation is not overdetermined

I want to solve a system of two differential equations ...
0
votes
1answer
109 views

Animate spherical pendulum

I want to plot the position of Spherical pendulum. First I tried to plot the simple pendulum: ...
1
vote
1answer
68 views

Numerical integration plot

I am trying to do a numerical integration plot in Mathematica. Specifically, I want to plot $r$ versus $\tau$ where, $$\frac{dr}{d\tau}=\sqrt{2F(r)}$$ $F$ being a quartic polynomial. The limit of $r$ ...
3
votes
1answer
113 views

how to speed up FEM-NDSolve for similar boundary conditions?

I need to solve several stoke's equations for very similar boundary conditions, which represent particles in a fluid with different velocities U[i]. I'm interested in computing the force acting on the ...
2
votes
1answer
72 views

Heat transportation equation 2

I would like to thank you for all the suggestions for the problem with the heat transportation equation I have posted time ago! I rewrited the code along with the implementation of the Neumann ...
3
votes
2answers
80 views

Animating 3D Parametric Plot

I'm trying to create an animation of Lorenz equation as a function of t. Here is the numerical solution of the equation: ...
1
vote
1answer
80 views

Problem with Finding Minimum of a Function Solved With NDSolve

I am trying to find the minimum of a particular function, which is solved using NDSolve. What is making it a little more complicated is that it is not the only function solved within the NDSolve and ...
0
votes
0answers
28 views

Use a Poincare map and Poincare section to show that the equation has a T-periodic solution [duplicate]

The equation is $x'(t)=-x(t)^3+h(t)$, where $h(t)$ is any smooth T-periodic function. Give an argument using a suitably constructed Poincare section/map that shows the equation has a T-periodic ...
3
votes
3answers
159 views

Can I plot the solution to an ODE without using dsolve?

I have the following ODE $\qquad x'= -x^3+\sin(t)$ I want to find numerically the initial condition $x(0)$ corresponding to the $2\pi$ periodic solution. I tried to use ...
1
vote
0answers
48 views

NDSolve cannot consistently handle a simple 2D reaction-diffusion PDE with spiky (and random) initial conditions

I am still fighting with a modified version of the 2D reaction-diffusion PDE that I have already posted in a previous question ( Reaction-diffusion PDE with NDSolve: either very slow or very ...
1
vote
1answer
113 views
0
votes
1answer
38 views

DSolve “Equation or list of equations expected instead of True in the first \ argument ”

I am getting the error: Equation or list of equations expected instead of True in the first argument when I try to evaluate the cell. ...
1
vote
1answer
46 views

Polar and Parametric Plots not showing up

I bet this has asked before, but I can't seem to find an answer anywhere. I have x and y functions which, individually plotted, are showing up fine, but when I try to plot the parametrically, the ...
6
votes
2answers
152 views

NDSolve and matrix formulation do not produce the same result. Why and which is more precise?

Let's assume a vector of variables $\vec a=(a_1(t),a_2(t),a_3(t))$ and let's assume a following DE$$C\ddot{\vec a}=B\dot{\vec a}+A\vec{a}+\vec F.$$ If $A$ and $C$ are both reversible, than the DE can ...
20
votes
2answers
363 views

Numerical optimal control

I was hoping to tackle optimal control using Mathematica in order to learn how I can use Mathematica's built in numerical integration and optimization functions together in order to solve an optimal ...
-2
votes
1answer
54 views

Nonlinear system of ODEs need phase-plane

I have this system of equations : I need to solve this equation for different parameters, for now theta is 1, alpha is 25 and ...
0
votes
1answer
44 views

Failed to extract parameters by fitting data with ODE

I'm trying to extract 4 parameters by fitting one set of data with an ODE via NonlinearModelFit. I originally thought the process should be straightforward, but later on I found I was totally wrong. ...
0
votes
1answer
52 views

NDSolve Unable to continue with complex values or beyond floating-point exceptions error

Initially tried to solve following first degree of differential equations using following method ...
0
votes
0answers
42 views

Solving system of ordinary differential equations where coefficients are matrices

I am attempting to use NDSolve to solve a system of ordinary differential equations defined by matrices. Here's the setup of my problem and what I have tried so far. I'm interested in applying the ...
2
votes
0answers
47 views

Unusual error message from Solve when called by DSolve

Question 111150 presents an equation nonlinear in its leading derivative ...
3
votes
2answers
217 views

Is there is a decreasing function that fits my data?

Cross-posted in http://tieba.baidu.com/p/4440556337 http://tieba.baidu.com/p/4440547130 I was wondering if there is a function that fits the following data perfectly or properly. There is a ...
0
votes
1answer
63 views

Sum the several solution functions returned by DSolve [duplicate]

Could you please explain why this code does not work. ...
0
votes
0answers
38 views

fit multiple dataset with a single ODE to extract parameters

I found my problem very similar to the one discussed in here. But I'm having difficulties to fit multiple dataset (with varied temperatures) to a single ODE. The hurdle for me now is since the rate ...
8
votes
1answer
121 views

Beam deformation due to localized force distribution

I would like to simulate the deformation of a beam which is fixed on one end and has some localized forces applied along the top. I have been inspired by this Mathematica example: ...
1
vote
0answers
65 views

Solving two PDEs (for f[x,y] and g[x,y]) coupled through a non-linear function |Grad[f[x,y]]|^2

I am new to Mathematica, trying to solve two PDEs (for f[x,y] and g[x,y]) coupled through Grad[f[x,y]]. In the following code when n is set to 1, NDSolve seems to be able to solve for f[x,y] and ...
2
votes
1answer
72 views

Problem with 3rd order ODE solution from DSolve

I Dsolved the following 3rd order ODE: x'[t]==(3x''[t]^2-Sqrt[9x''[t]^4+4x'''[t]^2])/(2x'''[t]) Mathematica (version 8) gave me three solutions including the ...
1
vote
1answer
57 views

How to rescale a NDSolve curve to fit in a box?

I'm numerically solving some differential equations and the 3D curve I get has a very complicated shape and size that cannot be predicted from the initial conditions. I need to rescale the whole ...
0
votes
1answer
72 views

Solve the n ODE in the loop, add the solutions to create one function

I need to solve the n ODE and create a function which is a sum of all the solutions. I wanted to use the code below, but only part (solve) and (create solution) works well. How can I store the ...
2
votes
1answer
125 views

The numerical solution of a nonlinear ODE: boundary value problem

I want to solve the following nonlinear ordinary differential equation with boundary values $f(0)=1$ and $f(1)=1$: $${\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}f \left( x \right) ={\dfrac { \left( 51-51\,x ...
2
votes
2answers
72 views
1
vote
2answers
83 views

A way to show all the 3D curve, with view still centered on {0, 0, 0}

I frequently have this graphics output problem : I draw a 3D curve using NDSolve, and depending on the initial data, there's no way I could predict the shape and ...
4
votes
3answers
649 views

Does it mean that there no analytical solution if Mathematica cannot find analytical solution?

I'm solving simple but coupled ODEs recently. I use both MATLAB symbolic computation and Mathematica. For example, my coupled ODE is the following \begin{align*} &\dot{x}(t)=y(t)-\rho ...
1
vote
1answer
70 views

Variable coefficient Differential Equation

Can anyone help with the following differential equation. The errors i am getting are also given. ...
1
vote
1answer
92 views

Solve a system of nonlinear ODE

I have a set of coupled non-linear ODEs, they look symmetric. But when I use matlab symbolic calculation to solve, matlab told me that no explicit solution found. I was recommended by my friend to try ...
0
votes
1answer
58 views

NDSolve and a vector function [duplicate]

I'm trying to solve a set of differential equations which include a numerical function that takes a vector and returns a vector. ...
2
votes
1answer
117 views

Improving the speed of evaluating numerical solution to a sytem of PDE

I am running a simulation for a kinetic reaction inside a pellet modeled by a system fo reaction-diffusion equations. ...
1
vote
1answer
55 views

How to define “DiffusionCoefficients” in InitializePDECoefficients that depend on ElementMarker?

Modelling a Poisson Equation on piecewise constant material, ElementMarkers on the PointElements are used for the boundary ...
3
votes
2answers
248 views

Why can't FindFit get the proper result in this problem?

Cross-posted in http://tieba.baidu.com/p/4425634583 The mathematica helps me solve one differential equation in numerial method and got the data of the answer,but the taylor series fit is not good in ...
0
votes
1answer
52 views

Reduce in NDSolve [duplicate]

Second order ode; neither the trig RHS simplifies nor is the solution obtained readily: ...
1
vote
0answers
50 views

How to use interpolating functions in NSolve with a piecewise function?

I have a system of two non-linear differential equations that are piecewise, the solutions of which I would like to use to obtain the values of a third variable that is also piecewise. I have trouble ...
1
vote
1answer
114 views

Problem with solving a system of differential equations [closed]

I am trying to reproduce a result that is part of a derivation of the flow due to a rotating disk. I have as given this system of differential equations. (In these equations primes indicate ...
0
votes
0answers
47 views

ndsz : step size is effectively zero; singularity or stiff system suspected

This is the first time I ask a question. I have seen many solutions to this error and tried but they are not working. Here is the code. Needs["DifferentialEquations...
0
votes
0answers
46 views

Output of NDSolve is useless when the system gets bigger

I have a bridge, which I want to discretize on $n$ stiff rods. To find a solution I wrote the following code in mathematica: ...
3
votes
0answers
53 views

EvaluationMonitor slows down NDSolve considerably

Hi guys I'm solving a nonlinear system of ODEs with NDSolve. Since the dimension of the system ranges from 2700 to 20000 equations, depending on the parameters, the time required to solve the entire ...
-3
votes
1answer
78 views

How to implement Bifurcation of delayed Lorenz system? [closed]

The delayed Lorenz system is as follows: (Bifurcation parameter as $\tau$) $$\begin{align*} D^\alpha x_1(t) = & a_1(x_2(t-\tau)-x_1(t)) \\ D^\alpha x_2(t) = & a_2x_1(t)-x_2(t)-x_1(t)x_3(t) ...