Thanks a lot for the hints. The sapply method may work for me. But how can I extract just the ".grad" expression from the return value of the deriv function? (and secondly store in a matrix of expressions?) Thanks again.