http://axcas.net I'm working on ode/pde/bvp solvers this week. I hope to use them to evaluate special functions.