R = compconst (nodes,Nnodes,elements,Nelements,D,C);