supplbasis <-
function(A,B,tol=sqrt(.Machine$double.eps)) {
if (qr(A,tol=tol)$rank < qr(cbind(A,B),tol=tol)$rank) return(basis(B-proj(A,tol=tol)%*%B))
else return(NULL)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.