Matlab - Ayuda urgente con algoritmo Bit Filling

Hola a todos, por favor necesito de su ayuda para programar el algoritmo de Bit Filling, diseñado por J. Campello. Este se puede ver en http://www.almaden.ibm.com/cs/people/dmodha/designing-ldpc-published.pdf
pero consiste básicamente en construir una matriz LDPC teniendo en cuenta de que no se debe infringir las restricciones de girth y de grado de chequeo...yo no se mucho de programación, asi que cualquier ayuda se agradece :)...

El programa es el siguiente:

Entradas: m, n, aj, br, gmáx, gmin.
1_ Set n=0, A=|m|, y U1=0, g=gmáx
2_ Para c ∈ |m|, set Grado(c)=0 y Nc=0
3_ Haga {
4_ ∀ c ∈ U1, setee Hc,n = 1 e incremente Grado(c) en uno
5_ Setee i = 0, U1 = 0, U = 0
6_ Haga {
7_ Calcule Fo =A\U
8_ Si (Fo 0) {
9_ Escoja c* desde Fo de acuerdo a la heurística {

9A_ Setee j = 0, Eo = Fo
9B_ ∀ c ∈ Fo, sea U1(c) = {c}
9C_ Mientras (|Fj|>1) y (Ej = Fj) {
9D_ j =j+1
9E_ ∀ c ∈ Fj-1, Calcule Wj(c) =
9F_ Calcule Fj = {c1 ∈ Fj-1: ∀ c2 ∈ Fj-1 Wj(c1) Wj(c2)}
9G_ ∀ c ∈ Fj, Calcule Vj+1(c) Usando V2(c)= Uj(c)
9H_ Calcule Ej = { c ∈ Fj : Vj+1(c) Vj (c)}
9I_ } Fin Mientras
9J_ Si (|Fj|=1)
9K_ Seleccione c* desde Fj
9L_ De otra manera
9M_ Seleccione c* desde Fj\Ej
9N_ Fin Si
}

10_ ∀ c ∈ U1, descargue Nc = Nc {c*} y descargue Nc* = Nc* U1
11_ Descargue U1 = U1 {c*}, U = U
12_ i = i+1
13_ } De otra manera {
14_ g = g – 2
15_ Recalcule U =
16_ }
17_ } Mientras (i<
18_ n = n+1
19_ } Mientras [(n<N) y (g
Salidas n, g, Hn, H.