c_rmultinom

#define MATHLIB_STANDALONE 1
#include "Rmath.h"
#include <stdio.h>


int main(int argc, char** argv) {
    int draws = 100;
    int classes = 12;
    int vals[classes];
    double probs[5] = {0.1,0.2,1.3,0.2,0.3};

    set_seed(123, 456);
    rmultinom(draws, probs, classes, vals);

    for(int j=0; j < classes; j++) {
        printf("Count of class %i drawn: %i\n", j, vals[j]);
    }
    return 0;
}



// Outputs :
// (gdb) p vals
// $3 = {0, 0, -1, 32758, 0, 0, 0, 0, 6290992, 0, 1635326112, 32758}