function[q] = profile(V,n,C, V_max)
    if (0 <= V && V <= n)
        q = C;
        return q;
    elseif (n < V && V <= V_max)
        q = (V_max - V)*C/(V_max - n);
        return q;
    else
        q = 0;
        return q;
    end
end