I want to learn how solve simple ax+by=c with matrices (assuming that's the fasted method?),
but it's difficult to find correct learning material.
I've been through this process:
4386x + 89744y = 2
NB: Matrices are unknown terrotorium, go easy, ty.
(links are also appreciated)

You don't need matrices really, only use the Extended Euclidean algorithm. Here is a layout for it. I continued one line further, in order to obtain the quotient of the two given numbers by their gcd:
Thus you have a solution: $\;(x,y)=(-22\,303,1090)$. The theory tells us all other solutions have the form: $$(x,y)=\Bigl(-22\,303+k44\,872,1090-k2193),\quad k\in\mathbf Z.$$
Some explanations:
Not only the gcd is a linear combination of the two given numbers, but all intermediate remainders (easy induction). Si if $r_i+1$ denotes the remainder and $q_i$ the quotient in the $i$-th division (initialisation for the remainders: $r_0=89\,744$, $r_1=4386$), we have $$r_{i-1}=q_ir_i+r_{i+1}\iff r_{i+1}=r_{i-1}-q_ir_i,$$ and of course the same relation is valid for the coefficients $u_i,v_i$, by linearity.