global c
hold off
options = odeset('OutputFcn', 'odephas2', 'Vectorized', 'on');

c = 0.25


ode45('pendu', [-20 50], [-9 4.115], options)
axis([-10 15 -3 3])
hold on

ode45('pendu', [-20 50], [12 -3.741], options)
ode45('pendu', [-20 50], [pi pi], options)
ode45('pendu', [-20 50], [pi -pi], options)

pause, hold off

c = 0.5

ode45('pendu', [-20 50], [0 3.085], options)
axis([-2 8 -3 3])
hold
ode45('pendu', [-20 50], [6 -2.93], options)
ode45('pendu', [-20 50], [pi pi], options)
ode45('pendu', [-20 50], [pi -pi], options)

