Weibull distribution

weibull.c

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

float randomf() {
	return (float)random() / (float)RAND_MAX;
}

float weibull(float alpha, float beta) {
	float u = 1.0 - randomf();
	return alpha * pow(-log(u), 1.0/beta);
}