dx/dt=ax+sinx => find the general solution of this eqn.

Results 1 to 9 of 9

- Nov 29th 2008, 09:48 AM #1

- Nov 29th 2008, 02:29 PM #2

- Joined
- Aug 2008
- Posts
- 903

- Nov 29th 2008, 09:54 PM #3

- Nov 29th 2008, 10:50 PM #4
Quite right. It should be

.

If you have reason for thinking otherwise, please don't be backward in coming forward.

By the way, if this is meant to be some sort of challenge question, please say so clearly in the question itself. Members are happy to help people who are genuinely stuck and need help. However, some members have better things to do with their time if that's not the case.

- Nov 30th 2008, 01:07 AM #5

- Nov 30th 2008, 02:30 AM #6

- Nov 30th 2008, 04:28 AM #7

- Nov 30th 2008, 05:05 AM #8

- Joined
- Aug 2008
- Posts
- 903

Hey mathemanyak. I do not feel that's an answer for one thing where's the t? Not sure what you did but allow me to make a suggestion which I have found very useful in helping me with differential equations: Solve it first numerically, then superimpose any analytical expression I get over the numerical solution. If they don't exactly agree within some tolerance, then it's very likely the analytical expression is wrong baring singularities or high oscillating problems or other non-linear effects. So I numerically solved the IVP:

using Mathematica's numerical solver "NDSolve":

Code:sol = NDSolve[{Derivative[1][x][t] == x[t] + Sin[x[t]], x[0] == 1}, x, {t, 0, 5}] Plot[Evaluate[x[t] /. sol], {t, 0, 5}]

- Nov 30th 2008, 09:43 AM #9

- Joined
- Aug 2008
- Posts
- 903

I was curious about this so I constructed a 10'th degree polynomial fit to the inverse data. That is, I numerically integrated the integral from 1 to 200 (increment by 1) and then formed a table of (t,x) pairs: and fit it to the polynomial. The plot below shows a overlay of the NDSolve results (in red) with the polynomial fit (green). I think red+green=brown.

Code:f[x_] := NIntegrate[1/(u + Sin[u]), {u, 1, x}]; tlist = Table[{f[x], x}, {x, 1, 200, 1}]; lp = ListPlot[tlist] clist = Table[Subscript[a, n], {n, 0, 10}] p[x_] = Sum[Subscript[a, n]*x^n, {n, 0, 10}]; alist = FindFit[tlist, p[x], clist, x]; pic1 = Plot[p[x] /. alist, {x, 0, 5}, PlotStyle -> Red] sol = NDSolve[{Derivative[1][y][t] == y[t] + Sin[y[t]], y[0] == 1}, y, {t, 0, 5}]; pic2 = Plot[Evaluate[y[t] /. sol], {t, 0, 5}, PlotStyle -> Green] Show[{pic2, pic1}, PlotRange -> {{0, 5}, {0, 200}}]