#include <stdio.h>
#include <cstdlib>
int main(int argc, char **argv)
{
    try {
        int device = argc > 1 ? atoi(argv[1]) : 0;
        int days = 9, sites = 4;
        int n = 10; 
        float day_[] =         {0, 0, 1, 2, 5, 5, 6, 6, 7, 8 }; 
        float site_[] =        {2, 3, 0, 1, 1, 2, 0, 1, 2, 1 };
        float measurement_[] = {9, 5, 6, 3, 3, 8, 2, 6, 5, 10}; 
        array measurement(n,measurement_);
 
            rainfall(s) = 
sum(measurement * (site == s));
        }
        printf("total rainfall at each site:\n");
        array is_between = 1 <= day && day <= 5; 
 
        float rain_between = sum<float>(measurement * is_between);
        printf("rain between days: %g\n", rain_between);
        printf(
"number of days with rain: %g\n", sum<float>(
diff1(day) > 0) + 1);
            per_day(d) = 
sum(measurement * (day == d));
        printf("total rainfall each day:\n");
        printf("number of days over five: %g\n", sum<float>(per_day > 5));
        fprintf(stderr, 
"%s\n", e.
what());
        throw;
    }
    #ifdef WIN32 // pause in Windows
    if (!(argc == 2 && argv[1][0] == '-')) {
        std::cout << "hit [enter]...";
        fflush(stdout);
        getchar();
    }
    #endif
    return 0;
}