% function y=compeq([a,c],f) % calculates our favorite comparametric equation function y=compeq([a,c],f) num=f^(1/c)*k^a; den=f^(1/c)*(k^a-1)+1; y = num/den; y=y^c; end