1 
X 

Z 

ip 

tol 
1 2 3 4 5 6 7 8 9 10 11 12 13 14  ## Should be DIRECTLY executable !! 
## ==> Define data, use random,
## or do help(data=index) for the standard data sets.
## The function is currently defined as
function (X, Z = diag(nrow(X)), ip = diag(nrow(X)), tol = 1e07)
{
help < "\n ConjComp returns a basis for the conjugate complement of the\n conjugate projection of X into span(Z) with respect to inner product with\n matrix ip.\n Note: Z is assumed to be of full column rank but not necessarily X.\n "
xq < qr(t(Z) %*% ip %*% X, tol = tol)
if (xq$rank == 0)
return(Z)
a < qr.Q(xq, complete = T)[, (1:xq$rank)]
Z %*% a
}

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.