Graphics Programs Reference
In-Depth Information
10
8
6
4
2
0
−
2
−
4
−
6
−
2
0
2
4
6
8
10
12
14
16
Initial velocities 5, 6, 7 are not large enoughfor the angle to increase past
π
,
but initial velocities 8, 9, 10 are enoughto make the pendulum swing over-
head. Let's see what happens between 7 and 8.
hold on
for a = 7.0:0.2:8.0
[t, xa] = ode45(g, [0 20], [0 a]);
plot(xa(:, 1), xa(:, 2))
end
hold off
10
8
6
4
2
0
−
2
−
4
−
6
−
2
0
2
4
6
8
10
12
14
16
We see that the cutoff is somewhere between 7.2 and 7.4. Let's make one
more refinement.
hold on
for a = 7.2:0.05:7.4
[t, xa] = ode45(g, [0 20], [0 a]);
Search WWH ::
Custom Search