home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / modellib / jjarray_def.c < prev    next >
Encoding:
Text File  |  1990-01-16  |  1.6 KB  |  87 lines
  1.  
  2. int jjarray_init()
  3. {
  4.     title_label = "Josephson Junction- Resistive Load";
  5.  
  6.     mapping_on = 0;
  7.     inverse_on = 1;
  8.     fderiv_on = 0;
  9.     enable_polar = 0;
  10.     enable_period = 1;
  11.  
  12.     var_dim = 6;
  13.     param_dim = 2;
  14.     func_dim = 2;
  15.  
  16.     (void) malloc_init();
  17.  
  18.     period_len[0] = twopi;
  19.     period_len[2] = twopi;
  20.     period_len[4] = twopi;
  21.     var_label[0] = "phi1";
  22.     var_label[1] = "psi1";
  23.     var_label[2] = "phi2";
  24.     var_label[3] = "psi2";
  25.     var_label[4] = "phi3";
  26.     var_label[5] = "psi3";
  27.     param_label[0] = "IB";
  28.     param_label[1] = "beta";
  29.     func_label[0] = "t";
  30.     func_label[1] = "Load";
  31.  
  32.     param[0] = 1.2;
  33.     param[1] = 1.6;
  34.  
  35.     var_i[0] = 0.3;
  36.     var_i[1] = 0.2;
  37.     var_i[2] = 0.324;
  38.     var_i[3] = 0.12;
  39.     var_i[4] = 0.41;
  40.     var_i[5] = 0;
  41.  
  42.     param_min[0]= -5; param_max[0]= 5;
  43.     param_min[1]= -5; param_max[1]= 5;
  44.  
  45.     var_min[0]= 0; var_max[0]= twopi;
  46.     var_min[1]= -.5; var_max[1]=2.5;
  47.     var_min[2]= 0; var_max[2]= twopi;
  48.     var_min[3]= -.5; var_max[3]=2.5;
  49.     var_min[4]= 0; var_max[4]= twopi;
  50.     var_min[5]= -.5; var_max[5]=2.5;
  51.     func_min[0]= -5; func_max[0]= 5;
  52.     func_min[1]= -3; func_max[0]= 6;
  53.  
  54.     f_p = jjarray_f;
  55.     func_p = jjarray_func;
  56. }
  57.  
  58. int jjarray_f(f,index,x,p,t,dim)
  59. int index,dim;
  60. double f[],x[],p[],t;
  61. {
  62.     int i;
  63.     double ib,beta,q;
  64.     ib = p[0];
  65.     beta = p[1];
  66.     q = 0;
  67.     for(i=1;i<dim;i += 2)
  68.         q += x[i];
  69.     f[0] = x[1];
  70.     f[1] = (ib - x[1] - sin(x[0]) - 2. / dim * q ) /beta;
  71.     f[2] = x[3];
  72.     f[3] = (ib - x[3] - sin(x[2]) - 2. / dim * q ) /beta;
  73.     f[4] = x[5];
  74.     f[5] = (ib - x[5] - sin(x[4]) - 2. / dim * q ) /beta;
  75. }
  76. int jjarray_func(f,x,p,t,dim)
  77. double f[],x[],p[],t;
  78. int dim;
  79. {
  80.     int i;
  81.  
  82.     f[0] = t;
  83.     f[1] = 0;
  84.     for(i=1;i<dim;i += 2)
  85.         f[1] += x[i];
  86. }
  87.