#include <iostream>
#include <complex>
#include <blitz/array.h>

int main(int argc, char** argv) 
{
    const int N = 100;
    blitz::Array<double,1> A(N);
    blitz::Range even(2,N-2,2), odd(1,N-3,2);
    
    A(even) = 0.5 * (A(even+1) - A(even-1));
    A(odd ) = 0.5 * (A(odd+1)  - A(odd-1));
    
    blitz::Array<double,1> Aeven = A( even );
    blitz::Array<double,1> Aodd  = A( odd  );
    
    Aeven = blitz::central12( Aeven );
    Aodd  = blitz::central12( Aodd  );     
    return 0;
}