1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
#include <math.h>
/*
some example functions
*/
lv(double *in,double *out,int nin,int nout,double *var,double *con)
{
double x=in[0],y=in[1];
double a=in[2],b=in[3],c=in[4],d=in[5];
double t=in[6];
out[0]=a*x*(b-y);
out[1]=c*y*(-d+x);
}
vdp(double *in,double *out,int nin,int nout,double *var,double *con)
{
double x=in[0],y=in[1];
double a=in[2],b=in[3],c=in[4],d=in[5];
double t=in[6];
out[0]=y;
out[1]=-x+a*y*(1-x*x);
}
duff(double *in,double *out,int nin,int nout,double *var,double *con)
{
double x=in[0],y=in[1];
double a=in[2],b=in[3],c=in[4],d=in[5];
double t=in[6];
out[0]=y;
out[1]=x*(1-x*x)+a*sin(b*t)-c*y;
}
|
