Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

How to solve the ordinary differential equation $\frac{d^2y}{dx^2}+\sin(x+y)=\sin x,y(0)=0,y'(0)=1$

Then its possible to solve it by numerical methods?

share|cite|improve this question
I doubt an exact solution exists. – David H Aug 17 '14 at 23:56
Given that it's a highly non-linear ODE, the dependence on boundary conditions / initial values is quite sensitive. – Semiclassical Aug 18 '14 at 0:17
@Semiclassical Why do you consider it highly nonlinear? The nonlinear part does not involve derivatives of unknown function, and is uniformly Lipschitz. This is as close to linear as one might hope. – Bookend Aug 18 '14 at 0:26
up vote 3 down vote accepted

Yes, you can easily solve this using available numerical solvers. As an example, I will use Sage since it's available online and is free. First step is to convert to first-order linear system by introducing another unknown $u=y'$: $$ y'=u,\quad u'=\sin x-\sin(x+y) $$ The rest is up to Sage:

var('x y u')
list_plot([ [i,j] for i,j,k in P])

Pretty neat plot:


It does some other weird things later on, as you can see by increasing the end_points parameter.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.