Just in case anyone doesn't want to figure out E[X^2]
here is how I did it.
E[x] = int( x * fX(x),x,a,b)
E[x^2] = int( x^2 * fX(x),x,a,b)
then I just go Var(x) = E[x^2] - (E[x])^2
