/*
 *  MorrisLecar_taylor15.c
 *
 *  C file with functions for computing the Taylor series approximate solution
 *  for the vector field named: MorrisLecar
 *
 *  This file was generated by the program VFGEN (Version:2.4.0)
 *  Generated on 10-Jul-2008 at 18:00
 */

#include <math.h>
#include "MorrisLecar_taylor15.h"

/*
 *  The vector field.
 */

void MorrisLecar_vf(double t, const double y_[], double f_[], double params[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;
    double minf, winf, tauw;
    double *p_;

    p_ = (double *) params;

    v          = y_[0];
    w          = y_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    minf =  tanh(-( v1-v)/v2)/2.0+(1.0/2.0);
    winf =  tanh(-( v3-v)/v4)/2.0+(1.0/2.0);
    tauw = 1.0/(cosh(-(1.0/2.0)*( v3-v)/v4));

    f_[0] = -1.0/c*( w*( v-vk)*gk+gl*( v-vl)-gca*minf*( vca-v)-ic);
    f_[1] = -1.0/tauw*( w-winf)*phi;

    }


/*
 *     deriv = Df(x)[v1]
 */
void MorrisLecar_diff1(double deriv[],double x_[], double p_[],double v1_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v = -(1.0/2.0)*1.0/c*( 2.0*gl+2.0*w*gk+gca*( tanh(-( v1-v)/v2)+1.0)+gca*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/v2);
    double vf_w = -1.0/c*( v-vk)*gk;
    deriv[0] = 0.0;
    deriv[0] += vf_v * v1_[0];
    deriv[0] += vf_w * v1_[1];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v = -( pow(tanh(-( v3-v)/v4),2.0)-1.0)/v4*cosh(-(1.0/2.0)*( v3-v)/v4)*phi/2.0-sinh(-(1.0/2.0)*( v3-v)/v4)/v4*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/4.0;
    double vf_w = -cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v * v1_[0];
    deriv[1] += vf_w * v1_[1];
    }

    return;
    }

/*
 *     deriv = D^2f(x)[v1,v2]
 */
void MorrisLecar_diff2(double deriv[],double x_[], double p_[],double v1_[],double v2_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v = 1.0/c*( gca*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/v2+gca*tanh(-( v1-v)/v2)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*v2));
    double vf_v_w = -1.0/c*gk;
    deriv[0] = 0.0;
    deriv[0] += vf_v_v * v1_[0]*v2_[0];
    deriv[0] += vf_v_w * v1_[0]*v2_[1];
    deriv[0] += vf_v_w * v1_[1]*v2_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v = -sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*v4)*phi/2.0+( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*v4)*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)-1.0/(v4*v4)*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/8.0;
    double vf_v_w = -(1.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/v4*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v * v1_[0]*v2_[0];
    deriv[1] += vf_v_w * v1_[0]*v2_[1];
    deriv[1] += vf_v_w * v1_[1]*v2_[0];
    }

    return;
    }

/*
 *     deriv = D^3f(x)[v1,v2,v3]
 */
void MorrisLecar_diff3(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v = -1.0/c*( gca*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*v2*v2)+3.0*gca*tanh(-( v1-v)/v2)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*v2)+2.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*v2*v2));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v * v1_[0]*v2_[0]*v3_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v =  -2.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*v4*v4)*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-(3.0/8.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*v4*v4)*cosh(-(1.0/2.0)*( v3-v)/v4)*phi-sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*v4*v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/16.0+(3.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*v4*v4)*phi*tanh(-( v3-v)/v4)-pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*v4*v4)*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    double vf_v_v_w = -(1.0/4.0)*1.0/(v4*v4)*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v * v1_[0]*v2_[0]*v3_[0];
    deriv[1] += vf_v_v_w * v1_[0]*v2_[0]*v3_[1];
    deriv[1] += vf_v_v_w * v1_[0]*v2_[1]*v3_[0];
    deriv[1] += vf_v_v_w * v1_[1]*v2_[0]*v3_[0];
    }

    return;
    }

/*
 *     deriv = D^4f(x)[v1,...,v4]
 */
void MorrisLecar_diff4(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v = 4.0*1.0/c*( 2.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*v2*v2)+gca*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*v2*v2)+gca*pow(tanh(-( v1-v)/v2),3.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*v2)*(v2*v2))+2.0*gca*tanh(-( v1-v)/v2)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*v2)*(v2*v2)));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v =  -2.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*v4)*(v4*v4))*phi-1.0/((v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/32.0+-4.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4)*(v4*v4))*phi*pow(tanh(-( v3-v)/v4),2.0)+4.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+8.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)-sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4)*(v4*v4))*phi/4.0+(3.0/2.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4);
    double vf_v_v_v_w = -(1.0/8.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*v4*v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0];
    deriv[1] += vf_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1];
    deriv[1] += vf_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0];
    deriv[1] += vf_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0];
    deriv[1] += vf_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0];
    }

    return;
    }

/*
 *     deriv = D^5f(x)[v1,...,v5]
 */
void MorrisLecar_diff5(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v = -4.0*1.0/c*( 2.0*gca*pow(tanh(-( v1-v)/v2),4.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*(v2*v2)*(v2*v2))+10.0*gca*tanh(-( v1-v)/v2)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*v2)*(v2*v2))+11.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*(v2*v2)*(v2*v2))+2.0*gca*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*(v2*v2)*(v2*v2))+5.0*gca*pow(tanh(-( v1-v)/v2),3.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*v2)*(v2*v2)));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v =  -5.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+10.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4)*(v4*v4))*phi*pow(tanh(-( v3-v)/v4),3.0)+-44.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-8.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-(5.0/2.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+(5.0/4.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4)*(v4*v4))*phi*tanh(-( v3-v)/v4)-sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*(v4*v4)*(v4*v4))*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/64.0+-(5.0/32.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-8.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+20.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4)*(v4*v4))*phi*tanh(-( v3-v)/v4);
    double vf_v_v_v_v_w = -(1.0/16.0)*1.0/((v4*v4)*(v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0];
    deriv[1] += vf_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1];
    deriv[1] += vf_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0];
    deriv[1] += vf_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0];
    deriv[1] += vf_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0];
    deriv[1] += vf_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0];
    }

    return;
    }

