Unsolvable first order nonlinear differential equation? $frac{dx(t)}{dt} + sin(x(t)) = sin(omega t)$
I do want to solve the following differential equation analytically:
$$frac{dx(t)}{dt} + sin(x(t)) = sin(omega t)$$
I tried several methods to solve this equation, unfortunately without any success. In the last week, I have read a lot of papers related to that kind of prototype and have looked into all the books that deal with first order differential equations. It still seems to be impossible. Neither Wolfram-alpha, nor Matlab's symbolic toolbox can give me a solution for it. Because of that, I am wondering if there is actually a solution for that kind of differential equation?
Thanks for your help!
calculus differential-equations
asked Nov 20 at 17:09
R. Caloudis
444
|
show 3 more comments
I do want to solve the following differential equation analytically:
$$frac{dx(t)}{dt} + sin(x(t)) = sin(omega t)$$
I tried several methods to solve this equation, unfortunately without any success. In the last week, I have read a lot of papers related to that kind of prototype and have looked into all the books that deal with first order differential equations. It still seems to be impossible. Neither Wolfram-alpha, nor Matlab's symbolic toolbox can give me a solution for it. Because of that, I am wondering if there is actually a solution for that kind of differential equation?
Thanks for your help!
calculus differential-equations
asked Nov 20 at 17:09
R. Caloudis
444
2
No CAS being able to solve is a very good sign that you can't neither. Due to the non-linearity (sine), I doubt that there is a closed-form expression. For small angles, use $sin xapprox x$.
– Yves Daoust
Nov 20 at 17:21
What about assumptions? Could that help you to find a solution? Actually I can not believe that such an equation is not solvable.
– R. Caloudis
Nov 20 at 17:24
3
It's very likely that there are no closed-form solutions to this differential equation. Most differential equations are like that. Of course there are solutions. You can solve the differential equation numerically, or find arbitrarily many terms of a series.
– Robert Israel
Nov 20 at 17:24
Thanks. However, I cannot use the small-angle approximation, because I do need that term in order to describe a certain physical phenomena. Any more guesses?
– R. Caloudis
Nov 20 at 17:26
The problem is the following: You do not have saturation nor harmonics without the sinus term. If you excite with a sinusoid at high amplitude the response would contain only the fundamental frequency without the term.
– R. Caloudis
Nov 20 at 17:30
|
show 3 more comments
I do want to solve the following differential equation analytically:
$$frac{dx(t)}{dt} + sin(x(t)) = sin(omega t)$$
I tried several methods to solve this equation, unfortunately without any success. In the last week, I have read a lot of papers related to that kind of prototype and have looked into all the books that deal with first order differential equations. It still seems to be impossible. Neither Wolfram-alpha, nor Matlab's symbolic toolbox can give me a solution for it. Because of that, I am wondering if there is actually a solution for that kind of differential equation?
Thanks for your help!
calculus differential-equations
asked Nov 20 at 17:09
R. Caloudis
444
I do want to solve the following differential equation analytically:
$$frac{dx(t)}{dt} + sin(x(t)) = sin(omega t)$$
I tried several methods to solve this equation, unfortunately without any success. In the last week, I have read a lot of papers related to that kind of prototype and have looked into all the books that deal with first order differential equations. It still seems to be impossible. Neither Wolfram-alpha, nor Matlab's symbolic toolbox can give me a solution for it. Because of that, I am wondering if there is actually a solution for that kind of differential equation?
Thanks for your help!
calculus differential-equations
calculus differential-equations
asked Nov 20 at 17:09
R. Caloudis
444
asked Nov 20 at 17:09
R. Caloudis
444
edited Nov 20 at 18:10
asked Nov 20 at 17:09
R. Caloudis
444
asked Nov 20 at 17:09
R. Caloudis
444
asked Nov 20 at 17:09
R. Caloudis
444
444
2
No CAS being able to solve is a very good sign that you can't neither. Due to the non-linearity (sine), I doubt that there is a closed-form expression. For small angles, use $sin xapprox x$.
– Yves Daoust
Nov 20 at 17:21
What about assumptions? Could that help you to find a solution? Actually I can not believe that such an equation is not solvable.
– R. Caloudis
Nov 20 at 17:24
3
It's very likely that there are no closed-form solutions to this differential equation. Most differential equations are like that. Of course there are solutions. You can solve the differential equation numerically, or find arbitrarily many terms of a series.
– Robert Israel
Nov 20 at 17:24
Thanks. However, I cannot use the small-angle approximation, because I do need that term in order to describe a certain physical phenomena. Any more guesses?
– R. Caloudis
Nov 20 at 17:26
The problem is the following: You do not have saturation nor harmonics without the sinus term. If you excite with a sinusoid at high amplitude the response would contain only the fundamental frequency without the term.
– R. Caloudis
Nov 20 at 17:30
|
show 3 more comments
2
No CAS being able to solve is a very good sign that you can't neither. Due to the non-linearity (sine), I doubt that there is a closed-form expression. For small angles, use $sin xapprox x$.
– Yves Daoust
Nov 20 at 17:21
What about assumptions? Could that help you to find a solution? Actually I can not believe that such an equation is not solvable.
– R. Caloudis
Nov 20 at 17:24
3
It's very likely that there are no closed-form solutions to this differential equation. Most differential equations are like that. Of course there are solutions. You can solve the differential equation numerically, or find arbitrarily many terms of a series.
– Robert Israel
Nov 20 at 17:24
Thanks. However, I cannot use the small-angle approximation, because I do need that term in order to describe a certain physical phenomena. Any more guesses?
– R. Caloudis
Nov 20 at 17:26
The problem is the following: You do not have saturation nor harmonics without the sinus term. If you excite with a sinusoid at high amplitude the response would contain only the fundamental frequency without the term.
– R. Caloudis
Nov 20 at 17:30
2
2
No CAS being able to solve is a very good sign that you can't neither. Due to the non-linearity (sine), I doubt that there is a closed-form expression. For small angles, use $sin xapprox x$.
– Yves Daoust
Nov 20 at 17:21
No CAS being able to solve is a very good sign that you can't neither. Due to the non-linearity (sine), I doubt that there is a closed-form expression. For small angles, use $sin xapprox x$.
– Yves Daoust
Nov 20 at 17:21
What about assumptions? Could that help you to find a solution? Actually I can not believe that such an equation is not solvable.
– R. Caloudis
Nov 20 at 17:24
What about assumptions? Could that help you to find a solution? Actually I can not believe that such an equation is not solvable.
– R. Caloudis
Nov 20 at 17:24
3
3
It's very likely that there are no closed-form solutions to this differential equation. Most differential equations are like that. Of course there are solutions. You can solve the differential equation numerically, or find arbitrarily many terms of a series.
– Robert Israel
Nov 20 at 17:24
It's very likely that there are no closed-form solutions to this differential equation. Most differential equations are like that. Of course there are solutions. You can solve the differential equation numerically, or find arbitrarily many terms of a series.
– Robert Israel
Nov 20 at 17:24
Thanks. However, I cannot use the small-angle approximation, because I do need that term in order to describe a certain physical phenomena. Any more guesses?
– R. Caloudis
Nov 20 at 17:26
Thanks. However, I cannot use the small-angle approximation, because I do need that term in order to describe a certain physical phenomena. Any more guesses?
– R. Caloudis
Nov 20 at 17:26
The problem is the following: You do not have saturation nor harmonics without the sinus term. If you excite with a sinusoid at high amplitude the response would contain only the fundamental frequency without the term.
– R. Caloudis
Nov 20 at 17:30
The problem is the following: You do not have saturation nor harmonics without the sinus term. If you excite with a sinusoid at high amplitude the response would contain only the fundamental frequency without the term.
– R. Caloudis
Nov 20 at 17:30
|
show 3 more comments
1 Answer
1
active
oldest
votes
Of course solutions do exist, we might just be unable to find closed form representations for them. This doesn't mean we can't compute them numerically: here are some IVP trajectories with $omega=1$
You can plot just a fundamental patch $tin[0,2pi]$, $xin[-pi,pi]$ thanks to the periodicity (I've highlighted in red the stable periodic solution and in blue the unstable periodic solution):
Alternative visualization:
Edit. Code for the plots, in Mathematica:
sol = ParametricNDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x0},
x, {t, 0, 4 [Pi]}, {{x0, -3 [Pi], 3 [Pi]}}];
Plot[Evaluate@Table[
x[x0][t] /. sol, {x0, -3 [Pi], 2.5 [Pi], .25}], {t, 0, 3.5 [Pi]}]
per = NDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x[2 [Pi]]},
x, {t, 0, 2 [Pi]}];
Show[StreamPlot[{1, Sin[t] - Sin[x]}, {t, 0, 2 [Pi]}, {x, -[Pi], [Pi]}],
Plot[Evaluate[x[t] /. per], {t, 0, 2 [Pi]}, PlotStyle -> Red]]
answered Nov 20 at 17:39
Federico
4,649514
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
1
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
1
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3006597%2funsolvable-first-order-nonlinear-differential-equation-fracdxtdt-sin%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
Of course solutions do exist, we might just be unable to find closed form representations for them. This doesn't mean we can't compute them numerically: here are some IVP trajectories with $omega=1$
You can plot just a fundamental patch $tin[0,2pi]$, $xin[-pi,pi]$ thanks to the periodicity (I've highlighted in red the stable periodic solution and in blue the unstable periodic solution):
Alternative visualization:
Edit. Code for the plots, in Mathematica:
sol = ParametricNDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x0},
x, {t, 0, 4 [Pi]}, {{x0, -3 [Pi], 3 [Pi]}}];
Plot[Evaluate@Table[
x[x0][t] /. sol, {x0, -3 [Pi], 2.5 [Pi], .25}], {t, 0, 3.5 [Pi]}]
per = NDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x[2 [Pi]]},
x, {t, 0, 2 [Pi]}];
Show[StreamPlot[{1, Sin[t] - Sin[x]}, {t, 0, 2 [Pi]}, {x, -[Pi], [Pi]}],
Plot[Evaluate[x[t] /. per], {t, 0, 2 [Pi]}, PlotStyle -> Red]]
answered Nov 20 at 17:39
Federico
4,649514
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
1
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
1
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
add a comment |
Of course solutions do exist, we might just be unable to find closed form representations for them. This doesn't mean we can't compute them numerically: here are some IVP trajectories with $omega=1$
You can plot just a fundamental patch $tin[0,2pi]$, $xin[-pi,pi]$ thanks to the periodicity (I've highlighted in red the stable periodic solution and in blue the unstable periodic solution):
Alternative visualization:
Edit. Code for the plots, in Mathematica:
sol = ParametricNDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x0},
x, {t, 0, 4 [Pi]}, {{x0, -3 [Pi], 3 [Pi]}}];
Plot[Evaluate@Table[
x[x0][t] /. sol, {x0, -3 [Pi], 2.5 [Pi], .25}], {t, 0, 3.5 [Pi]}]
per = NDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x[2 [Pi]]},
x, {t, 0, 2 [Pi]}];
Show[StreamPlot[{1, Sin[t] - Sin[x]}, {t, 0, 2 [Pi]}, {x, -[Pi], [Pi]}],
Plot[Evaluate[x[t] /. per], {t, 0, 2 [Pi]}, PlotStyle -> Red]]
answered Nov 20 at 17:39
Federico
4,649514
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
1
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
1
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
add a comment |
Of course solutions do exist, we might just be unable to find closed form representations for them. This doesn't mean we can't compute them numerically: here are some IVP trajectories with $omega=1$
You can plot just a fundamental patch $tin[0,2pi]$, $xin[-pi,pi]$ thanks to the periodicity (I've highlighted in red the stable periodic solution and in blue the unstable periodic solution):
Alternative visualization:
Edit. Code for the plots, in Mathematica:
sol = ParametricNDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x0},
x, {t, 0, 4 [Pi]}, {{x0, -3 [Pi], 3 [Pi]}}];
Plot[Evaluate@Table[
x[x0][t] /. sol, {x0, -3 [Pi], 2.5 [Pi], .25}], {t, 0, 3.5 [Pi]}]
per = NDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x[2 [Pi]]},
x, {t, 0, 2 [Pi]}];
Show[StreamPlot[{1, Sin[t] - Sin[x]}, {t, 0, 2 [Pi]}, {x, -[Pi], [Pi]}],
Plot[Evaluate[x[t] /. per], {t, 0, 2 [Pi]}, PlotStyle -> Red]]
answered Nov 20 at 17:39
Federico
4,649514
Of course solutions do exist, we might just be unable to find closed form representations for them. This doesn't mean we can't compute them numerically: here are some IVP trajectories with $omega=1$
You can plot just a fundamental patch $tin[0,2pi]$, $xin[-pi,pi]$ thanks to the periodicity (I've highlighted in red the stable periodic solution and in blue the unstable periodic solution):
Alternative visualization:
Edit. Code for the plots, in Mathematica:
sol = ParametricNDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x0},
x, {t, 0, 4 [Pi]}, {{x0, -3 [Pi], 3 [Pi]}}];
Plot[Evaluate@Table[
x[x0][t] /. sol, {x0, -3 [Pi], 2.5 [Pi], .25}], {t, 0, 3.5 [Pi]}]
per = NDSolve[{
x'[t] + Sin[x[t]] == Sin[t],
x[0] == x[2 [Pi]]},
x, {t, 0, 2 [Pi]}];
Show[StreamPlot[{1, Sin[t] - Sin[x]}, {t, 0, 2 [Pi]}, {x, -[Pi], [Pi]}],
Plot[Evaluate[x[t] /. per], {t, 0, 2 [Pi]}, PlotStyle -> Red]]
answered Nov 20 at 17:39
Federico
4,649514
edited Nov 21 at 15:33
answered Nov 20 at 17:39
Federico
4,649514
answered Nov 20 at 17:39
Federico
4,649514
answered Nov 20 at 17:39
Federico
4,649514
4,649514
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
1
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
1
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
add a comment |
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
1
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
1
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
Thanks! What exactly does the fundamental patch tell me? It does not seem to be just another representation of the solution x(t) with respect to t.
– R. Caloudis
Nov 20 at 18:01
1
1
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
Well, you know that the entire plane is filled with this repeating pattern, so this picture gives a pretty complete idea of the behavior of the solutions.
– Federico
Nov 20 at 18:03
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
I've now added a plot of the periodic solution overlaid on top of the stream plot. From the first picture you can see that this periodic solutions are attractive. You might be interested in studying this phenomenon, related to the stability of the system
– Federico
Nov 20 at 18:05
1
1
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
A bit harder to see in that plot, but there's also another (unstable) periodic solution with $y(0) approx 2.604965594$.
– Robert Israel
Nov 20 at 19:48
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
@RobertIsrael Updated to the plot, thanks
– Federico
Nov 21 at 15:24
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3006597%2funsolvable-first-order-nonlinear-differential-equation-fracdxtdt-sin%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
2
No CAS being able to solve is a very good sign that you can't neither. Due to the non-linearity (sine), I doubt that there is a closed-form expression. For small angles, use $sin xapprox x$.
– Yves Daoust
Nov 20 at 17:21
What about assumptions? Could that help you to find a solution? Actually I can not believe that such an equation is not solvable.
– R. Caloudis
Nov 20 at 17:24
3
It's very likely that there are no closed-form solutions to this differential equation. Most differential equations are like that. Of course there are solutions. You can solve the differential equation numerically, or find arbitrarily many terms of a series.
– Robert Israel
Nov 20 at 17:24
Thanks. However, I cannot use the small-angle approximation, because I do need that term in order to describe a certain physical phenomena. Any more guesses?
– R. Caloudis
Nov 20 at 17:26
The problem is the following: You do not have saturation nor harmonics without the sinus term. If you excite with a sinusoid at high amplitude the response would contain only the fundamental frequency without the term.
– R. Caloudis
Nov 20 at 17:30