gsl 1.12.108

data/ext/geometry.c

@@ -0,0 +1,139 @@
+#include "rb_gsl.h"
+
+/*!
+  Counter-clockwise rotation around the X-axis 
+
+  / x' \    /  1     0      0    \ / x \
+  | y' |  = |  0  cos_th -sin_th | | y |
+  \ z' /    \  0  sin_th  cos_th / \ z /
+
+ */
+void vector3_rotateX(const double x[3] /*!< Input */, 
+                        double theta /*!< Rotation angle */, 
+                        double xout[3]  /*!< Output */)
+{
+  double a, b, c;
+  double costheta, sintheta;
+  
+  costheta = cos(theta);
+  sintheta = sin(theta);
+
+  a = x[0];
+  b = x[1]*costheta - x[2]*sintheta;
+  c = x[1]*sintheta + x[2]*costheta;
+
+  xout[0] = a;
+  xout[1] = b;
+  xout[2] = c;
+}
+
+/*!
+  Counter-clockwise rotation around the Y-axis 
+  (Note the sign of the matrix element.)
+
+  / x' \    /  cos_th     0   sin_th \ / x \
+  | y' |  = |    0        1      0   | | y |
+  \ z' /    \ -sin_th     0   cos_th / \ z /
+
+ */
+void vector3_rotateY(const double x[3] /*!< Input */, 
+                        double theta /*!< Rotation angle */,
+                        double xout[3] /*!< Output */)
+{
+  double a, b, c;
+  double costheta, sintheta;
+
+  costheta = cos(theta);
+  sintheta = sin(theta);
+
+  a =  x[0]*costheta + x[2]*sintheta;
+  b =  x[1];
+  c = -x[0]*sintheta + x[2]*costheta;
+
+  xout[0] = a;
+  xout[1] = b;
+  xout[2] = c;
+}
+
+/*!
+  Counter-clockwise rotation around the Z-axis 
+
+  / x' \    /  cos_th  -sin_th  0 \ / x \
+  | y' |  = |  sin_th   cos_th  0 | | y |
+  \ z' /    \    0        0     1 / \ z /
+
+ */
+void vector3_rotateZ(const double x[3] /*!< Input */, 
+                        double theta  /*!< Rotation angle */, 
+                        double xout[3] /*!< Output */)
+{
+  double a, b, c;
+  double costheta, sintheta;
+
+  costheta = cos(theta);
+  sintheta = sin(theta);
+
+  a = x[0]*costheta - x[1]*sintheta;
+  b = x[0]*sintheta + x[1]*costheta;
+  c = x[2];
+
+  xout[0] = a;
+  xout[1] = b;
+  xout[2] = c;
+}
+
+/*!
+  Rotate a 3-vector.
+
+  If flag != 1, cos/sin values stored in the static variables are used 
+  for efficiency.
+  The input and the output vectors (x, xout) can be the same pointer.
+
+  / x' \     /   cos_phi -sin_phi   0  \ /  cos_th   0   sin_th \  / x \
+  | y' |  =  |   sin_phi  cos_phi   0  | |     0     1          |  | y |
+  \ z' /     \      0        0      1  / \ -sin_th   0   cos_th /  \ z /
+
+ */
+void vector3_rotate(const double x[3], double theta, double phi,
+ double xout[3])
+{
+  vector3_rotateY(x, theta,xout);
+  vector3_rotateZ(xout, phi, xout);
+}
+
+static VALUE rb_gsl_vector_rotateX(VALUE obj, VALUE angle)
+{
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_vector, v);
+  vector3_rotateX(v->data, NUM2DBL(angle), v->data);
+  return obj;
+}
+static VALUE rb_gsl_vector_rotateY(VALUE obj, VALUE angle)
+{
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_vector, v);
+  vector3_rotateY(v->data, NUM2DBL(angle), v->data);
+  return obj;
+}
+static VALUE rb_gsl_vector_rotateZ(VALUE obj, VALUE angle)
+{
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_vector, v);
+  vector3_rotateZ(v->data, NUM2DBL(angle), v->data);
+  return obj;
+}
+static VALUE rb_gsl_vector_rotate(VALUE obj, VALUE theta, VALUE phi)
+{
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_vector, v);
+  vector3_rotate(v->data, NUM2DBL(theta), NUM2DBL(phi), v->data);
+  return obj;
+}
+
+void Init_geometry(VALUE module)
+{
+  rb_define_method(cgsl_vector, "rotateX", rb_gsl_vector_rotateX, 1);
+  rb_define_method(cgsl_vector, "rotateY", rb_gsl_vector_rotateY, 1);
+  rb_define_method(cgsl_vector, "rotateZ", rb_gsl_vector_rotateZ, 1);
+  rb_define_method(cgsl_vector, "rotate", rb_gsl_vector_rotate, 2);
+}

data/ext/graph.c

