Allora, il problema si lascia risolvere abbastanza bene in coordinate baricentriche rispetto al triangolo ortico $H_aH_bH_c$ di $ABC$; in quanto:
1) Ovviamente per la scelta del triangolo di riferimento valgono $H_a=[1:0:0]$ e cicliche.
2) Grazie ad una ben nota proprietà che si dimostra per angle chasing, $A$ è l'excentro opposto al vertice $H_a$ del triangolo ortico (e cicliche), e valgono dunque $A=[-a

c]$ e cicliche.
3) Ovviamente la circonferenza di Feuerbach di $ABC$ è (per definizione) la circoscritta ad $H_aH_bH_c$, ed ha dunque la semplice equazione $a^2yz+b^2xz+c^2xy=0$.
4)I punti medi $M_a,M_b,M_c$ si trovano sommando le coordinate normalizzate dei due excentri "giusti", ovvero
$$M_a=\left[\frac{a}{a-b+c}+\frac{a}{a+b-c}:\frac{-b}{a-b+c}+\frac{b}{a+b-c}:\frac{c}{a-b+c}+\frac{-c}{a+b-c} \right]$$
Che può essere scritto meglio moltiplicando tutto per $(a-b+c)(a+b-c)$ eliminando i denominatori:
$$M_a=[a(a+b-c)+a(a-b+c): -b(a+b-c)+b(a-b+c):c(a+b-c)-c(a-b+c)]\Rightarrow$$
$$\Rightarrow\ \ M_a=[a^2:bc-b^2:bc-c^2]\ \ \ \ ; \ \ \ \ M_b=[ac-a^2:b^2:ac-c^2]\ \ \ \ ; \ \ \ \ M_c=[ab-a^2:ab-b^2:c^2]$$
Adesso si tratta di scrivere le rette $l_a, l'_a$ e cicliche; che in effetti sono le
polari di punti abbastanza belli rispetto alla circoscritta
grazie Drago!, e quindi possiamo usare un lemma che non conoscevo ma che semplifica molte cose
Lemma overkill sulle coniche: Dato un punto $(u

w)$ e una qualsiasi circumconica di equazione $pyz+qxz+rxy=0$, la sua polare rispetto alla conica è la retta di equazione
$$p(wy+vz)+q(uz+wx)+r(vx+uy)=0$$
L'enunciato di questo lemma l'ho trovato su
"Introduction to the Geometry of the Triangle, Paul Yiu", quindi boh, spero sia vero... di sicuro con Strong EFFT (che però non so usare troppo bene; dunque non mi fidavo troppo) mi veniva la stessa cosa, quindi almeno un po' ci spero che funzioni

.
Dal lemma è immediato ricavare che
$$l_a: \ a^2(0\cdot y+0\cdot z)+b^2(1\cdot z+0\cdot x)+c^2(0\cdot x+1\cdot y)=0\ \ \Rightarrow \ l_a:\ b^2z+c^2y=0$$
$$l'_a:\ a^2[(bc-c^2)\cdot y+(bc-b^2)\cdot z]+b^2[a^2\cdot z+(bc-c^2)\cdot x]+c^2[(bc-b^2)\cdot x+a^2\cdot y]=0\ \ \Rightarrow \ l'_a:\ (b-c)^2x+a^2y+a^2z=0$$
E quindi $A'$ si trova abbastanza facilmente come intersezione delle due rette (si semplifica pure un fattore $b-c$...)
$$A'=[a^2b+a^2c:-b^3+b^2c:bc^2-c^3]$$
Adesso per la retta $AA'$ operiamo allo stesso modo e troviamo
$$AA':\ \ [-bc(b+c)(b-c): -ac(bc-c^2+ab+ac): ab(ab+ac-b^2+bc)]$$
$$BB':\ \ [bc(bc+ab-c^2+ac): -ac(c+a)(c-a):-ab(ac-a^2+bc+ab)]$$
$$CC'=\ \ [-bc(ab-b^2+ac+bc): ac(ac+bc-a^2+ab): -ab(a+b)(a-b)]$$
Se adesso chiamiamo $P=ab+bc+ca$ e troviamo l'intersezione delle prime due, otteniamo (con un ultimo sforzo)
$$J=[aP^2-2a^3P+a^3c^2+a^3b^2-ab^2c^2: bP^2+b^3c^2+b^3a^2-2b^3P-a^2bc^2: cP^2+b^2c^3+a^2c^3-2Pc-a^2b^2c]$$
Che ha le coordinate cicliche ($P$ non cambia ciclando), e quindi le tre rette concorrono