function pendulum_double_small ( )

%*****************************************************************************80
%
%% pendulum_double_small runs the small perturbation test.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    16 November 2020
%
%  Author:  
%
%    John Burkardt.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'pendulum_double_small\n' );
  fprintf ( 1, '  MATLAB/Octave version %s\n', version ( ) );
  fprintf ( 1, '  Animate the double pendulum small perturbation problem.\n' );
%
%  Set initial conditions.
%
  theta1 = 0.25;
  dtheta1dt = 0.0;
  theta2 = 0.0;
  dtheta2dt = 0.0;
%
%  Set problem parameters.
%
  [ g, m1, m2, l1, l2 ] = pendulum_double_parameters ( );
%
%  Set movie parameters.
%
  t0 = 0.0;
  tstop = 100;
  tspan = [ t0, tstop ];

  fps = 10;
  movie = true;
  label = 'pendulum_double_small';
%
%  Pack initial conditions and parameters into one array.
%
  y0 = [ theta1; dtheta1dt; theta2; dtheta2dt ];
%
%  Call for solution and movie.
%
  pendulum_double_movie ( tspan, y0, fps, movie, label );

  return
end

