function y = prod_cd(k, alpha)
    %PROD_CD Cobb-Douglas production y = k^alpha
    %   y = PROD_CD(k, alpha) computes k.^alpha elementwise
        y = k.^alpha;
    end