Differential equations: y’= sin(cx) sin(yx), y(0)=1

Faema
Messages
2
Reaction score
0
It’s a non-linear first order ode, the trig form of which is y’= sin(cx) sin(yx), y(0)=1. It’s possible that I’ve forgotten an important detail in the construction of the equation. If there’s someone who has some ideas about solving or approximating this, or who knows how to get mathematica or MATLAB to do something useful with it, I would be so grateful. Thank you in advance!
 
Physics news on Phys.org
Here's one way. Since y(0)= 1 and y'(x)= sin(cx)sin(xy), y'(0)= sin(0)sin(0)= 0.
Differentiating, y''(x)= ccos(cx)sin(xy)+ ysin(cx)cos(xy) and so y''(0)= 0. Differentiating again, y'''(x)= -c^2sin(cx)sin(xy)+ 2cycos(cx)cox(xy)- y^2sin(cx)sin(xy) and y'''(0)= 2cy(0)= 2c. So far, we have, for the Taylor's series about x=0, y(x)= 1+ 0x+ 0x^2+ 2c/3! x^3= 1+ (c/3)x^3. Although, obviously, the calculations are going to get more complicated, continuing in that way, you will have a power series approximatiom for y.
 
Faema said:
It’s a non-linear first order ode, the trig form of which is y’= sin(cx) sin(yx), y(0)=1. It’s possible that I’ve forgotten an important detail in the construction of the equation. If there’s someone who has some ideas about solving or approximating this, or who knows how to get mathematica or MATLAB to do something useful with it, I would be so grateful. Thank you in advance!

This is what I do: give it a few minutes of a little of this, little of that myself. Ok, no go. Then I shamelessly turn to Mathematica's DSolve which attempts to find a symbolic solution. Still no go. I'm feelin' good by now cus' Mathematica can't do it neither. Finally I stop trying to be a purist and just turn to Mathematica's NDSolve which computes a numeric solution of the IVP:

myc=1;
mysol=NDSolve[{y'[x]==Sin[myc x] Sin[y[x] x],y[0]=0},y,{x,0,1}];
Plot[y[x]/.mysol,{x,0,1}]

. . . yeah, I know. I'm not proud. Hall's got a good idea up there. If you're a Mathematica junkie, it wouldn't be hard to get maybe say 20-30 of those derivatives, maybe more, before Mathematica starts choking. Then form the series and compare it to the numerical results. Be good practice for when you really need to do something with the software.
 
Last edited:
This looks promising! Thankyou! Looking into a Taylor expansion..
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top