
H=[0:0.25:14];

Fe=3.2-3*H;
FeOH=1.0-2*H;
FeOH2=-2.5-H;
FeOH4=-18.4+H;
Fe2OH2=3.5-4*H;
Fe3OH4=3.3-5*H;
soluble=10.^Fe + 10.^FeOH +10.^FeOH2 +10.^FeOH4 +10.^Fe2OH2 +10.^Fe3OH4;
solid=log10(soluble);
plot(H,Fe,'k--',H,FeOH,'b--',H,FeOH2,'k--',H,FeOH4,'k--',H,Fe2OH2,'b--',H,Fe3OH4,'k--',H,solid,'r')
axis([0 14 -14 0]);
Title('Solubility diagram for Iron (III) Hydroxide in Pure Water')
xlabel('pH');
ylabel('log [C]');
gtext('Fe(OH)3(s)')
gtext('FeOH 2+')
gtext('Fe(OH)2 +')
gtext('Fe(OH)4 -')
gtext('Fe2(OH)2 +4')
gtext('Fe3(OH)4 +5')
gtext('Fe +3')