/*
 *     deriv = D^6f(x)[v1,...,v6]
 */
void MorrisLecar_diff6(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v = 8.0*1.0/c*( 6.0*gca*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*(v2*v2)*(v2*v2))+26.0*gca*pow(tanh(-( v1-v)/v2),3.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*v2*v2)*(v2*v2*v2))+6.0*gca*pow(tanh(-( v1-v)/v2),4.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*(v2*v2)*(v2*v2))+17.0*gca*tanh(-( v1-v)/v2)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/((v2*v2*v2)*(v2*v2*v2))+33.0*gca*pow(tanh(-( v1-v)/v2),2.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*(v2*v2)*(v2*v2))+2.0*gca*pow(tanh(-( v1-v)/v2),5.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*v2*v2)*(v2*v2*v2)));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v =  136.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-24.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4*v4)*(v4*v4*v4))*phi*pow(tanh(-( v3-v)/v4),4.0)+-132.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*v4*v4)*(v4*v4*v4))*phi*pow(tanh(-( v3-v)/v4),2.0)+16.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)-1.0/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/128.0+-(5.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*v4*v4)*(v4*v4*v4))*phi+208.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-(3.0/32.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4*v4)*(v4*v4*v4))*phi+15.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-5.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4*v4)*(v4*v4*v4))*phi*pow(tanh(-( v3-v)/v4),2.0)+-24.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*v4*v4)*(v4*v4*v4))*phi+30.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+(15.0/16.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4);
    double vf_v_v_v_v_v_w = -(1.0/32.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*(v4*v4)*(v4*v4))*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0];
    deriv[1] += vf_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1];
    deriv[1] += vf_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0];
    deriv[1] += vf_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0];
    deriv[1] += vf_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0];
    deriv[1] += vf_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0];
    deriv[1] += vf_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0];
    }

    return;
    }

/*
 *     deriv = D^7f(x)[v1,...,v7]
 */
void MorrisLecar_diff7(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v = -8.0*1.0/c*( 180.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*(v2*v2*v2)*(v2*v2*v2))+182.0*gca*pow(tanh(-( v1-v)/v2),3.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*v2*v2)*(v2*v2*v2))+4.0*gca*pow(tanh(-( v1-v)/v2),6.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*(v2*v2*v2)*(v2*v2*v2))+119.0*gca*tanh(-( v1-v)/v2)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/((v2*v2*v2)*(v2*v2*v2))+14.0*gca*pow(tanh(-( v1-v)/v2),5.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*v2*v2)*(v2*v2*v2))+17.0*gca*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*(v2*v2*v2)*(v2*v2*v2))+114.0*gca*pow(tanh(-( v1-v)/v2),4.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*(v2*v2*v2)*(v2*v2*v2)));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v =  728.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi*pow(tanh(-( v3-v)/v4),3.0)+476.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi*tanh(-( v3-v)/v4)+-231.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+35.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi*tanh(-( v3-v)/v4)+-(35.0/8.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-136.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+56.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi*pow(tanh(-( v3-v)/v4),5.0)+-(7.0/128.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-912.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-1440.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-42.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-32.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+(35.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi*pow(tanh(-( v3-v)/v4),3.0)+-42.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+(21.0/32.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi*tanh(-( v3-v)/v4)-sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*(v4*v4*v4)*(v4*v4*v4))*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/256.0+-(35.0/16.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    double vf_v_v_v_v_v_v_w = -(1.0/64.0)*1.0/((v4*v4*v4)*(v4*v4*v4))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0];
    deriv[1] += vf_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0];
    }

    return;
    }

/*
 *     deriv = D^8f(x)[v1,...,v8]
 */
