** This is part of the CTSim program
** Copyright (c) 1983-2001 Kevin Rosenberg
**
-** $Id: cubicinterp.cpp,v 1.3 2001/02/09 14:34:16 kevin Exp $
+** $Id: cubicinterp.cpp,v 1.4 2001/02/09 21:27:51 kevin Exp $
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
double* temp = new double [n - 1];
temp[0] = 0;
- for (int i = 1; i < n - 1; i++) {
+ int i;
+ for (i = 1; i < n - 1; i++) {
double t = 2 + (0.5 * m_pdY2[i-1]);
temp[i] = y[i+1] + y[i-1] - y[i] - y[i];
temp[i] = (3 * temp[i] - 0.5 * temp[i-1]) / t;