To illustrate the process we shall find
from the equation of the hyperbola
Sage can be made to do a lot of this work for you (though the notation doesn't get any prettier):
[fontsize=\small,fontfamily=courier,fontshape=tt,frame=single,label=\sage]
sage: x = var("x")
sage: y = function("y",x)
sage: a = var("a")
sage: b = var("b")
sage: F = x^2/a^2 - y^2/b^2 - 1
sage: F.diff(x)
2*x/a^2 - 2*y(x)*diff(y(x), x, 1)/b^2
sage: F.diff(x,2)
-2*y(x)*diff(y(x), x, 2)/b^2 - 2*diff(y(x), x, 1)^2/b^2 + 2/a^2
sage: solve(F.diff(x) == 0, diff(y(x), x, 1))
[diff(y(x), x, 1) == b^2*x/(a^2*y(x))]
sage: solve(F.diff(x,2) == 0, diff(y(x), x, 2))
[diff(y(x), x, 2) == (b^2 - a^2*diff(y(x), x, 1)^2)/(a^2*y(x))]
This basically says
david joyner 2008-11-22