float sf_quantile(int q /* quantile */, int n /* array length */, float* a /* array [n] */) /*< find quantile (caution: a is changed) >*/ { ... }