How to analyze a acid and salt mixture using titrimetry?
Prologue

The following will show you not only how to sketch a titration curve but how to produce an analytical form of a titration curve. So this might not be an easy solution even if I try to keep it as easy as possible. There is a wonderful yet german book that has (to my knowledge) never been translated to english which describes everything of the following in detail:


Bliefert, Claus; Linek, Alfred; Morawietz, Gerd (1978): PH-Wert-Berechnungen. Weinheim: Verlag Chemie.


There might also be some imprecise terms because I don't know many of the english counterparts to the german terms I am familiar with.

Answer

Let's think a bit about your problem first to simplify it from scratch:


As you mentioned, reaction (a) will take place right after all protons from hydrochloric acid have reacted and before reaction (b) and (c) would start.
The pK$_a$ of (b) and (c) are so close to another that they are unable to be separeted by titration and
the pK$_a$ of (b) and (c) are so high that they are unable to be determined by titration in diluted aqueous solutions.


If we take a look at the titration curve of phosphoric acid, we can see what is meant by "too high to be determined". You can clearly identify the two inflection points for the first two protons but there is no third inflection point that could be determined. What is likewise to be seen, is that one can neglect this third deprotonation step in titration curves of diluted solutions (yellow vs. blue titration curve).



Right now we can reduce your problem to a mixture of hydrochloric acid and M$^{3+}$. To continue we need to establish a proton condition (Protonen-Bedingung):

$$c(H_3O^+) = c(OH^-) + c(A^-) + c(X^-) - c(HB^+) \tag{1}$$

It means that protons come from the autoprotolysis of water (c(OH$^-$)), from hydrochloric acid (c(A$^-$)), from your metal ion (c(X$^-$)) and that the protons react with the used base (c(HB$^+$)).

To define an equation for the acid residue anions a small equation system (where small c's denote actual concentrations in solution and the big C denotes the origin concentration)
$$\begin{align}c(HA) \cdot k_a - c(H_3O^+) \cdot c(A^-) =&amp;~0\\c(HA)+c(A^-) =&amp;~C(HA)\end{align} \tag{2.1}$$
needs to be solved to give
$$c(A^-) = \frac{C(HA) \cdot k_a}{c(H_3O^+)+k_a}. \tag{2.2}$$

This can be done analogously for the other acid and also for the base.

What remains is the concentration of the hydroxide ions from the autoprotolysis of water that is simply determined by rearranging the equation for the autoprotolysis itself:

$$c(H_3O^+) \cdot c(OH^-) = k_W \Leftrightarrow c(OH^-) = \frac{k_W}{c(H_3O^+)} \tag{3}$$

Combining everything so far we get the following proton condition:
$$x = \frac{k_W}{x} + \frac{C(HA) \cdot k_{a,HA}}{x+k_{a,HA}} + \frac{C(HX) \cdot k_{a,HX}}{x+k_{a,HX}} - \frac{x \cdot C(B) \cdot k_b}{x \cdot k_b + k_W} \tag{4}$$

where $c(H_3O^+)$ has been replaced by $x$ for the sake of simplicity.

The plot of titration curves usually shows the pH over the base volume (see eq. 8) or the degree of neutralisation ($\tau$ or $\gamma$) that is defined as $$\tau = \frac{c(B)}{C(HA)} \tag{5.1}$$

In the case of your mixture of two acids - and as always neglecting the increase in volume - it can be written as
$$\tau = \frac{c(B)}{C(HA)+C(HX)} \tag{5.2}$$

To sketch the titration curve let's say that we know the two concentrations of the acids. To calculate (5.2) we also need to know c(B) that can easily be obtained by rearranging (4) to get:
$$c(B) = \left(1+\frac{k_w}{x \cdot k_b}\right)\left(\frac{k_w}{x}-x+\frac{C(HA) \cdot k_{a,HA}}{x+k_{a,HA}} + \frac{C(HX) \cdot k_{a,HX}}{x+k_{a,HX}}\right) \tag{6}$$

If one decides to include the increasing volume the concentrations of the acids and the base have to be extended to give
$$\begin{align}C(HA) \rightarrow&amp;~\frac{C(HA) \cdot V(HA)}{V(HA)+V(B)}\\C(B) \rightarrow&amp;~\frac{C(B) \cdot V(B)}{V(HA)+V(B)}\end{align} \tag{7}$$

Inserting (7) in (4) and rearranging for V(B) we finally get a neat equation

$$V(B)=\frac{V(HA)~\left(\frac{k_w}{x}-x+\frac{C(HA) \cdot k_{a,HA}}{x+k_{a,HA}} + \frac{C(HX) \cdot k_{a,HX}}{x+k_{a,HX}}\right)}{\left(x-\frac{k_w}{x}+\frac{C(B) \cdot k_b}{\frac{k_w}{x}+k_b}\right)} \tag{8}$$

As the pH is defined as $pH = -log_{10}(c(H_3O^+))$ and $c(H_3O^+)=x$ you can replace every $x$ in equation 8 with $10^{-x}$.

Now you are able to calculate with e.g. Excel enough (V(B),pH)-pairs to draw a titration curve and get some information about your system of two acids during the titration with your desired base.

The following figure shows some examples of your mixture using equation 8.



Please let me know if someone is interested in the full equation including also the second and third deprotonation step with pK$_a$ values 13.3 and 13.5 and I will add it here.