function [z] = zquant(quantile)
z = 4.91 * (power(quantile,.14) - power(1-quantile,.14));