home *** CD-ROM | disk | FTP | other *** search
- /*
- ### compute a periodic orbit of period r for vector fields ###
-
- Input: r: 0: equilibrium point, >0: periodic orbits
- time: time
- time_step: (r=0) time_step for a next orbit
- x[n]: x(time)
- Output: x[n]: r=0: value of a vetor field, f(x)
- r>0: r'th image of an orbit x(time+period*r)
- */
-
- fode_user(x,time,r,dim)
- int r,dim;
- double x[],time;
- {
- int i,j;
- extern int model;
- extern double *t_v,*param;
- extern int (*f_p)();
-
-
- if(r==0){
- (int) f_p(t_v,0,x,param,time,dim);
- for(j=0;j<dim;j++) x[j] = t_v[j];
- }
- else {
- /*
- for(i=0;i<r;i++){
- fode_one_user(t_v,x);
- for(j=0;j<dim;j++) x[j] = t_v[j];
- }
- */
- }
- }
-