Ecco mostrato il testo dell'algoritmo probabilistico che impiega Matlab. I simboli sottolineati sono dei link.

L = 100

nm = 3

nM = 14

f = inline(’floor((n-1)/2)’)

N = inline(’factorial(n+1)/(factorial(k+1)*factorial(n-k))-1’)

S = inline(’ceil(((N+1)/((n-k)*(k+1)+1))-1)’)

Smax = S(N(f(nM),nM),f(nM),nM)

TA = []

L'output TA dà le dimensioni delle varietà s-secanti alle Grassmanniane G(k,n) e le rispettive difettività, per n fra i due estremi nm e nM .