@@ -0,0 +1,1638 @@
+/*
+  graph.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY.
+*/
+#include "rb_gsl_config.h"
+#include "rb_gsl_graph.h"
+
+void make_graphcommand(char *command, VALUE hash)
+{
+  VALUE val;
+  if (TYPE(hash) == T_STRING) {
+    sprintf(command, "graph -T X -g 3 %s", STR2CSTR(hash));
+    return;
+  }
+    
+  strcpy(command, "graph");
+  if (TYPE(hash) != T_HASH) rb_raise(rb_eTypeError, 
+				     "wrong argument type %s (Hash expected)",
+				     rb_class2name(CLASS_OF(hash)));
+  if ((val = rb_hash_aref(hash, rb_str_new2("T"))) != Qnil)
+    sprintf(command, "%s -T %s", command, STR2CSTR(val));
+  else
+    sprintf(command, "%s -T X", command);
+
+  val = rb_hash_aref(hash, rb_str_new2("C"));
+  if (val == Qtrue)
+    sprintf(command, "%s -C", command);
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("g"))) != Qnil)
+    sprintf(command, "%s -g %d", command, (int) FIX2INT(val));
+  else
+    sprintf(command, "%s -g 3", command);
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("B"))) == Qtrue)
+    sprintf(command, "%s -B", command);
+  if ((val = rb_hash_aref(hash, rb_str_new2("E"))) != Qnil)
+    sprintf(command, "%s -E %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("f"))) != Qnil)
+    sprintf(command, "%s -f %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("F"))) != Qnil)
+    sprintf(command, "%s -F %s", command, STR2CSTR(val));
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("h"))) != Qnil)
+    sprintf(command, "%s -h %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("k"))) != Qnil)
+    sprintf(command, "%s -k %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("K"))) != Qnil)
+    sprintf(command, "%s -K %d", command, (int) FIX2INT(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("l"))) != Qnil) {
+    if (str_tail_grep(STR2CSTR(val), "xy") || str_tail_grep(STR2CSTR(val), "x/y"))
+      sprintf(command, "%s -l x -l y", command);
+    else
+      sprintf(command, "%s -l %s", command, STR2CSTR(val));
+  }
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("L"))) != Qnil)
+    sprintf(command, "%s -L \"%s\"", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("N"))) != Qnil)
+    sprintf(command, "%s -N %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("r"))) != Qnil)
+    sprintf(command, "%s -r %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("R"))) != Qnil)
+    sprintf(command, "%s -R %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("s"))) == Qtrue)
+    sprintf(command, "%s -s", command);
+  if ((val = rb_hash_aref(hash, rb_str_new2("t"))) == Qtrue)
+    sprintf(command, "%s -t", command);
+  if ((val = rb_hash_aref(hash, rb_str_new2("u"))) != Qnil)
+    sprintf(command, "%s -u %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("w"))) != Qnil)
+    sprintf(command, "%s -w %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("x"))) != Qnil)
+    sprintf(command, "%s -x %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("X"))) != Qnil)
+    sprintf(command, "%s -X \"%s\"", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("y"))) != Qnil)
+    sprintf(command, "%s -y %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("Y"))) != Qnil)
+    sprintf(command, "%s -Y \"%s\"", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("bg-color"))) != Qnil)
+    sprintf(command, "%s --bg-color %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("bitmap-size"))) != Qnil)
+    sprintf(command, "%s --bitmap-size %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("frame-color"))) != Qnil)
+    sprintf(command, "%s --frame-color %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("frame-line-width"))) != Qnil)
+    sprintf(command, "%s --frame-line-width %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("max-line-length"))) != Qnil)
+    sprintf(command, "%s --max-line-length %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("page-size"))) != Qnil)
+    sprintf(command, "%s --page-size %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("pen-colors"))) != Qnil)
+    sprintf(command, "%s --pen-colors %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("rotation"))) != Qnil)
+    sprintf(command, "%s --rotation %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("title-font-name"))) != Qnil)
+    sprintf(command, "%s --title-font-name %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("title-font-size"))) != Qnil)
+    sprintf(command, "%s --title-font-size %f", command, NUM2DBL(val));
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("toggle-rotate-y-label"))) == Qtrue)
+    sprintf(command, "%s --toggle-rotate-y-label", command);
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("m"))) != Qnil)
+    sprintf(command, "%s -m %d", command, (int) FIX2INT(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("S"))) != Qnil)
+    sprintf(command, "%s -S %d", command, (int) FIX2INT(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("W"))) != Qnil)
+    sprintf(command, "%s -W %f", command, NUM2DBL(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("q"))) != Qnil)
+    sprintf(command, "%s -q %f", command, NUM2DBL(val));
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("symbol-font-name"))) != Qnil)
+    sprintf(command, "%s --symbol-font-name %s", command, STR2CSTR(val));
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("reposition"))) != Qnil)
+    sprintf(command, "%s --reposition %s", command, STR2CSTR(val));
+  if ((val = rb_hash_aref(hash, rb_str_new2("blankout"))) != Qnil)
+    sprintf(command, "%s --blankout %s", command, STR2CSTR(val));
+
+  if ((val = rb_hash_aref(hash, rb_str_new2("O"))) == Qtrue)
+    sprintf(command, "%s -O", command);
+}
+
+static void gsl_graph_init(gsl_graph *g);
+gsl_graph* gsl_graph_new()
+{
+  gsl_graph *g = NULL;
+  g = ALLOC(gsl_graph);
+  gsl_graph_init(g);
+  return g;
+}
+
+static void gsl_graph_init(gsl_graph *g)
+{
+  g->xdata = Qnil;
+  g->ydata = Qnil;
+
+  g->T = Qnil;
+  g->E = Qnil;
+  g->f = Qnil;
+  g->F = Qnil;
+  g->g = Qnil;
+  g->h = Qnil;
+  g->k = Qnil;
+  g->K = Qnil;
+  g->l = Qnil;
+  g->L = Qnil;
+  g->N = Qnil;
+  g->r = Qnil;
+  g->R = Qnil;
+  g->u = Qnil;
+  g->w = Qnil;
+  g->x = Qnil;
+  g->y = Qnil;
+  g->X = Qnil;
+  g->Y = Qnil;
+
+  g->bg = Qnil;
+  g->bitmap_size = Qnil;
+  g->frame = Qnil;
+  g->frame_line_width = Qnil;
+  g->max_line_length = Qnil;
+  g->page_size = Qnil;
+  g->pen_colors = Qnil;
+  g->rotation = Qnil;
+  g->title_font_size = Qnil;
+  g->title_font_name = Qnil;
+  g->rotate_y_label = Qfalse;
+
+  g->I = Qnil;
+  g->m = Qnil;
+  g->S = Qnil;
+  g->W = Qnil;
+  g->q = Qnil;
+  g->symbol_font_name = Qnil;
+  g->reposition = Qnil;
+  g->blankout = Qnil;
+
+  g->s = Qfalse;
+  g->t = Qfalse;
+
+  g->B = Qfalse;
+  g->C = Qfalse;
+  g->O = Qfalse;
+}
+
+static void gsl_graph_mark(gsl_graph *g)
+{
+  rb_gc_mark(g->xdata);
+  rb_gc_mark(g->ydata);
+  rb_gc_mark(g->T);
+  rb_gc_mark(g->E);
+  rb_gc_mark(g->f);
+  rb_gc_mark(g->F);
+  rb_gc_mark(g->g);
+  rb_gc_mark(g->h);
+  rb_gc_mark(g->k);
+  rb_gc_mark(g->K);
+  rb_gc_mark(g->l);
+  rb_gc_mark(g->L);
+  rb_gc_mark(g->N);
+  rb_gc_mark(g->r);
+  rb_gc_mark(g->R);
+  rb_gc_mark(g->s);
+  rb_gc_mark(g->t);
+  rb_gc_mark(g->u);
+  rb_gc_mark(g->w);
+  rb_gc_mark(g->x);
+  rb_gc_mark(g->y);
+  rb_gc_mark(g->X);
+  rb_gc_mark(g->Y);
+  rb_gc_mark(g->bg);
+  rb_gc_mark(g->bitmap_size);
+  rb_gc_mark(g->frame);
+  rb_gc_mark(g->frame_line_width);
+  rb_gc_mark(g->max_line_length);
+  rb_gc_mark(g->page_size);
+  rb_gc_mark(g->pen_colors);
+  rb_gc_mark(g->rotation);
+  rb_gc_mark(g->title_font_name);
+  rb_gc_mark(g->title_font_size);
+  rb_gc_mark(g->rotate_y_label);
+  rb_gc_mark(g->I);
+  rb_gc_mark(g->B);
+  rb_gc_mark(g->m);
+  rb_gc_mark(g->S);
+  rb_gc_mark(g->W);
+  rb_gc_mark(g->q);
+  rb_gc_mark(g->C);
+  rb_gc_mark(g->symbol_font_name);
+  rb_gc_mark(g->reposition);
+  rb_gc_mark(g->blankout);
+  rb_gc_mark(g->O);
+}
+
+void gsl_graph_free(gsl_graph *g)
+{
+  free((gsl_graph *) g);
+}
+
+static VALUE rb_gsl_graph_set_xdata(VALUE obj, VALUE xx);
+static VALUE rb_gsl_graph_set_ydata(VALUE obj, VALUE yy);
+static VALUE rb_gsl_graph_set_xydata(VALUE obj, VALUE xx, VALUE yy);
+static VALUE rb_gsl_graph_new(int argc, VALUE *argv, VALUE klass)
+{
+  gsl_graph *g = NULL;
+  VALUE obj;
+  g = gsl_graph_new();
+  obj = Data_Wrap_Struct(klass, gsl_graph_mark, gsl_graph_free, g);
+  switch (argc) {
+  case 1:
+    rb_gsl_graph_set_xdata(obj, argv[0]);
+    break;
+  case 2:
+    rb_gsl_graph_set_xydata(obj, argv[0], argv[1]);
+    break;
+  }
+  return obj;
+}
+
+static VALUE rb_gsl_graph_init(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  gsl_graph_init(g);
+  return obj;
+}
+
+static VALUE rb_gsl_graph_xdata(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->xdata;
+}
+
+static VALUE rb_gsl_graph_ydata(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->ydata;
+}
+
+static VALUE rb_gsl_graph_xydata(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return rb_ary_new3(2, g->xdata, g->ydata);
+}
+
+static VALUE rb_gsl_graph_set_xdata(VALUE obj, VALUE xx)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  CHECK_VECTOR(xx);
+  g->xdata = xx;
+  return obj;
+}
+
+static VALUE rb_gsl_graph_set_ydata(VALUE obj, VALUE yy)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  CHECK_VECTOR(yy);
+  g->ydata = yy;
+  return obj;
+}
+
+static VALUE rb_gsl_graph_set_xydata(VALUE obj, VALUE xx, VALUE yy)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  CHECK_VECTOR(xx); CHECK_VECTOR(yy);
+  g->xdata = xx;
+  g->ydata = yy;
+  return obj;
+}
+
+static VALUE rb_gsl_graph_set_T(VALUE obj, VALUE T)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  Check_Type(T, T_STRING);
+  g->T = T;
+  return T;
+}
+
+static VALUE rb_gsl_graph_T(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->T;
+}
+
+static VALUE rb_gsl_graph_set_E(VALUE obj, VALUE E)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  Check_Type(E, T_STRING);
+  g->E = E;
+  return E;
+}
+
+static VALUE rb_gsl_graph_E(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->E;
+}
+
+static VALUE rb_gsl_graph_set_f(VALUE obj, VALUE f)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->f = f;
+  return f;
+}
+
+static VALUE rb_gsl_graph_f(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->f;
+}
+
+static VALUE rb_gsl_graph_set_F(VALUE obj, VALUE F)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->F = F;
+  return F;
+}
+
+static VALUE rb_gsl_graph_F(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->F;
+}
+
+static VALUE rb_gsl_graph_set_g(VALUE obj, VALUE gg)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  Check_Type(gg, T_FIXNUM);
+  g->g = gg;
+  return gg;
+}
+
+static VALUE rb_gsl_graph_g(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->g;
+}
+
+static VALUE rb_gsl_graph_set_h(VALUE obj, VALUE h)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->h = h;
+  return h;
+}
+
+static VALUE rb_gsl_graph_h(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->h;
+}
+
+static VALUE rb_gsl_graph_set_k(VALUE obj, VALUE k)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->k = k;
+  return k;
+}
+
+static VALUE rb_gsl_graph_k(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->k;
+}
+
+static VALUE rb_gsl_graph_set_K(VALUE obj, VALUE K)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->K = K;
+  return K;
+}
+
+static VALUE rb_gsl_graph_K(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->K;
+}
+
+static VALUE rb_gsl_graph_set_l(VALUE obj, VALUE l)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->l = l;
+  return l;
+}
+
+static VALUE rb_gsl_graph_l(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->l;
+}
+
+static VALUE rb_gsl_graph_set_L(VALUE obj, VALUE L)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->L = L;
+  return L;
+}
+
+static VALUE rb_gsl_graph_L(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->L;
+}
+
+static VALUE rb_gsl_graph_set_N(VALUE obj, VALUE N)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->N = N;
+  return N;
+}
+
+static VALUE rb_gsl_graph_N(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->N;
+}
+
+static VALUE rb_gsl_graph_set_r(VALUE obj, VALUE r)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->r = r;
+  return r;
+}
+
+static VALUE rb_gsl_graph_r(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->r;
+}
+
+static VALUE rb_gsl_graph_set_R(VALUE obj, VALUE R)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->R = R;
+  return R;
+}
+
+static VALUE rb_gsl_graph_R(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->R;
+}
+
+static VALUE rb_gsl_graph_set_s(VALUE obj, VALUE s)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->s = s;
+  return s;
+}
+
+static VALUE rb_gsl_graph_s(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->s;
+}
+
+static VALUE rb_gsl_graph_set_t(VALUE obj, VALUE t)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->t = t;
+  return t;
+}
+
+static VALUE rb_gsl_graph_t(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->t;
+}
+
+static VALUE rb_gsl_graph_set_u(VALUE obj, VALUE u)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->u = u;
+  return u;
+}
+
+static VALUE rb_gsl_graph_u(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->u;
+}
+
+static VALUE rb_gsl_graph_set_w(VALUE obj, VALUE w)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->w = w;
+  return w;
+}
+
+static VALUE rb_gsl_graph_w(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->w;
+}
+
+static VALUE rb_gsl_graph_set_x(VALUE obj, VALUE x)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->x = x;
+  return x;
+}
+
+static VALUE rb_gsl_graph_x(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->x;
+}
+
+static VALUE rb_gsl_graph_set_y(VALUE obj, VALUE y)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->y = y;
+  return y;
+}
+
+static VALUE rb_gsl_graph_y(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->y;
+}
+static VALUE rb_gsl_graph_set_X(VALUE obj, VALUE X)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  Check_Type(X, T_STRING);
+  g->X = X;
+  return X;
+}
+
+static VALUE rb_gsl_graph_set_Y(VALUE obj, VALUE Y)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  Check_Type(Y, T_STRING);
+  g->Y = Y;
+  return Y;
+}
+
+static VALUE rb_gsl_graph_X(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->X;
+}
+
+static VALUE rb_gsl_graph_Y(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->Y;
+}
+
+static VALUE rb_gsl_graph_set_bg(VALUE obj, VALUE bg)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->bg = bg;
+  return bg;
+}
+
+static VALUE rb_gsl_graph_bg(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->bg;
+}
+
+static VALUE rb_gsl_graph_set_bitmap_size(VALUE obj, VALUE bitmap_size)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->bitmap_size = bitmap_size;
+  return bitmap_size;
+}
+
+static VALUE rb_gsl_graph_bitmap_size(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->bitmap_size;
+}
+
+static VALUE rb_gsl_graph_set_frame(VALUE obj, VALUE frame)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->frame = frame;
+  return frame;
+}
+
+static VALUE rb_gsl_graph_frame(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->frame;
+}
+
+static VALUE rb_gsl_graph_set_frame_line_width(VALUE obj, VALUE frame_line_width)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->frame_line_width = frame_line_width;
+  return frame_line_width;
+}
+
+static VALUE rb_gsl_graph_frame_line_width(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->frame_line_width;
+}
+
+static VALUE rb_gsl_graph_set_max_line_length(VALUE obj, VALUE max_line_length)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->max_line_length = max_line_length;
+  return max_line_length;
+}
+
+static VALUE rb_gsl_graph_max_line_length(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->max_line_length;
+}
+
+static VALUE rb_gsl_graph_set_page_size(VALUE obj, VALUE page_size)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->page_size = page_size;
+  return page_size;
+}
+
+static VALUE rb_gsl_graph_page_size(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->page_size;
+}
+
+static VALUE rb_gsl_graph_set_pen_colors(VALUE obj, VALUE pen_colors)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->pen_colors = pen_colors;
+  return pen_colors;
+}
+
+static VALUE rb_gsl_graph_pen_colors(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->pen_colors;
+}
+
+static VALUE rb_gsl_graph_set_rotation(VALUE obj, VALUE rotation)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->rotation = rotation;
+  return rotation;
+}
+
+static VALUE rb_gsl_graph_rotation(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->rotation;
+}
+
+static VALUE rb_gsl_graph_set_title_font_name(VALUE obj, VALUE title_font_name)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->title_font_name = title_font_name;
+  return title_font_name;
+}
+
+static VALUE rb_gsl_graph_title_font_name(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->title_font_name;
+}
+
+static VALUE rb_gsl_graph_set_title_font_size(VALUE obj, VALUE title_font_size)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->title_font_size = title_font_size;
+  return title_font_size;
+}
+
+static VALUE rb_gsl_graph_title_font_size(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->title_font_size;
+}
+
+static VALUE rb_gsl_graph_set_rotate_y_label(VALUE obj, VALUE rotate_y_label)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->rotate_y_label = rotate_y_label;
+  return rotate_y_label;
+}
+
+static VALUE rb_gsl_graph_rotate_y_label(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->rotate_y_label;
+}
+
+static VALUE rb_gsl_graph_set_I(VALUE obj, VALUE I)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->I = I;
+  return I;
+}
+
+static VALUE rb_gsl_graph_I(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->I;
+}
+
+static VALUE rb_gsl_graph_set_B(VALUE obj, VALUE B)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->B = B;
+  return B;
+}
+
+static VALUE rb_gsl_graph_B(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->B;
+}
+
+static VALUE rb_gsl_graph_set_m(VALUE obj, VALUE m)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->m = m;
+  return m;
+}
+
+static VALUE rb_gsl_graph_m(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->m;
+}
+
+static VALUE rb_gsl_graph_set_S(VALUE obj, VALUE S)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->S = S;
+  return S;
+}
+
+static VALUE rb_gsl_graph_S(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->S;
+}
+
+static VALUE rb_gsl_graph_set_W(VALUE obj, VALUE W)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->W = W;
+  return W;
+}
+
+static VALUE rb_gsl_graph_W(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->W;
+}
+
+static VALUE rb_gsl_graph_set_q(VALUE obj, VALUE q)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->q = q;
+  return q;
+}
+
+static VALUE rb_gsl_graph_q(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->q;
+}
+
+static VALUE rb_gsl_graph_set_C(VALUE obj, VALUE C)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->C = C;
+  return C;
+}
+
+static VALUE rb_gsl_graph_C(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->C;
+}
+
+static VALUE rb_gsl_graph_set_symbol_font_name(VALUE obj, VALUE symbol_font_name)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->symbol_font_name = symbol_font_name;
+  return symbol_font_name;
+}
+
+static VALUE rb_gsl_graph_symbol_font_name(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->symbol_font_name;
+}
+
+static VALUE rb_gsl_graph_set_reposition(VALUE obj, VALUE r)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->reposition = r;
+  return r;
+}
+
+static VALUE rb_gsl_graph_reposition(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->reposition;
+}
+
+static VALUE rb_gsl_graph_set_blankout(VALUE obj, VALUE r)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->blankout = r;
+  return r;
+}
+
+static VALUE rb_gsl_graph_blankout(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->blankout;
+}
+
+static VALUE rb_gsl_graph_set_O(VALUE obj, VALUE O)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  g->O = O;
+  return O;
+}
+
+static VALUE rb_gsl_graph_O(VALUE obj)
+{
+  gsl_graph *g = NULL;
+  Data_Get_Struct(obj, gsl_graph, g);
+  return g->O;
+}
+
+static void gsl_graph_set_command(gsl_graph *g, char *command)
+{
+  char str[256];
+  size_t i, len;
+  VALUE val;
+  strcpy(command, "graph");
+  if (g->T == Qnil)
+    sprintf(command, "%s -T X", command);
+  else
+    sprintf(command, "%s -T %s", command, STR2CSTR(g->T));
+
+  if (TYPE(g->E) == T_STRING) {
+    strcpy(str, STR2CSTR(g->E));
+    if (strcmp(str, "x") == 0)
+      sprintf(command, "%s -E x", command);
+    else if (strcmp(str, "y") == 0)
+      sprintf(command, "%s -E y", command);
+    else if (strcmp(str, "xy")*strcmp(str, "x/y") == 0)
+      sprintf(command, "%s -E x -E y", command);
+    else 
+      rb_raise(rb_eRuntimeError, "unrecognized -E option %s", str);
+  }
+
+  if (g->f != Qnil)
+    sprintf(command, "%s -f %f", command, NUM2DBL(g->f));
+
+  if (TYPE(g->F) == T_STRING)
+    sprintf(command, "%s -F %s", command, STR2CSTR(g->F));
+
+  if (TYPE(g->g) == T_FIXNUM)
+    sprintf(command, "%s -g %d", command, (int) FIX2INT(g->g));
+
+  if (g->h != Qnil)
+    sprintf(command, "%s -h %f", command, NUM2DBL(g->h));
+
+  if (g->k != Qnil)
+    sprintf(command, "%s -k %f", command, NUM2DBL(g->k));
+
+  if (TYPE(g->K) == T_FIXNUM)
+    sprintf(command, "%s -K %d", command, (int) FIX2INT(g->K));
+
+  if (TYPE(g->l) == T_STRING) {
+    strcpy(str, STR2CSTR(g->l));
+    if (strcmp(str, "x") == 0)
+      sprintf(command, "%s -l x", command);
+    else if (strcmp(str, "y") == 0)
+      sprintf(command, "%s -l y", command);
+    else if (strcmp(str, "xy")*strcmp(str, "x/y") == 0)
+      sprintf(command, "%s -l x -l y", command);
+    else 
+      rb_raise(rb_eRuntimeError, "unrecognized -l option %s", str);
+  }
+
+  if (TYPE(g->L) == T_STRING)
+    sprintf(command, "%s -L \"%s\"", command, STR2CSTR(g->L));
+
+  if (TYPE(g->N) == T_STRING) {
+    strcpy(str, STR2CSTR(g->N));
+    if (strcmp(str, "x") == 0)
+      sprintf(command, "%s -N x", command);
+    else if (strcmp(str, "y") == 0)
+      sprintf(command, "%s -N y", command);
+    else if (strcmp(str, "xy")*strcmp(str, "x/y") == 0)
+      sprintf(command, "%s -N x -N y", command);
+    else 
+      rb_raise(rb_eRuntimeError, "unrecognized -N option %s", str);
+  }
+
+  if (g->r != Qnil)
+    sprintf(command, "%s -r %f", command, NUM2DBL(g->r));
+
+  if (TYPE(g->R) == T_STRING) {
+    strcpy(str, STR2CSTR(g->R));
+    if (strcmp(str, "x") == 0)
+      sprintf(command, "%s -R x", command);
+    else if (strcmp(str, "y") == 0)
+      sprintf(command, "%s -R y", command);
+    else if (strcmp(str, "xy")*strcmp(str, "x/y") == 0)
+      sprintf(command, "%s -R x -R y", command);
+    else 
+      rb_raise(rb_eRuntimeError, "unrecognized -R option %s", str);
+  }
+
+  if (g->u != Qnil)
+    sprintf(command, "%s -u %f", command, NUM2DBL(g->u));
+
+  if (g->w != Qnil)
+    sprintf(command, "%s -w %f", command, NUM2DBL(g->w));
+
+  switch (TYPE(g->x)) {
+  case T_STRING:
+    sprintf(command, "%s -x %s", command, STR2CSTR(g->x));
+    break;
+  case T_ARRAY:
+    sprintf(command, "%s -x", command);
+    //    len = RARRAY(g->x)->len;
+    len = RARRAY_LEN(g->x);
+    for (i = 0; i < len; i++) {
+      val = rb_ary_entry(g->x, i);
+      Need_Float(val);
+      sprintf(command, "%s %f", command, NUM2DBL(val));
+    }
+    break;
+  default:
+    /* do nothing */
+    break;
+  }
+
+  switch (TYPE(g->y)) {
+  case T_STRING:
+    sprintf(command, "%s -y %s", command, STR2CSTR(g->y));
+    break;
+  case T_ARRAY:
+    sprintf(command, "%s -y", command);
+    //    len = RARRAY(g->y)->len;
+    len = RARRAY_LEN(g->y);
+    for (i = 0; i < len; i++) {
+      val = rb_ary_entry(g->y, i);
+      Need_Float(val);
+      sprintf(command, "%s %f", command, NUM2DBL(val));
+    }
+    break;
+  default:
+    /* do nothing */
+    break;
+  }
+
+  if (g->X != Qnil)
+    sprintf(command, "%s -X \"%s\"", command, STR2CSTR(g->X));
+  if (g->Y != Qnil)
+    sprintf(command, "%s -Y \"%s\"", command, STR2CSTR(g->Y));
+
+  if (TYPE(g->bg) == T_STRING)
+    sprintf(command, "%s --bg-color %s", command, STR2CSTR(g->bg));
+
+  if (TYPE(g->bitmap_size) == T_STRING)
+    sprintf(command, "%s --bitmap-size %s", command, STR2CSTR(g->bitmap_size));
+
+  if (TYPE(g->frame) == T_STRING)
+    sprintf(command, "%s --frame-color %s", command, STR2CSTR(g->frame));
+
+  if (g->frame_line_width != Qnil)
+    sprintf(command, "%s --frame-line-width %f", command, NUM2DBL(g->frame_line_width));
+
+ if (g->max_line_length != Qnil)
+   sprintf(command, "%s --max_line_length %d", command, 
+	   (int) FIX2INT(g->max_line_length));
+ 
+  if (g->page_size != Qnil)
+    sprintf(command, "%s --page-size %s", command, STR2CSTR(g->page_size));
+
+  if (g->pen_colors != Qnil)
+    sprintf(command, "%s --pen-colors %s", command, STR2CSTR(g->pen_colors));
+
+  if (g->rotation != Qnil)
+    sprintf(command, "%s --rotation %d", command, (int) FIX2INT(g->rotation));
+
+  if (g->title_font_name != Qnil)
+    sprintf(command, "%s --title-font-name %s", command, STR2CSTR(g->title_font_name));  if (g->title_font_size != Qnil)
+    sprintf(command, "%s --title-font-size %f", command, NUM2DBL(g->title_font_size));
+
+  if (g->rotate_y_label == Qtrue)
+    sprintf(command, "%s --toggle-rotate-y-label", command);
+
+  if (g->I != Qnil)
+    sprintf(command, "%s -I %s", command, STR2CSTR(g->I));
+
+  if (g->s == Qtrue)
+    sprintf(command, "%s -s", command);
+  if (g->t == Qtrue)
+    sprintf(command, "%s -t", command);
+  if (g->B == Qtrue)
+    sprintf(command, "%s -B", command);
+
+  if (g->m != Qnil)
+    sprintf(command, "%s -m %d", command, (int) FIX2INT(g->m));
+
+  switch (TYPE(g->S)) {
+  case T_STRING:
+    sprintf(command, "%s -S %s", command, STR2CSTR(g->S));
+    break;
+  case T_ARRAY:
+    //    if (RARRAY(g->S)->len == 2) 
+    if (RARRAY_LEN(g->S) == 2) 
+      sprintf(command, "%s -S %d %f", command, (int) FIX2INT(rb_ary_entry(g->S, 0)),
+	      NUM2DBL(rb_ary_entry(g->S, 1)));
+    break;
+  default:
+    /* do nothing */
+    break;
+  }
+
+  if (g->W != Qnil)
+    sprintf(command, "%s -W %f", command, NUM2DBL(g->W));
+
+  if (g->q != Qnil)
+    sprintf(command, "%s -q %f", command, NUM2DBL(g->q));
+
+  if (g->C == Qtrue)
+    sprintf(command, "%s -C", command);
+
+  if (g->symbol_font_name != Qnil)
+    sprintf(command, "%s --symbol_font_name %s", command, STR2CSTR(g->symbol_font_name));
+
+  switch (TYPE(g->reposition)) {
+  case T_STRING:
+    sprintf(command, "%s --reposition %s", command, STR2CSTR(g->reposition));
+    break;
+  case T_ARRAY:
+    sprintf(command, "%s --reposition", command);
+    //    len =  RARRAY(g->reposition)->len;
+    len =  RARRAY_LEN(g->reposition);
+    for (i = 0; i <len; i++) {
+      val = rb_ary_entry(g->reposition, i);
+      Need_Float(val);
+      sprintf(command, "%s %f", command, NUM2DBL(val));
+    }
+    break;
+  default:
+    /* do nothing */
+    break;
+  }
+
+  if (g->blankout != Qnil)
+    sprintf(command, "%s --blankout %f", command, NUM2DBL(g->blankout));
+
+  if (g->O == Qtrue)
+    sprintf(command, "%s -O", command);
+}
+
+static VALUE rb_gsl_graph_graph(int argc, VALUE *argv, VALUE obj)
+{
+#ifdef HAVE_GNU_GRAPH
+  gsl_graph *g = NULL;
+  gsl_histogram *h = NULL;
+  gsl_vector *x = NULL, *y = NULL;
+  size_t i, size;
+  FILE *fp;
+  char command[1024];
+  Data_Get_Struct(obj, gsl_graph, g);
+
+  gsl_graph_set_command(g, command);
+  switch (argc) {
+  case 3:
+    Check_Type(argv[2], T_STRING);
+    sprintf(command, "%s %s", command, STR2CSTR(argv[2]));
+    /* no break */
+  case 2:
+    if (TYPE(argv[1]) == T_STRING) {
+      sprintf(command, "%s %s", command, STR2CSTR(argv[1]));
+    } else if (VECTOR_P(argv[1])) {
+      g->ydata = argv[1];
+    } else {
+      rb_raise(rb_eTypeError, "wrong argument type %s (Vector or String expected)",
+	       rb_class2name(CLASS_OF(argv[1])));
+    }
+    /* no break */
+  case 1:
+    if (TYPE(argv[0]) == T_STRING) {
+      sprintf(command, "%s %s", command, STR2CSTR(argv[0]));
+    } else if (VECTOR_P(argv[0])) {
+      g->xdata = argv[0]; 
+    } else if (HISTOGRAM_P(argv[0])) {
+      Data_Get_Struct(argv[0], gsl_histogram, h);
+    } else {
+      rb_raise(rb_eTypeError, "wrong argument type %s (Vector or String expected)",
+	       rb_class2name(CLASS_OF(argv[0])));
+    }
+    break;
+  default:
+    rb_raise(rb_eArgError, "wrong number of argumeuts (%d for 1-3)", argc);
+    break;
+  }
+  if (VECTOR_P(g->xdata)) Data_Get_Struct(g->xdata, gsl_vector, x);
+  if (VECTOR_P(g->ydata)) Data_Get_Struct(g->ydata, gsl_vector, y);
+
+  if (x == NULL && h == NULL)
+    rb_raise(rb_eRuntimeError, "data is not given");
+
+  if (h) size = h->n;
+  else size = x->size;
+
+  fp = popen(command, "w");
+  if (fp == NULL) rb_raise(rb_eIOError, "GNU graph not found.");
+  for (i = 0; i < size; i++) {
+    if (h) 
+      fprintf(fp, "%g %g\n%g %g\n", h->range[i], h->bin[i], h->range[i+1], h->bin[i]);
+    else if (y == NULL) 
+      fprintf(fp, "%d %g\n", (int) i, gsl_vector_get(x, i));
+    else
+      fprintf(fp, "%g %g\n", gsl_vector_get(x, i), gsl_vector_get(y, i));
+  }
+  fflush(fp);
+  pclose(fp);
+  fp = NULL;
+  return Qtrue;
+#else
+  rb_raise(rb_eNoMethodError, "GNU plotutils required");
+  return Qfalse;
+#endif
+}
+
+static VALUE rb_gsl_graph_step(int argc, VALUE *argv, VALUE obj)
+{
+#ifdef HAVE_GNU_GRAPH
+  gsl_graph *g = NULL;
+  gsl_vector *x = NULL, *y = NULL;
+  size_t i, size;
+  FILE *fp;
+  char command[1024];
+  Data_Get_Struct(obj, gsl_graph, g);
+
+  gsl_graph_set_command(g, command);
+  switch (argc) {
+  case 3:
+    Check_Type(argv[2], T_STRING);
+    sprintf(command, "%s %s", command, STR2CSTR(argv[2]));
+    /* no break */
+  case 2:
+    if (TYPE(argv[1]) == T_STRING) {
+      sprintf(command, "%s %s", command, STR2CSTR(argv[1]));
+    } else if (VECTOR_P(argv[1])) {
+      g->ydata = argv[1];
+    } else {
+      rb_raise(rb_eTypeError, "wrong argument type %s (Vector or String expected)",
+	       rb_class2name(CLASS_OF(argv[1])));
+    }
+    /* no break */
+  case 1:
+    if (TYPE(argv[0]) == T_STRING) {
+      sprintf(command, "%s %s", command, STR2CSTR(argv[0]));
+    } else if (VECTOR_P(argv[0])) {
+      g->xdata = argv[0]; 
+    } else {
+      rb_raise(rb_eTypeError, "wrong argument type %s (Vector or String expected)",
+	       rb_class2name(CLASS_OF(argv[0])));
+    }
+    break;
+  default:
+    rb_raise(rb_eArgError, "wrong number of argumeuts (%d for 1-3)", argc);
+    break;
+  }
+  if (VECTOR_P(g->xdata)) Data_Get_Struct(g->xdata, gsl_vector, x);
+  if (VECTOR_P(g->ydata)) Data_Get_Struct(g->ydata, gsl_vector, y);
+
+  if (x == NULL)
+    rb_raise(rb_eRuntimeError, "data is not given");
+
+  size = x->size;
+
+  fp = popen(command, "w");
+  if (fp == NULL) rb_raise(rb_eIOError, "GNU graph not found.");
+  for (i = 0; i < size; i++) {
+    if (y == NULL) {
+      fprintf(fp, "%d %g\n%d %g\n", (int) i, gsl_vector_get(x, i),
+	      (int) (i+1), gsl_vector_get(x, i));
+    } else {
+      if (i != size-1)
+	fprintf(fp, "%g %g\n%g %g\n", gsl_vector_get(x, i), gsl_vector_get(y, i),
+		gsl_vector_get(x, i+1), gsl_vector_get(y, i));
+      else 	
+	fprintf(fp, "%g %g\n%g %g", gsl_vector_get(x, i), gsl_vector_get(y, i),
+		2.0*gsl_vector_get(x, i)-gsl_vector_get(x, i-1), gsl_vector_get(y, i));
+    }
+  }
+  fflush(fp);
+  pclose(fp);
+  fp = NULL;
+  return Qtrue;
+#else
+  rb_raise(rb_eNoMethodError, "GNU plotutils required");
+  return Qfalse;
+#endif
+}
+
+void Init_gsl_graph(VALUE module)
+{
+  VALUE cgsl_graph;
+
+  cgsl_graph = rb_define_class_under(module, "Graph", cGSL_Object);
+
+  rb_define_singleton_method(cgsl_graph, "new", rb_gsl_graph_new, -1);
+  rb_define_singleton_method(cgsl_graph, "alloc", rb_gsl_graph_new, -1);
+
+  /*****/
+
+  rb_define_method(cgsl_graph, "init", rb_gsl_graph_init, 0);
+
+  rb_define_method(cgsl_graph, "set_xdata", rb_gsl_graph_set_xdata, 1);
+  rb_define_method(cgsl_graph, "set_ydata", rb_gsl_graph_set_ydata, 1);
+  rb_define_method(cgsl_graph, "set_xydata", rb_gsl_graph_set_xydata, 2);
+
+  rb_define_method(cgsl_graph, "xdata", rb_gsl_graph_xdata, 0);
+  rb_define_method(cgsl_graph, "ydata", rb_gsl_graph_ydata, 0);
+  rb_define_method(cgsl_graph, "xydata", rb_gsl_graph_xydata, 0);
+
+  rb_define_method(cgsl_graph, "graph", rb_gsl_graph_graph, -1);
+  rb_define_alias(cgsl_graph, "draw", "graph");
+  rb_define_alias(cgsl_graph, "plot", "graph");
+
+  rb_define_method(cgsl_graph, "graph_step", rb_gsl_graph_step, -1);
+  rb_define_alias(cgsl_graph, "step", "graph_step");
+
+  /*****/
+  rb_define_method(cgsl_graph, "set_T", rb_gsl_graph_set_T, 1);
+  rb_define_alias(cgsl_graph, "T=", "set_T");
+  rb_define_alias(cgsl_graph, "display_type=", "set_T");
+  rb_define_method(cgsl_graph, "T", rb_gsl_graph_T, 0);
+  rb_define_alias(cgsl_graph, "display_type", "T");
+
+  rb_define_method(cgsl_graph, "set_E", rb_gsl_graph_set_E, 1);
+  rb_define_alias(cgsl_graph, "E=", "set_E");
+  rb_define_alias(cgsl_graph, "axis_end=", "set_E");
+  rb_define_alias(cgsl_graph, "toggle_axis_end=", "set_E");
+  rb_define_method(cgsl_graph, "E", rb_gsl_graph_E, 0);
+  rb_define_alias(cgsl_graph, "axis_end", "E");
+  rb_define_alias(cgsl_graph, "toggle_axis_end", "E");
+
+  rb_define_method(cgsl_graph, "set_f", rb_gsl_graph_set_f, 1);
+  rb_define_alias(cgsl_graph, "f=", "set_f");
+  rb_define_alias(cgsl_graph, "font_size=", "set_f");
+  rb_define_method(cgsl_graph, "f", rb_gsl_graph_f, 0);
+  rb_define_alias(cgsl_graph, "font_size", "f");
+
+  rb_define_method(cgsl_graph, "set_F", rb_gsl_graph_set_F, 1);
+  rb_define_alias(cgsl_graph, "F=", "set_F");
+  rb_define_alias(cgsl_graph, "font_name=", "set_F");
+  rb_define_method(cgsl_graph, "F", rb_gsl_graph_F, 0);
+  rb_define_alias(cgsl_graph, "font_name", "F");
+
+  rb_define_method(cgsl_graph, "set_g", rb_gsl_graph_set_g, 1);
+  rb_define_alias(cgsl_graph, "g=", "set_g");
+  rb_define_alias(cgsl_graph, "grid_style=", "set_g");
+  rb_define_method(cgsl_graph, "g", rb_gsl_graph_g, 0);
+  rb_define_alias(cgsl_graph, "grid_style", "g");
+
+  rb_define_method(cgsl_graph, "set_h", rb_gsl_graph_set_h, 1);
+  rb_define_alias(cgsl_graph, "h=", "set_h");
+  rb_define_alias(cgsl_graph, "height=", "set_h");
+  rb_define_alias(cgsl_graph, "height_of_plot=", "set_h");
+  rb_define_method(cgsl_graph, "h", rb_gsl_graph_h, 0);
+  rb_define_alias(cgsl_graph, "height", "h");
+  rb_define_alias(cgsl_graph, "height_of_plot", "h");
+
+  rb_define_method(cgsl_graph, "set_k", rb_gsl_graph_set_k, 1);
+  rb_define_alias(cgsl_graph, "k=", "set_k");
+  rb_define_alias(cgsl_graph, "tick_size=", "set_k");
+  rb_define_method(cgsl_graph, "k", rb_gsl_graph_k, 0);
+  rb_define_alias(cgsl_graph, "tick_size", "k");
+
+  rb_define_method(cgsl_graph, "set_K", rb_gsl_graph_set_K, 1);
+  rb_define_alias(cgsl_graph, "K=", "set_K");
+  rb_define_alias(cgsl_graph, "clip_mode=", "set_K");
+  rb_define_method(cgsl_graph, "K", rb_gsl_graph_K, 0);
+  rb_define_alias(cgsl_graph, "clip_mode", "K");
+
+  rb_define_method(cgsl_graph, "set_l", rb_gsl_graph_set_l, 1);
+  rb_define_alias(cgsl_graph, "l=", "set_l");
+  rb_define_alias(cgsl_graph, "log_axis=", "set_l");
+  rb_define_alias(cgsl_graph, "toggle_log_axis=", "set_l");
+  rb_define_method(cgsl_graph, "l", rb_gsl_graph_l, 0);
+  rb_define_alias(cgsl_graph, "log_axis", "l");
+  rb_define_alias(cgsl_graph, "toggle_log_axis", "l");
+
+  rb_define_method(cgsl_graph, "set_L", rb_gsl_graph_set_L, 1);
+  rb_define_alias(cgsl_graph, "L=", "set_L");
+  rb_define_alias(cgsl_graph, "top_label=", "set_L");
+  rb_define_method(cgsl_graph, "L", rb_gsl_graph_L, 0);
+  rb_define_alias(cgsl_graph, "top_label", "L");
+
+  rb_define_method(cgsl_graph, "set_N", rb_gsl_graph_set_N, 1);
+  rb_define_alias(cgsl_graph, "N=", "set_N");
+  rb_define_alias(cgsl_graph, "no_tics=", "set_N");
+  rb_define_alias(cgsl_graph, "toggle_no_tics=", "set_N");
+  rb_define_method(cgsl_graph, "N", rb_gsl_graph_N, 0);
+  rb_define_alias(cgsl_graph, "no_tics", "N");
+  rb_define_alias(cgsl_graph, "toggle_no_tics", "N");
+
+  rb_define_method(cgsl_graph, "set_r", rb_gsl_graph_set_r, 1);
+  rb_define_alias(cgsl_graph, "r=", "set_r");
+  rb_define_alias(cgsl_graph, "right_shift=", "set_r");
+  rb_define_method(cgsl_graph, "r", rb_gsl_graph_r, 0);
+  rb_define_alias(cgsl_graph, "right_shift", "r");
+
+  rb_define_method(cgsl_graph, "set_R", rb_gsl_graph_set_R, 1);
+  rb_define_alias(cgsl_graph, "R=", "set_R");
+  rb_define_alias(cgsl_graph, "round-to-next-tick=", "set_R");
+  rb_define_alias(cgsl_graph, "toggle_round-to-next-tick=", "set_R");
+  rb_define_method(cgsl_graph, "R", rb_gsl_graph_R, 0);
+  rb_define_alias(cgsl_graph, "round-to-next-tick", "R");
+  rb_define_alias(cgsl_graph, "toggle_round-to-next-tick", "R");
+
+  rb_define_method(cgsl_graph, "set_s", rb_gsl_graph_set_s, 1);
+  rb_define_alias(cgsl_graph, "s=", "set_s");
+  rb_define_alias(cgsl_graph, "save_screen=", "set_s");
+  rb_define_method(cgsl_graph, "s", rb_gsl_graph_s, 0);
+  rb_define_alias(cgsl_graph, "save_screen", "s");
+
+  rb_define_method(cgsl_graph, "set_t", rb_gsl_graph_set_t, 1);
+  rb_define_alias(cgsl_graph, "t=", "set_t");
+  rb_define_alias(cgsl_graph, "transpose_axes=", "set_t");
+  rb_define_alias(cgsl_graph, "toggle_transpose_axes=", "set_t");
+  rb_define_method(cgsl_graph, "t", rb_gsl_graph_t, 0);
+  rb_define_alias(cgsl_graph, "transpose_axes", "t");
+  rb_define_alias(cgsl_graph, "toggle_transpose_axes", "t");
+
+  rb_define_method(cgsl_graph, "set_u", rb_gsl_graph_set_u, 1);
+  rb_define_alias(cgsl_graph, "u=", "set_u");
+  rb_define_alias(cgsl_graph, "upward_shift=", "set_u");
+  rb_define_method(cgsl_graph, "u", rb_gsl_graph_u, 0);
+  rb_define_alias(cgsl_graph, "upward_shift", "u");
+
+  rb_define_method(cgsl_graph, "set_w", rb_gsl_graph_set_w, 1);
+  rb_define_alias(cgsl_graph, "w=", "set_w");
+  rb_define_alias(cgsl_graph, "width=", "set_w");
+  rb_define_alias(cgsl_graph, "width_of_plot=", "set_w");
+  rb_define_method(cgsl_graph, "w", rb_gsl_graph_w, 0);
+  rb_define_alias(cgsl_graph, "width", "w");
+  rb_define_alias(cgsl_graph, "width_of_plot", "w");
+
+  rb_define_method(cgsl_graph, "set_x", rb_gsl_graph_set_x, 1);
+  rb_define_alias(cgsl_graph, "x=", "set_x");
+  rb_define_alias(cgsl_graph, "x_limits=", "set_x");
+  rb_define_method(cgsl_graph, "x", rb_gsl_graph_x, 0);
+  rb_define_alias(cgsl_graph, "x_limits", "x");
+
+  rb_define_method(cgsl_graph, "set_y", rb_gsl_graph_set_y, 1);
+  rb_define_alias(cgsl_graph, "y=", "set_y");
+  rb_define_alias(cgsl_graph, "y_limits=", "set_y");
+  rb_define_method(cgsl_graph, "y", rb_gsl_graph_y, 0);
+  rb_define_alias(cgsl_graph, "y_limits", "y");
+
+  rb_define_method(cgsl_graph, "set_X", rb_gsl_graph_set_X, 1);
+  rb_define_alias(cgsl_graph, "X=", "set_X");
+  rb_define_alias(cgsl_graph, "x_label=", "set_X");
+  rb_define_method(cgsl_graph, "X", rb_gsl_graph_X, 0);
+  rb_define_alias(cgsl_graph, "x_label", "X");
+  rb_define_method(cgsl_graph, "set_Y", rb_gsl_graph_set_Y, 1);
+  rb_define_alias(cgsl_graph, "Y=", "set_Y");
+  rb_define_alias(cgsl_graph, "y_label=", "set_Y");
+  rb_define_method(cgsl_graph, "Y", rb_gsl_graph_Y, 0);
+  rb_define_alias(cgsl_graph, "y_label", "Y");
+
+  rb_define_method(cgsl_graph, "set_bg", rb_gsl_graph_set_bg, 1);
+  rb_define_alias(cgsl_graph, "bg=", "set_bg");
+  rb_define_alias(cgsl_graph, "bg_color=", "set_bg");
+  rb_define_alias(cgsl_graph, "set_bg_color", "set_bg");
+  rb_define_method(cgsl_graph, "bg", rb_gsl_graph_bg, 0);
+  rb_define_alias(cgsl_graph, "bg_color", "bg");
+
+  rb_define_method(cgsl_graph, "set_bitmap_size", rb_gsl_graph_set_bitmap_size, 1);
+  rb_define_alias(cgsl_graph, "bitmap_size=", "set_bitmap_size");
+  rb_define_alias(cgsl_graph, "set_bitmap_size", "set_bitmap_size");
+  rb_define_method(cgsl_graph, "bitmap_size", rb_gsl_graph_bitmap_size, 0);
+
+  rb_define_method(cgsl_graph, "set_frame", rb_gsl_graph_set_frame, 1);
+  rb_define_alias(cgsl_graph, "frame=", "set_frame");
+  rb_define_alias(cgsl_graph, "frame_color=", "set_frame");
+  rb_define_alias(cgsl_graph, "set_frame_color", "set_frame");
+  rb_define_method(cgsl_graph, "frame", rb_gsl_graph_frame, 0);
+  rb_define_alias(cgsl_graph, "frame_color", "frame");
+
+  rb_define_method(cgsl_graph, "set_frame_line_width", rb_gsl_graph_set_frame_line_width, 1);
+  rb_define_alias(cgsl_graph, "frame_line_width=", "set_frame_line_width");
+  rb_define_method(cgsl_graph, "frame_line_width", rb_gsl_graph_frame_line_width, 0);
+
+  rb_define_method(cgsl_graph, "set_max_line_length", rb_gsl_graph_set_max_line_length, 1);
+  rb_define_alias(cgsl_graph, "max_line_length=", "set_max_line_length");
+  rb_define_method(cgsl_graph, "max_line_length", rb_gsl_graph_max_line_length, 0);
+
+  rb_define_method(cgsl_graph, "set_page_size", rb_gsl_graph_set_page_size, 1);
+  rb_define_alias(cgsl_graph, "page_size=", "set_page_size");
+  rb_define_method(cgsl_graph, "page_size", rb_gsl_graph_page_size, 0);
+
+  rb_define_method(cgsl_graph, "set_pen_colors", rb_gsl_graph_set_pen_colors, 1);
+  rb_define_alias(cgsl_graph, "pen_colors=", "set_pen_colors");
+  rb_define_alias(cgsl_graph, "pen=", "set_pen_colors");
+  rb_define_method(cgsl_graph, "pen_colors", rb_gsl_graph_pen_colors, 0);
+  rb_define_alias(cgsl_graph, "pen", "pen_colors");
+
+  rb_define_method(cgsl_graph, "set_rotation", rb_gsl_graph_set_rotation, 1);
+  rb_define_alias(cgsl_graph, "rotation=", "set_rotation");
+  rb_define_method(cgsl_graph, "rotation", rb_gsl_graph_rotation, 0);
+
+  rb_define_method(cgsl_graph, "set_title_font_name", rb_gsl_graph_set_title_font_name, 1);
+  rb_define_alias(cgsl_graph, "title_font_name=", "set_title_font_name");
+  rb_define_method(cgsl_graph, "title_font_name", rb_gsl_graph_title_font_name, 0);
+
+  rb_define_method(cgsl_graph, "set_title_font_size", rb_gsl_graph_set_title_font_size, 1);
+  rb_define_alias(cgsl_graph, "title_font_size=", "set_title_font_size");
+  rb_define_method(cgsl_graph, "title_font_size", rb_gsl_graph_title_font_size, 0);
+
+  rb_define_method(cgsl_graph, "set_rotate_y_label", rb_gsl_graph_set_rotate_y_label, 1);
+  rb_define_alias(cgsl_graph, "rotate_y_label=", "set_rotate_y_label");
+  rb_define_alias(cgsl_graph, "toggle_rotate_y_label=", "set_rotate_y_label");
+  rb_define_method(cgsl_graph, "rotate_y_label", rb_gsl_graph_rotate_y_label, 0);
+  rb_define_alias(cgsl_graph, "toggle_rotate_y_label", "rotate_y_label");
+
+  rb_define_method(cgsl_graph, "set_B", rb_gsl_graph_set_B, 1);
+  rb_define_alias(cgsl_graph, "B=", "set_B");
+  rb_define_alias(cgsl_graph, "auto_dump=", "set_B");
+  rb_define_alias(cgsl_graph, "toggle_auto_dump=", "set_B");
+  rb_define_method(cgsl_graph, "B", rb_gsl_graph_B, 0);
+  rb_define_alias(cgsl_graph, "auto_dump", "B");
+  rb_define_alias(cgsl_graph, "toggle_auto_dump", "B");
+
+  rb_define_method(cgsl_graph, "set_I", rb_gsl_graph_set_I, 1);
+  rb_define_alias(cgsl_graph, "I=", "set_I");
+  rb_define_alias(cgsl_graph, "input_format=", "set_I");
+  rb_define_method(cgsl_graph, "I", rb_gsl_graph_I, 0);
+  rb_define_alias(cgsl_graph, "input_format", "I");
+
+  rb_define_method(cgsl_graph, "set_m", rb_gsl_graph_set_m, 1);
+  rb_define_alias(cgsl_graph, "m=", "set_m");
+  rb_define_alias(cgsl_graph, "line_mode=", "set_m");
+  rb_define_method(cgsl_graph, "m", rb_gsl_graph_m, 0);
+  rb_define_alias(cgsl_graph, "line_mode", "m");
+
+  rb_define_method(cgsl_graph, "set_S", rb_gsl_graph_set_S, 1);
+  rb_define_alias(cgsl_graph, "S=", "set_S");
+  rb_define_alias(cgsl_graph, "symbol=", "set_S");
+  rb_define_method(cgsl_graph, "S", rb_gsl_graph_S, 0);
+  rb_define_alias(cgsl_graph, "symbol", "S");
+
+  rb_define_method(cgsl_graph, "set_W", rb_gsl_graph_set_W, 1);
+  rb_define_alias(cgsl_graph, "W=", "set_W");
+  rb_define_alias(cgsl_graph, "line_width=", "set_W");
+  rb_define_method(cgsl_graph, "W", rb_gsl_graph_W, 0);
+  rb_define_alias(cgsl_graph, "line_width", "W");
+
+  rb_define_method(cgsl_graph, "set_q", rb_gsl_graph_set_q, 1);
+  rb_define_alias(cgsl_graph, "q=", "set_q");
+  rb_define_alias(cgsl_graph, "fill_fraction=", "set_q");
+  rb_define_method(cgsl_graph, "q", rb_gsl_graph_q, 0);
+  rb_define_alias(cgsl_graph, "fill_fraction", "q");
+
+  rb_define_method(cgsl_graph, "set_C", rb_gsl_graph_set_C, 1);
+  rb_define_alias(cgsl_graph, "C=", "set_C");
+  rb_define_alias(cgsl_graph, "use_color=", "set_C");
+  rb_define_alias(cgsl_graph, "toggle_use_color=", "set_C");
+  rb_define_method(cgsl_graph, "C", rb_gsl_graph_C, 0);
+  rb_define_alias(cgsl_graph, "use_color", "C");
+  rb_define_alias(cgsl_graph, "toggle_use_color", "C");
+
+  rb_define_method(cgsl_graph, "set_symbol_font_name", 
+		   rb_gsl_graph_set_symbol_font_name, 1);
+  rb_define_alias(cgsl_graph, "symbol_font_name=", "set_symbol_font_name");
+  rb_define_method(cgsl_graph, "symbol_font_name", rb_gsl_graph_symbol_font_name, 0);
+
+  rb_define_method(cgsl_graph, "set_reposition", 
+		   rb_gsl_graph_set_reposition, 1);
+  rb_define_alias(cgsl_graph, "reposition=", "set_reposition");
+  rb_define_method(cgsl_graph, "reposition", rb_gsl_graph_reposition, 0);
+
+  rb_define_method(cgsl_graph, "set_blankout", 
+		   rb_gsl_graph_set_blankout, 1);
+  rb_define_alias(cgsl_graph, "blankout=", "set_blankout");
+  rb_define_method(cgsl_graph, "blankout", rb_gsl_graph_blankout, 0);
+
+  rb_define_method(cgsl_graph, "set_O", rb_gsl_graph_set_O, 1);
+  rb_define_alias(cgsl_graph, "O=", "set_O");
+  rb_define_alias(cgsl_graph, "portable_output=", "set_O");
+  rb_define_method(cgsl_graph, "O", rb_gsl_graph_O, 0);
+  rb_define_alias(cgsl_graph, "portable_output", "O");
+}