// utility functions
#ifndef UTILH
#define UTILH

#include <cstdlib>
#include <iostream>

void error( char s[] );

int min( int i, int j);

inline double sq( double x )
{ return x*x; }
#endif