void MorrisLecar_diff8(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v = 64.0*1.0/c*( gca*pow(tanh(-( v1-v)/v2),7.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+4.0*gca*pow(tanh(-( v1-v)/v2),6.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*(v2*v2*v2)*(v2*v2*v2))+17.0*gca*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*(v2*v2*v2)*(v2*v2*v2))+60.0*gca*pow(tanh(-( v1-v)/v2),5.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+114.0*gca*pow(tanh(-( v1-v)/v2),4.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*(v2*v2*v2)*(v2*v2*v2))+192.0*gca*pow(tanh(-( v1-v)/v2),3.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+62.0*gca*tanh(-( v1-v)/v2)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+180.0*gca*pow(tanh(-( v1-v)/v2),2.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*(v2*v2*v2)*(v2*v2*v2)));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v =  -308.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-544.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi+(35.0/2.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+1456.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+(7.0/16.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-56.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+-(7.0/4.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi+3968.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-128.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+-5760.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+35.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)-sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi/32.0+-(7.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+952.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+3840.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+-3648.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+64.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+112.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+12288.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-56.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi-1.0/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/512.0;
    double vf_v_v_v_v_v_v_v_w = -(1.0/128.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*(v4*v4*v4)*(v4*v4*v4))*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0];
    }

    return;
    }

/*
 *     deriv = D^9f(x)[v1,...,v9]
 */
void MorrisLecar_diff9(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[],double v9_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v = -64.0*1.0/c*( 62.0*gca*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+540.0*gca*pow(tanh(-( v1-v)/v2),5.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+558.0*gca*tanh(-( v1-v)/v2)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+9.0*gca*pow(tanh(-( v1-v)/v2),7.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+1072.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+247.0*gca*pow(tanh(-( v1-v)/v2),6.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+1728.0*gca*pow(tanh(-( v1-v)/v2),3.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+2.0*gca*pow(tanh(-( v1-v)/v2),8.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+1452.0*gca*pow(tanh(-( v1-v)/v2),4.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2))));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v =  17856.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+-15808.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)-sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/1024.0+(9.0/32.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+-63.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-12960.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-8208.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-(21.0/16.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-3968.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-288.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+17280.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+(63.0/4.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-63.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+288.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),7.0)+(63.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+168.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+-(9.0/512.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-128.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),8.0)+-1224.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+55296.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-(21.0/8.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+2184.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+1428.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+-68608.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-92928.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-(693.0/2.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0);
    double vf_v_v_v_v_v_v_v_v_w = -(1.0/256.0)*1.0/(((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0];
    }

    return;
    }

/*
 *     deriv = D^10f(x)[v1,...,v10]
 */
void MorrisLecar_diff10(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[],double v9_[],double v10_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v = 128.0*1.0/c*( 1235.0*gca*pow(tanh(-( v1-v)/v2),6.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+7192.0*gca*pow(tanh(-( v1-v)/v2),3.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+5360.0*gca*pow(tanh(-( v1-v)/v2),2.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+2.0*gca*pow(tanh(-( v1-v)/v2),9.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+310.0*gca*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+10.0*gca*pow(tanh(-( v1-v)/v2),8.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+1382.0*gca*tanh(-( v1-v)/v2)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+502.0*gca*pow(tanh(-( v1-v)/v2),7.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+7260.0*gca*pow(tanh(-( v1-v)/v2),4.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*((v2*v2)*(v2*v2))*((v2*v2)*(v2*v2)))+5097.0*gca*pow(tanh(-( v1-v)/v2),5.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2))));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v =  -13680.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+-63.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi+1785.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-640.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),8.0)+-79040.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+-(15.0/8.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+210.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+138240.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)-1.0/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/2048.0+-343040.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+64256.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-464640.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+720.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-(5.0/512.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi+43200.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+920576.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-63.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+256.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),9.0)+2730.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-19840.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi+-(15.0/16.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi+652416.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+(105.0/8.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-480.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+-21600.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+(105.0/4.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+176896.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+44640.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-(693.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-2040.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi+(45.0/256.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4);
    double vf_v_v_v_v_v_v_v_v_v_w = -(1.0/512.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*((v4*v4)*(v4*v4))*((v4*v4)*(v4*v4)))*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[1]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0];
    }

    return;
    }

/*
 *     deriv = D^11f(x)[v1,...,v11]
 */
void MorrisLecar_diff11(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[],double v9_[],double v10_[],double v11_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v = -128.0*1.0/c*( 56067.0*gca*pow(tanh(-( v1-v)/v2),5.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+2026.0*gca*pow(tanh(-( v1-v)/v2),8.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+15202.0*gca*tanh(-( v1-v)/v2)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+83021.0*gca*pow(tanh(-( v1-v)/v2),4.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+79112.0*gca*pow(tanh(-( v1-v)/v2),3.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+1382.0*gca*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+4.0*gca*pow(tanh(-( v1-v)/v2),10.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+35396.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+22.0*gca*pow(tanh(-( v1-v)/v2),9.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+5522.0*gca*pow(tanh(-( v1-v)/v2),7.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+34096.0*gca*pow(tanh(-( v1-v)/v2),6.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2))));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v =  (165.0/16.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-259328.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),8.0)+3588288.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+-(231.0/4.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+5063168.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-1760.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),8.0)+1408.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),9.0)+81840.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+-(165.0/128.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+(55.0/512.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)-sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/4096.0+353408.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),7.0)+-4364288.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+-943360.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-54560.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-(165.0/256.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-(2541.0/8.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-4530688.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+79200.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+-1277760.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-217360.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+-(231.0/4.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+3003.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+(3927.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+-18810.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-29700.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-10626688.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-660.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+(165.0/8.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+1320.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),7.0)+972928.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*tanh(-( v3-v)/v4)+-512.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),10.0)+-176896.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+231.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+-(11.0/2048.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+253440.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-2805.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    double vf_v_v_v_v_v_v_v_v_v_v_w = -(1.0/1024.0)*1.0/((v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[1]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[1]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0];
    }

    return;
    }

/*
 *     deriv = D^12f(x)[v1,...,v12]
 */
void MorrisLecar_diff12(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[],double v9_[],double v10_[],double v11_[],double v12_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v_v = 512.0*1.0/c*( 217186.0*gca*pow(tanh(-( v1-v)/v2),5.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+2.0*gca*pow(tanh(-( v1-v)/v2),11.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+249063.0*gca*pow(tanh(-( v1-v)/v2),4.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+106188.0*gca*pow(tanh(-( v1-v)/v2),2.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+102288.0*gca*pow(tanh(-( v1-v)/v2),6.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+6078.0*gca*pow(tanh(-( v1-v)/v2),8.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+171511.0*gca*pow(tanh(-( v1-v)/v2),3.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+12.0*gca*pow(tanh(-( v1-v)/v2),10.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2)))+2036.0*gca*pow(tanh(-( v1-v)/v2),9.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+55196.0*gca*pow(tanh(-( v1-v)/v2),7.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+21844.0*gca*tanh(-( v1-v)/v2)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+4146.0*gca*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/(v2*(v2*(v2*v2)*(v2*v2))*(v2*(v2*v2)*(v2*v2))));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v_v =  3003.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+4224.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),9.0)+-26185728.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+-(55.0/64.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+87813632.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+15189504.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+1024.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),11.0)+28260352.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-1886720.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-2555520.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+(495.0/64.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-(99.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi+10764864.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+1042432.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),9.0)+-22572.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+-(3.0/1024.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi+122760.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-1061376.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi+2918784.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+(495.0/32.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-434720.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+-3520.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),8.0)+(33.0/512.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+11184128.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-(55.0/128.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi+-109120.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi+(3927.0/2.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-27184128.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-(1089.0/4.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-3072.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),10.0)+-3366.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi+1060224.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-63760128.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)-1.0/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/8192.0+380160.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+231.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+-(99.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+118800.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+1980.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-1555968.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),8.0)+111199232.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+-35640.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-792.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0);
    double vf_v_v_v_v_v_v_v_v_v_v_v_w = -(1.0/2048.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*(v4*(v4*v4)*(v4*v4))*(v4*(v4*v4)*(v4*v4)))*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[1]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[1]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[1]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0];
    }

    return;
    }

/*
 *     deriv = D^13f(x)[v1,...,v13]
 */
void MorrisLecar_diff13(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[],double v9_[],double v10_[],double v11_[],double v12_[],double v13_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v_v_v = -512.0*1.0/c*( 2801040.0*gca*pow(tanh(-( v1-v)/v2),4.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+26468.0*gca*pow(tanh(-( v1-v)/v2),9.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+2123860.0*gca*pow(tanh(-( v1-v)/v2),6.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+21844.0*gca*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,7.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+283972.0*gca*tanh(-( v1-v)/v2)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+349500.0*gca*pow(tanh(-( v1-v)/v2),8.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+2229643.0*gca*pow(tanh(-( v1-v)/v2),3.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+8166.0*gca*pow(tanh(-( v1-v)/v2),10.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+26.0*gca*pow(tanh(-( v1-v)/v2),11.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+717548.0*gca*pow(tanh(-( v1-v)/v2),7.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+776661.0*gca*pow(tanh(-( v1-v)/v2),2.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+2823418.0*gca*pow(tanh(-( v1-v)/v2),5.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+4.0*gca*pow(tanh(-( v1-v)/v2),12.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2))));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v_v_v =  -1.0874163200000000e+09*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+(715.0/64.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*tanh(-( v3-v)/v4)+2574.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),7.0)+-(1287.0/32.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+6324032.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*tanh(-( v3-v)/v4)+-706420.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+-207220416.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-858.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+-(13.0/8192.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-3449472.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+(429.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+494208.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-24453.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+-(7293.0/2.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-1.4341324800000000e+09*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+7.2279500800000000e+08*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+-38610.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-9984.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),10.0)+-2048.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),12.0)+-178944000.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),8.0)+5.7078860800000000e+08*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+(5577.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+6775808.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),9.0)+(7293.0/4.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*tanh(-( v3-v)/v4)+-4180992.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),10.0)+-5056896.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),8.0)+159588.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*tanh(-( v3-v)/v4)+72696832.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*tanh(-( v3-v)/v4)+(39.0/1024.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*tanh(-( v3-v)/v4)+2297152.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),7.0)-sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/16384.0+-5720.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),8.0)+154440.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+-177320.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+-(143.0/512.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+(715.0/128.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-397650432.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-85103616.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),6.0)+-(1287.0/32.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+32910592.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),3.0)+-(143.0/256.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-11184128.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,7.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi+23323872.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),5.0)+6656.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),11.0)+-(14157.0/64.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-88348416.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+-4152720.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),4.0)+9152.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),9.0)+-3065920.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),2.0)+183692288.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),7.0);
    double vf_v_v_v_v_v_v_v_v_v_v_v_v_w = -(1.0/4096.0)*1.0/(((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[1]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[1]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[1]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[1]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0];
    }

    return;
    }

/*
 *     deriv = D^14f(x)[v1,...,v14]
 */
void MorrisLecar_diff14(double deriv[],double x_[], double p_[],double v1_[],double v2_[],double v3_[],double v4_[],double v5_[],double v6_[],double v7_[],double v8_[],double v9_[],double v10_[],double v11_[],double v12_[],double v13_[],double v14_[])
    {
    double v, w;
    double gca, gk, gl, vca, vk, vl, c, phi, ic, v1, v2, v3, v4;

    v          = x_[0];
    w          = x_[1];

    gca        = p_[0];
    gk         = p_[1];
    gl         = p_[2];
    vca        = p_[3];
    vk         = p_[4];
    vl         = p_[5];
    c          = p_[6];
    phi        = p_[7];
    ic         = p_[8];
    v1         = p_[9];
    v2         = p_[10];
    v3         = p_[11];
    v4         = p_[12];

    {
    /*
     *  Partial derivatives of vf[0]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v_v_v_v = 1024.0*1.0/c*( 10262046.0*gca*pow(tanh(-( v1-v)/v2),3.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+20376780.0*gca*pow(tanh(-( v1-v)/v2),5.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+5436627.0*gca*pow(tanh(-( v1-v)/v2),2.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,6.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+28.0*gca*pow(tanh(-( v1-v)/v2),12.0)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+1089330.0*gca*pow(tanh(-( v1-v)/v2),9.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+2446500.0*gca*pow(tanh(-( v1-v)/v2),8.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,3.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+19607280.0*gca*pow(tanh(-( v1-v)/v2),4.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,5.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+9893440.0*gca*pow(tanh(-( v1-v)/v2),7.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+14867020.0*gca*pow(tanh(-( v1-v)/v2),6.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,4.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+57162.0*gca*pow(tanh(-( v1-v)/v2),10.0)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2)))+4.0*gca*pow(tanh(-( v1-v)/v2),13.0)*( vca-v)*( pow(tanh(-( v1-v)/v2),2.0)-1.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+929569.0*gca*tanh(-( v1-v)/v2)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,7.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+16356.0*gca*pow(tanh(-( v1-v)/v2),11.0)*( vca-v)*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,2.0)/((v2*(v2*v2*v2)*(v2*v2*v2))*(v2*(v2*v2*v2)*(v2*v2*v2)))+152908.0*gca*pow( pow(tanh(-( v1-v)/v2),2.0)-1.0,7.0)/(v2*((v2*v2*v2)*(v2*v2*v2))*((v2*v2*v2)*(v2*v2*v2))));
    deriv[0] = 0.0;
    deriv[0] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    }
    {
    /*
     *  Partial derivatives of vf[1]. 
     *  Any derivative not listed here is zero.
     */
    double vf_v_v_v_v_v_v_v_v_v_v_v_v_v_v =  3003.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-38610.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-858.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+-11799424.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),8.0)+2.5297825280000000e+09*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+2.0865822720000000e+10*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+-7.6119142400000000e+09*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+180180.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+1.1154739200000000e+09*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),9.0)-1.0/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*( 2.0*w-tanh(-( v3-v)/v4)-1.0)*phi/32768.0+-(1001.0/32.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+1.0130882560000000e+10*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+576576.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+(3003.0/16.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+4020016.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+-8048768.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+-(7.0/8192.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+-483514304.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+-78288896.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,7.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+-(1001.0/32.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+-1.2526080000000000e+09*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),8.0)+-23296.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),10.0)+-206146304.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+9.5187865600000000e+08*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,7.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-(11011.0/64.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-14336.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),12.0)+(51051.0/32.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+254438912.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-8008.0*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),8.0)+-988988.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+(91.0/4096.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+11067056.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+(1001.0/128.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+186186.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*tanh(-( v3-v)/v4)+-1.0038927360000000e+10*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+-(7293.0/2.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+-24453.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+(1001.0/256.0)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+23715328.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),9.0)+40816776.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),5.0)+-5813808.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),4.0)+16748544.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),11.0)+4096.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),13.0)+1.0508335104000000e+10*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-4292288.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+57593536.0*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,4.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+23296.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),11.0)+6.4292300800000000e+08*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),7.0)+(39039.0/16.0)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-198575104.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,3.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),6.0)+16016.0*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),9.0)+-(91.0/256.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*( pow(tanh(-( v3-v)/v4),2.0)-1.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-(91.0/512.0)*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+-2.7835530240000000e+09*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,6.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),2.0)+-248248.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi+1.9977601280000000e+09*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,5.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*cosh(-(1.0/2.0)*( v3-v)/v4)*phi*pow(tanh(-( v3-v)/v4),3.0)+-29266944.0*sinh(-(1.0/2.0)*( v3-v)/v4)*pow( pow(tanh(-( v3-v)/v4),2.0)-1.0,2.0)/((v4*(v4*v4*v4)*(v4*v4*v4))*(v4*(v4*v4*v4)*(v4*v4*v4)))*phi*pow(tanh(-( v3-v)/v4),10.0);
    double vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w = -(1.0/8192.0)*sinh(-(1.0/2.0)*( v3-v)/v4)/(v4*((v4*v4*v4)*(v4*v4*v4))*((v4*v4*v4)*(v4*v4*v4)))*phi;
    deriv[1] = 0.0;
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_v * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[1];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[1]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[1]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[1]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[1]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[1]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[1]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[1]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[1]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[0]*v5_[1]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[0]*v4_[1]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[0]*v3_[1]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[0]*v2_[1]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    deriv[1] += vf_v_v_v_v_v_v_v_v_v_v_v_v_v_w * v1_[1]*v2_[0]*v3_[0]*v4_[0]*v5_[0]*v6_[0]*v7_[0]*v8_[0]*v9_[0]*v10_[0]*v11_[0]*v12_[0]*v13_[0]*v14_[0];
    }

    return;
    }

