% matlab script to solve pendulum equation
% Ascher and Petzold p.5

clear all, close all

global m g

% problem parameters
m = 1;
g = 9.8;
theta0 = 1;
omega0 = 1;
y0 = [theta0,omega0]';

t0 = 0;
tf = 3;

[t,y] = ode45('pendulum_rhs',[t0,tf],y0);

plot(t,y(:,1))