//
// pendulum.sci
//
// Scilab vector field functions for: pendulum
//
// This file was generated by the program VFGEN (Version:2.4.0)
// Generated on 10-Jul-2008 at 17:12
//

//
// pendulum_vf
//
// The vector field 
//
function vf_ = pendulum_vf(t,x_,p_)
    Pi = %pi;
    theta      = x_(1);
    v          = x_(2);
    g          = p_(1);
    b          = p_(2);
    L          = p_(3);
    m          = p_(4);
    vf_ = zeros(2,1);
    vf_(1) = v;
    vf_(2) = -m^(-1)*v*L^(-2)*b-sin(theta)*g*L^(-1);
endfunction

//
// pendulum_jac
//
// The Jacobian of the vector field
//
function jac_ = pendulum_jac(t,x_,p_)
    Pi = %pi;
    theta      = x_(1);
    v          = x_(2);
    g          = p_(1);
    b          = p_(2);
    L          = p_(3);
    m          = p_(4);
    jac_ = zeros(2,2);
    jac_(1,2) = 1;
    jac_(2,1) = -g*L^(-1)*cos(theta);
    jac_(2,2) = -m^(-1)*L^(-2)*b;
endfunction

//
// pendulum_energy
//
// This function implements the user-defined function "energy"
//
function r_ = pendulum_energy(t,x_,p_)
    Pi = %pi;
    theta      = x_(1);
    v          = x_(2);
    g          = p_(1);
    b          = p_(2);
    L          = p_(3);
    m          = p_(4);
    r_ = -m*g*L*cos(theta)+1/2*m*v^2*L^2;
endfunction