/*
 *  MorrisLecar_derivs15
 *
 *  Compute the coefficients in the Taylor polynomial at X.
 *  These are just the derivatives; they have not been scaled
 *  by the appropriate factorial.
 *
 */

void MorrisLecar_derivs15(double Xderiv[15][2], double X[], double params[])
    {
    int i;
    double s;
    double Q[2];

    MorrisLecar_vf(0.0,X,Xderiv[0],params);

    for (i = 0; i < 2; ++i)
        Xderiv[1][i] = 0.0;
    /*    [ 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[1][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[2][i] = 0.0;
    /*    [ 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[2][i] += Q[i];
    /*    [ 2]  coeff = 1.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[2][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[3][i] = 0.0;
    /*    [ 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[3][i] += Q[i];
    /*    [ 1 1]  coeff = 3.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[3][i] += 3.*Q[i];
    /*    [ 3]  coeff = 1.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[3][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[4][i] = 0.0;
    /*    [ 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[4][i] += Q[i];
    /*    [ 1 0 1]  coeff = 4.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[4][i] += 4.*Q[i];
    /*    [ 0 2]  coeff = 3.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[4][i] += 3.*Q[i];
    /*    [ 2 1]  coeff = 6.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[4][i] += 6.*Q[i];
    /*    [ 4]  coeff = 1.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[4][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[5][i] = 0.0;
    /*    [ 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += Q[i];
    /*    [ 1 0 0 1]  coeff = 5.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += 5.*Q[i];
    /*    [ 0 1 1]  coeff = 10.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += 10.*Q[i];
    /*    [ 2 0 1]  coeff = 10.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += 10.*Q[i];
    /*    [ 1 2]  coeff = 15.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += 15.*Q[i];
    /*    [ 3 1]  coeff = 10.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += 10.*Q[i];
    /*    [ 5]  coeff = 1.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[5][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[6][i] = 0.0;
    /*    [ 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += Q[i];
    /*    [ 1 0 0 0 1]  coeff = 6.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 6.*Q[i];
    /*    [ 0 1 0 1]  coeff = 15.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 15.*Q[i];
    /*    [ 2 0 0 1]  coeff = 15.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 15.*Q[i];
    /*    [ 0 0 2]  coeff = 10.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 10.*Q[i];
    /*    [ 1 1 1]  coeff = 60.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 60.*Q[i];
    /*    [ 3 0 1]  coeff = 20.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 20.*Q[i];
    /*    [ 0 3]  coeff = 15.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 15.*Q[i];
    /*    [ 2 2]  coeff = 45.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 45.*Q[i];
    /*    [ 4 1]  coeff = 15.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += 15.*Q[i];
    /*    [ 6]  coeff = 1.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[6][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[7][i] = 0.0;
    /*    [ 0 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += Q[i];
    /*    [ 1 0 0 0 0 1]  coeff = 7.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 7.*Q[i];
    /*    [ 0 1 0 0 1]  coeff = 21.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 21.*Q[i];
    /*    [ 2 0 0 0 1]  coeff = 21.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 21.*Q[i];
    /*    [ 0 0 1 1]  coeff = 35.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 35.*Q[i];
    /*    [ 1 1 0 1]  coeff = 105.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 105.*Q[i];
    /*    [ 3 0 0 1]  coeff = 35.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 35.*Q[i];
    /*    [ 0 2 1]  coeff = 105.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 105.*Q[i];
    /*    [ 1 0 2]  coeff = 70.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 70.*Q[i];
    /*    [ 2 1 1]  coeff = 210.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 210.*Q[i];
    /*    [ 4 0 1]  coeff = 35.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 35.*Q[i];
    /*    [ 1 3]  coeff = 105.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 105.*Q[i];
    /*    [ 3 2]  coeff = 105.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 105.*Q[i];
    /*    [ 5 1]  coeff = 21.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += 21.*Q[i];
    /*    [ 7]  coeff = 1.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[7][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[8][i] = 0.0;
    /*    [ 0 0 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += Q[i];
    /*    [ 1 0 0 0 0 0 1]  coeff = 8.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 8.*Q[i];
    /*    [ 0 1 0 0 0 1]  coeff = 28.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 28.*Q[i];
    /*    [ 2 0 0 0 0 1]  coeff = 28.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 28.*Q[i];
    /*    [ 0 0 1 0 1]  coeff = 56.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 56.*Q[i];
    /*    [ 1 1 0 0 1]  coeff = 168.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 168.*Q[i];
    /*    [ 3 0 0 0 1]  coeff = 56.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 56.*Q[i];
    /*    [ 0 0 0 2]  coeff = 35.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 35.*Q[i];
    /*    [ 0 2 0 1]  coeff = 210.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 210.*Q[i];
    /*    [ 1 0 1 1]  coeff = 280.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 280.*Q[i];
    /*    [ 2 1 0 1]  coeff = 420.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 420.*Q[i];
    /*    [ 4 0 0 1]  coeff = 70.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 70.*Q[i];
    /*    [ 0 1 2]  coeff = 280.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 280.*Q[i];
    /*    [ 1 2 1]  coeff = 840.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 840.*Q[i];
    /*    [ 2 0 2]  coeff = 280.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 280.*Q[i];
    /*    [ 3 1 1]  coeff = 560.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 560.*Q[i];
    /*    [ 5 0 1]  coeff = 56.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 56.*Q[i];
    /*    [ 0 4]  coeff = 105.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 105.*Q[i];
    /*    [ 2 3]  coeff = 420.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 420.*Q[i];
    /*    [ 4 2]  coeff = 210.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 210.*Q[i];
    /*    [ 6 1]  coeff = 28.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += 28.*Q[i];
    /*    [ 8]  coeff = 1.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[8][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[9][i] = 0.0;
    /*    [ 0 0 0 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += Q[i];
    /*    [ 1 0 0 0 0 0 0 1]  coeff = 9.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 9.*Q[i];
    /*    [ 0 1 0 0 0 0 1]  coeff = 36.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 36.*Q[i];
    /*    [ 2 0 0 0 0 0 1]  coeff = 36.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 36.*Q[i];
    /*    [ 0 0 1 0 0 1]  coeff = 84.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 84.*Q[i];
    /*    [ 1 1 0 0 0 1]  coeff = 252.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 252.*Q[i];
    /*    [ 3 0 0 0 0 1]  coeff = 84.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 84.*Q[i];
    /*    [ 0 0 0 1 1]  coeff = 126.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[3],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 126.*Q[i];
    /*    [ 0 2 0 0 1]  coeff = 378.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 378.*Q[i];
    /*    [ 1 0 1 0 1]  coeff = 504.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 504.*Q[i];
    /*    [ 2 1 0 0 1]  coeff = 756.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 756.*Q[i];
    /*    [ 4 0 0 0 1]  coeff = 126.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 126.*Q[i];
    /*    [ 0 1 1 1]  coeff = 1260.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1260.*Q[i];
    /*    [ 1 0 0 2]  coeff = 315.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 315.*Q[i];
    /*    [ 1 2 0 1]  coeff = 1890.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1890.*Q[i];
    /*    [ 2 0 1 1]  coeff = 1260.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1260.*Q[i];
    /*    [ 3 1 0 1]  coeff = 1260.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1260.*Q[i];
    /*    [ 5 0 0 1]  coeff = 126.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 126.*Q[i];
    /*    [ 0 0 3]  coeff = 280.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[2],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 280.*Q[i];
    /*    [ 0 3 1]  coeff = 1260.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1260.*Q[i];
    /*    [ 1 1 2]  coeff = 2520.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 2520.*Q[i];
    /*    [ 2 2 1]  coeff = 3780.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 3780.*Q[i];
    /*    [ 3 0 2]  coeff = 840.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 840.*Q[i];
    /*    [ 4 1 1]  coeff = 1260.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1260.*Q[i];
    /*    [ 6 0 1]  coeff = 84.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 84.*Q[i];
    /*    [ 1 4]  coeff = 945.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 945.*Q[i];
    /*    [ 3 3]  coeff = 1260.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 1260.*Q[i];
    /*    [ 5 2]  coeff = 378.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 378.*Q[i];
    /*    [ 7 1]  coeff = 36.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += 36.*Q[i];
    /*    [ 9]  coeff = 1.  */
    MorrisLecar_diff9(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[9][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[10][i] = 0.0;
    /*    [ 0 0 0 0 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[9]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += Q[i];
    /*    [ 1 0 0 0 0 0 0 0 1]  coeff = 10.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 10.*Q[i];
    /*    [ 0 1 0 0 0 0 0 1]  coeff = 45.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 45.*Q[i];
    /*    [ 2 0 0 0 0 0 0 1]  coeff = 45.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 45.*Q[i];
    /*    [ 0 0 1 0 0 0 1]  coeff = 120.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 120.*Q[i];
    /*    [ 1 1 0 0 0 0 1]  coeff = 360.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 360.*Q[i];
    /*    [ 3 0 0 0 0 0 1]  coeff = 120.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 120.*Q[i];
    /*    [ 0 0 0 1 0 1]  coeff = 210.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[3],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 210.*Q[i];
    /*    [ 0 2 0 0 0 1]  coeff = 630.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 630.*Q[i];
    /*    [ 1 0 1 0 0 1]  coeff = 840.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 840.*Q[i];
    /*    [ 2 1 0 0 0 1]  coeff = 1260.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 1260.*Q[i];
    /*    [ 4 0 0 0 0 1]  coeff = 210.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 210.*Q[i];
    /*    [ 0 0 0 0 2]  coeff = 126.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[4],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 126.*Q[i];
    /*    [ 0 1 1 0 1]  coeff = 2520.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2520.*Q[i];
    /*    [ 1 0 0 1 1]  coeff = 1260.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[3],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 1260.*Q[i];
    /*    [ 1 2 0 0 1]  coeff = 3780.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 3780.*Q[i];
    /*    [ 2 0 1 0 1]  coeff = 2520.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2520.*Q[i];
    /*    [ 3 1 0 0 1]  coeff = 2520.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2520.*Q[i];
    /*    [ 5 0 0 0 1]  coeff = 252.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 252.*Q[i];
    /*    [ 0 0 2 1]  coeff = 2100.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[2],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2100.*Q[i];
    /*    [ 0 1 0 2]  coeff = 1575.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 1575.*Q[i];
    /*    [ 0 3 0 1]  coeff = 3150.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 3150.*Q[i];
    /*    [ 1 1 1 1]  coeff = 12600.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 12600.*Q[i];
    /*    [ 2 0 0 2]  coeff = 1575.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 1575.*Q[i];
    /*    [ 2 2 0 1]  coeff = 9450.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 9450.*Q[i];
    /*    [ 3 0 1 1]  coeff = 4200.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 4200.*Q[i];
    /*    [ 4 1 0 1]  coeff = 3150.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 3150.*Q[i];
    /*    [ 6 0 0 1]  coeff = 210.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 210.*Q[i];
    /*    [ 0 2 2]  coeff = 6300.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 6300.*Q[i];
    /*    [ 1 0 3]  coeff = 2800.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2800.*Q[i];
    /*    [ 1 3 1]  coeff = 12600.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 12600.*Q[i];
    /*    [ 2 1 2]  coeff = 12600.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 12600.*Q[i];
    /*    [ 3 2 1]  coeff = 12600.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 12600.*Q[i];
    /*    [ 4 0 2]  coeff = 2100.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2100.*Q[i];
    /*    [ 5 1 1]  coeff = 2520.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 2520.*Q[i];
    /*    [ 7 0 1]  coeff = 120.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 120.*Q[i];
    /*    [ 0 5]  coeff = 945.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 945.*Q[i];
    /*    [ 2 4]  coeff = 4725.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 4725.*Q[i];
    /*    [ 4 3]  coeff = 3150.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 3150.*Q[i];
    /*    [ 6 2]  coeff = 630.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 630.*Q[i];
    /*    [ 8 1]  coeff = 45.  */
    MorrisLecar_diff9(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += 45.*Q[i];
    /*    [ 10]  coeff = 1.  */
    MorrisLecar_diff10(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[10][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[11][i] = 0.0;
    /*    [ 0 0 0 0 0 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[10]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += Q[i];
    /*    [ 1 0 0 0 0 0 0 0 0 1]  coeff = 11.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[9]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 11.*Q[i];
    /*    [ 0 1 0 0 0 0 0 0 1]  coeff = 55.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 55.*Q[i];
    /*    [ 2 0 0 0 0 0 0 0 1]  coeff = 55.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 55.*Q[i];
    /*    [ 0 0 1 0 0 0 0 1]  coeff = 165.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 165.*Q[i];
    /*    [ 1 1 0 0 0 0 0 1]  coeff = 495.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 495.*Q[i];
    /*    [ 3 0 0 0 0 0 0 1]  coeff = 165.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 165.*Q[i];
    /*    [ 0 0 0 1 0 0 1]  coeff = 330.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[3],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 330.*Q[i];
    /*    [ 0 2 0 0 0 0 1]  coeff = 990.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 990.*Q[i];
    /*    [ 1 0 1 0 0 0 1]  coeff = 1320.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 1320.*Q[i];
    /*    [ 2 1 0 0 0 0 1]  coeff = 1980.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 1980.*Q[i];
    /*    [ 4 0 0 0 0 0 1]  coeff = 330.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 330.*Q[i];
    /*    [ 0 0 0 0 1 1]  coeff = 462.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[4],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 462.*Q[i];
    /*    [ 0 1 1 0 0 1]  coeff = 4620.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[2],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 4620.*Q[i];
    /*    [ 1 0 0 1 0 1]  coeff = 2310.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[3],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 2310.*Q[i];
    /*    [ 1 2 0 0 0 1]  coeff = 6930.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 2 0 1 0 0 1]  coeff = 4620.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 4620.*Q[i];
    /*    [ 3 1 0 0 0 1]  coeff = 4620.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 4620.*Q[i];
    /*    [ 5 0 0 0 0 1]  coeff = 462.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 462.*Q[i];
    /*    [ 0 0 2 0 1]  coeff = 4620.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[2],Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 4620.*Q[i];
    /*    [ 0 1 0 1 1]  coeff = 6930.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[3],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 0 3 0 0 1]  coeff = 6930.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 1 0 0 0 2]  coeff = 1386.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[4],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 1386.*Q[i];
    /*    [ 1 1 1 0 1]  coeff = 27720.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 27720.*Q[i];
    /*    [ 2 0 0 1 1]  coeff = 6930.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[3],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 2 2 0 0 1]  coeff = 20790.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 20790.*Q[i];
    /*    [ 3 0 1 0 1]  coeff = 9240.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 9240.*Q[i];
    /*    [ 4 1 0 0 1]  coeff = 6930.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 6 0 0 0 1]  coeff = 462.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[4]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 462.*Q[i];
    /*    [ 0 0 1 2]  coeff = 5775.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[2],Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 5775.*Q[i];
    /*    [ 0 2 1 1]  coeff = 34650.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 34650.*Q[i];
    /*    [ 1 0 2 1]  coeff = 23100.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 23100.*Q[i];
    /*    [ 1 1 0 2]  coeff = 17325.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 17325.*Q[i];
    /*    [ 1 3 0 1]  coeff = 34650.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 34650.*Q[i];
    /*    [ 2 1 1 1]  coeff = 69300.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 69300.*Q[i];
    /*    [ 3 0 0 2]  coeff = 5775.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[3],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 5775.*Q[i];
    /*    [ 3 2 0 1]  coeff = 34650.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 34650.*Q[i];
    /*    [ 4 0 1 1]  coeff = 11550.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 11550.*Q[i];
    /*    [ 5 1 0 1]  coeff = 6930.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 7 0 0 1]  coeff = 330.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[3]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 330.*Q[i];
    /*    [ 0 1 3]  coeff = 15400.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[2],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 15400.*Q[i];
    /*    [ 0 4 1]  coeff = 17325.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 17325.*Q[i];
    /*    [ 1 2 2]  coeff = 69300.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 69300.*Q[i];
    /*    [ 2 0 3]  coeff = 15400.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 15400.*Q[i];
    /*    [ 2 3 1]  coeff = 69300.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 69300.*Q[i];
    /*    [ 3 1 2]  coeff = 46200.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 46200.*Q[i];
    /*    [ 4 2 1]  coeff = 34650.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 34650.*Q[i];
    /*    [ 5 0 2]  coeff = 4620.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 4620.*Q[i];
    /*    [ 6 1 1]  coeff = 4620.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 4620.*Q[i];
    /*    [ 8 0 1]  coeff = 165.  */
    MorrisLecar_diff9(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[2]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 165.*Q[i];
    /*    [ 1 5]  coeff = 10395.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 10395.*Q[i];
    /*    [ 3 4]  coeff = 17325.  */
    MorrisLecar_diff7(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 17325.*Q[i];
    /*    [ 5 3]  coeff = 6930.  */
    MorrisLecar_diff8(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 6930.*Q[i];
    /*    [ 7 2]  coeff = 990.  */
    MorrisLecar_diff9(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 990.*Q[i];
    /*    [ 9 1]  coeff = 55.  */
    MorrisLecar_diff10(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += 55.*Q[i];
    /*    [ 11]  coeff = 1.  */
    MorrisLecar_diff11(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0]);
    for (i = 0; i < 2; ++i)
        Xderiv[11][i] += Q[i];

    for (i = 0; i < 2; ++i)
        Xderiv[12][i] = 0.0;
    /*    [ 0 0 0 0 0 0 0 0 0 0 0 1]  coeff = 1.  */
    MorrisLecar_diff1(Q,X,params,Xderiv[11]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += Q[i];
    /*    [ 1 0 0 0 0 0 0 0 0 0 1]  coeff = 12.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[0],Xderiv[10]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 12.*Q[i];
    /*    [ 0 1 0 0 0 0 0 0 0 1]  coeff = 66.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[1],Xderiv[9]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 66.*Q[i];
    /*    [ 2 0 0 0 0 0 0 0 0 1]  coeff = 66.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[9]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 66.*Q[i];
    /*    [ 0 0 1 0 0 0 0 0 1]  coeff = 220.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[2],Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 220.*Q[i];
    /*    [ 1 1 0 0 0 0 0 0 1]  coeff = 660.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 660.*Q[i];
    /*    [ 3 0 0 0 0 0 0 0 1]  coeff = 220.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[8]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 220.*Q[i];
    /*    [ 0 0 0 1 0 0 0 1]  coeff = 495.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[3],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 495.*Q[i];
    /*    [ 0 2 0 0 0 0 0 1]  coeff = 1485.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 1485.*Q[i];
    /*    [ 1 0 1 0 0 0 0 1]  coeff = 1980.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[2],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 1980.*Q[i];
    /*    [ 2 1 0 0 0 0 0 1]  coeff = 2970.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 2970.*Q[i];
    /*    [ 4 0 0 0 0 0 0 1]  coeff = 495.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[7]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 495.*Q[i];
    /*    [ 0 0 0 0 1 0 1]  coeff = 792.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[4],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 792.*Q[i];
    /*    [ 0 1 1 0 0 0 1]  coeff = 7920.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[2],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 7920.*Q[i];
    /*    [ 1 0 0 1 0 0 1]  coeff = 3960.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[3],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 3960.*Q[i];
    /*    [ 1 2 0 0 0 0 1]  coeff = 11880.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 11880.*Q[i];
    /*    [ 2 0 1 0 0 0 1]  coeff = 7920.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[2],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 7920.*Q[i];
    /*    [ 3 1 0 0 0 0 1]  coeff = 7920.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 7920.*Q[i];
    /*    [ 5 0 0 0 0 0 1]  coeff = 792.  */
    MorrisLecar_diff6(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[0],Xderiv[6]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 792.*Q[i];
    /*    [ 0 0 0 0 0 2]  coeff = 462.  */
    MorrisLecar_diff2(Q,X,params,Xderiv[5],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 462.*Q[i];
    /*    [ 0 0 2 0 0 1]  coeff = 9240.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[2],Xderiv[2],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 9240.*Q[i];
    /*    [ 0 1 0 1 0 1]  coeff = 13860.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[1],Xderiv[3],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 13860.*Q[i];
    /*    [ 0 3 0 0 0 1]  coeff = 13860.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[1],Xderiv[1],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 13860.*Q[i];
    /*    [ 1 0 0 0 1 1]  coeff = 5544.  */
    MorrisLecar_diff3(Q,X,params,Xderiv[0],Xderiv[4],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 5544.*Q[i];
    /*    [ 1 1 1 0 0 1]  coeff = 55440.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[1],Xderiv[2],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 55440.*Q[i];
    /*    [ 2 0 0 1 0 1]  coeff = 13860.  */
    MorrisLecar_diff4(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[3],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 13860.*Q[i];
    /*    [ 2 2 0 0 0 1]  coeff = 41580.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[1],Xderiv[1],Xderiv[5]);
    for (i = 0; i < 2; ++i)
        Xderiv[12][i] += 41580.*Q[i];
    /*    [ 3 0 1 0 0 1]  coeff = 18480.  */
    MorrisLecar_diff5(Q,X,params,Xderiv[0],Xderiv[0],Xderiv[