new file mode 100644
--- /dev/null
+++ b/liboctave/Quad.h
@@ -0,0 +1,114 @@
+// Quad.h                                           -*- C++ -*-
+/*
+
+Copyright (C) 1992, 1993 John W. Eaton
+
+This file is part of Octave.
+
+Octave is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 2, or (at your option) any
+later version.
+
+Octave is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with Octave; see the file COPYING.  If not, write to the Free
+Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
+
+*/
+
+#if !defined (_Quad_h)
+#define _Quad_h 1
+
+#ifdef __GNUG__
+#pragma interface
+#endif
+
+#include "Matrix.h"
+
+#ifndef Vector
+#define Vector ColumnVector
+#endif
+
+#ifndef _Quad_typedefs
+#define _Quad_typedefs 1
+
+typedef double (*integrand_fcn) (double x);
+
+#endif
+
+class Quad
+{
+ public:
+
+  Quad (integrand_fcn fcn);
+  Quad (integrand_fcn fcn, double abs, double rel);
+
+  virtual double integrate (void);
+  virtual double integrate (int& ier);
+  virtual double integrate (int& ier, int& neval);
+  virtual double integrate (int& ier, int& neval, double& abserr) = 0;
+
+ protected:
+
+  double absolute_tolerance;
+  double relative_tolerance;
+
+  integrand_fcn f;
+};
+
+class DefQuad : public Quad
+{
+ public:
+
+  DefQuad (integrand_fcn fcn);
+  DefQuad (integrand_fcn fcn, double ll, double ul);
+  DefQuad (integrand_fcn fcn, double ll, double ul, double abs, double rel);
+  DefQuad (integrand_fcn fcn, double ll, double ul, const Vector& sing);
+  DefQuad (integrand_fcn fcn, const Vector& sing, double abs, double rel);
+  DefQuad (integrand_fcn fcn, const Vector& sing);
+  DefQuad (integrand_fcn fcn, double ll, double ul, const Vector& sing,
+	   double abs, double rel);
+
+  double integrate (int& ier, int& neval, double& abserr);
+
+ private:
+
+  double lower_limit;
+  double upper_limit;
+
+  Vector singularities;
+};
+
+class IndefQuad : public Quad
+{
+ public:
+
+  enum IntegralType { bound_to_inf, neg_inf_to_bound, doubly_infinite };
+
+  IndefQuad (integrand_fcn fcn);
+  IndefQuad (integrand_fcn fcn, double b, IntegralType t);
+  IndefQuad (integrand_fcn fcn, double b, IntegralType t, double abs,
+	     double rel);
+  IndefQuad (integrand_fcn fcn, double abs, double rel);
+
+  double integrate (int& ier, int& neval, double& abserr);
+
+ private:
+
+  double bound;
+  IntegralType type;
+};
+
+#endif
+
+/*
+;;; Local Variables: ***
+;;; mode: C++ ***
+;;; page-delimiter: "^/\\*" ***
+;;; End: ***
+*/