import mir.test: shouldApprox;

cauchyInvCDF(0.0).shouldApprox == -double.infinity;
cauchyInvCDF(0.1).shouldApprox == -3.077684;
cauchyInvCDF(0.2).shouldApprox == -1.376382;
cauchyInvCDF(0.3).shouldApprox == -0.7265425;
cauchyInvCDF(0.4).shouldApprox == -0.3249197;
cauchyInvCDF(0.5).shouldApprox == 0.0;
cauchyInvCDF(0.6).shouldApprox == 0.3249197;
cauchyInvCDF(0.7).shouldApprox == 0.7265425;
cauchyInvCDF(0.8).shouldApprox == 1.376382;
cauchyInvCDF(0.9).shouldApprox == 3.077684;
cauchyInvCDF(1.0).shouldApprox == double.infinity;

// Can include location/scale
cauchyInvCDF(0.2, 1, 2).shouldApprox == -1.752764;
cauchyInvCDF(0.4, 1, 2).shouldApprox == 0.3501606;
cauchyInvCDF(0.6, 1, 2).shouldApprox == 1.649839;
cauchyInvCDF(0.8, 1, 2).shouldApprox == 3.752764;