thank you i was missing the "symbols" command at the beginning. Now it does work 2008/5/15 Przemek Klosowski <przemek at jazz.ncnr.nist.gov>: > symbols > x=sym("x") > f=x^2+3*x > differentiate(f,x) >