// Creates vect of size end+1, by Tim Jones @ Drexel University

double *vect(int end)
{
  double *temp;
  temp=(double *)malloc((size_t) ((end+1)*sizeof(double)));
  assert(temp != NULL);
  return temp;
}