Problems solving this differential equation for a Pendulum with Mathematica

In summary, The conversation discusses difficulties with solving a specific differential equation in Mathematica. The user shares their code and asks for help in identifying the issue. Another user suggests adjusting the initial conditions and provides a solution for differentiating between the two functions in a plot.
  • #1
Lambda96
161
60
TL;DR Summary
Mathematica does not solve the differential equations
Hi,

unfortunately, I have problems that Mathematica does not solve the differential equation. The task is as follows and it is about the task c

Bildschirmfoto 2023-06-28 um 18.54.02.png


In the Mathematica Notebook, the following was written for task c

"You can use the following two lines of code to produce the solutions of the exact and approximated differential equations with boundary conditions. After that, write a plot[ ] function to sketch the two together."

exact[t_] := Evaluate[phi[t] /. NDSolve[..., phi[t], {t, 0, 20}]];
approx[t_] := Evaluate[phi[t] /. DSolve[..., phi[t], t]];I then tried to calculate the differential equation for equation (4) using the two codes and got the following

Bildschirmfoto 2023-06-28 um 18.58.49.png

I then tried solving the differential equation using only NDSolve and get the following.

Bildschirmfoto 2023-06-28 um 19.04.34.png

I am a beginner in Mathematica, so I do not know what I have done wrong in the two formulas?
 
Physics news on Phys.org
  • #2
As the error messages say, your initial conditions are wrong. You have
Code:
phi[0]==0
and
Code:
phi[0]==Pi/2
which is not meaningful. You should have the first derivative of phi for one of those.
Code:
phi'[0]==...
 
  • #3
Thanks Dale for your help and for looking over my code 👍

I have now tried plotting the differential equations again, but this time with phi'[0]=0 and with ##\gamma=0.1## and got the following.

Bildschirmfoto 2023-06-28 um 21.44.11.png

Is there any way that I can show in the plot which function is which?
 
  • #4
Lambda96 said:
Is there any way that I can show in the plot which function is which?
There are different options. Try for instance
Code:
Plot[{exact[t],approx[t]},{t,0,20},PlotLegends->"Expressions"]
 
  • #5
Thanks DrClaude for your help, this is exactly what I was looking for 👍
 
  • Like
Likes berkeman

Similar threads

  • MATLAB, Maple, Mathematica, LaTeX
Replies
1
Views
173
  • MATLAB, Maple, Mathematica, LaTeX
Replies
1
Views
1K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
2
Views
297
Replies
3
Views
1K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
2
Views
1K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
4
Views
2K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
5
Views
1K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
1
Views
1K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
1
Views
250
  • MATLAB, Maple, Mathematica, LaTeX
Replies
4
Views
2K
Back
Top