Здравствуйте, Botsman.
Объём раствора V=V_NaOH+V_HCl
При этом разница количеств кислоты и основания (в данном случае также разность количеств ионов [H⁺] и [OH^-])
n(H⁺)-n(OH^-)=V_HCl*C_HCl-V_NaOH*C_NaOH
разница концентраций [H⁺] и [OH^-]
[$916$]C=[H⁺]-[OH^-]=(n(H⁺)-n(OH^-))/V=(V_HCl*C_HCl-V_NaOH*C_NaOH)/(V_NaOH+V_HCl)
до достижения точки эквивалентности это значение практически равно концентрации [H⁺] и pH=-lg([$916$]C), после точки эквивалентности это значение практически равно концентрации [OH^-] и pH=14+lg(-[$916$]C)
Заметные отклонения этих функций от pH наблюдается только вблизи точки эквивалентности (при отклонении от данной точки не превышающем 0,001%)

Либо выведем точную зависимость [H⁺] от [$916$]C
[$916$]C=[H⁺]-[OH^-]=[H⁺]-10^-14/[H⁺]
[H⁺]²-[$916$]C*[H⁺]-10^-14=0
[H⁺]=([$916$]C+[$8730$]([$916$]C²+4*10^-14))/2
pH=-lg(([$916$]C+[$8730$]([$916$]C²+4*10^-14))/2)

В данном случае (сильная кислота и сильное основание) рН вблизи точки эквивалентности очень резко изменяется от 4-5 до 9-10. Подойдёт любой индикатор, меняющий окраску в диапазоне 4-9 (метиловый красный, лакмоид, нейтральный красный, п-нитрофенол, м-нитрофенол, феноловый красный и т. д.). Даже метилоранж и фенолфталеин будут давать результат с минимальным отклонением от точки эквивалентности.

В приложении файл с рассчётом кривой титрования.