sistemi lineari con fuction matlab

[mpulcina]

mi potreste dare una mano con questo esercizio? non riesco neanche ad impostarlo.
la soluzione del sistema lineare a matrice triangolare sia superiore che inferiore si ottiene costruendo le seguenti sequenze ricorsive
$\{($f_0$=0),($f_k$=$(-a_k,k+1)/$(a_(k,k-1)*f_k-1+a_k,k)):}$ con k=1,...,n-1
$\{($g_0$=0),($g_k=$(b_k-a_k,k-1*g_k-1)/$(a_k,k-1*f_k-1+a_k,k)):}$ con k=1,...,n-1
e calcolando quindi x_n=g_n, x_i0f_j*xx_j+1, con j=n-1,n-2,...1.
scrivere una fuction matlab per implementare l'algoritmo descritto(riportare l'intera sequenza di istruzioni)

P.S. ho visto che le formule non sono uscite bene nonostante abbia usato il codice ascii, comunque sono 2 sistemi, il primo ha 2eq che sono: f di 0=0 e f di k=(-a di k,k+1)/(a di (k,k-1)*f di (k-1) + a di k,k) con k=0,..,n-1 e il secondo g di 0=0 e g di k= (b di k - a di (k,k-1)* g di (k-1))/(a di (k,k-1)*f di (k-1) + a di k,k) con k=0,..,n-1 e calcolando quindi x di n=g di n, x di j= f di j*x di (j+1)+g di j, con j=n-1,n-2,...1.
scrivere una fuction matlab per implementare l'algoritmo descritto(riportare l'intera sequenza di istruzioni)