General Relativity calculations



Hello!

I'd like to do the following computations. There is a mertic, which
has power series of h - constant scalar:

g_{ik}=g0_{ik}+h*g1_{ik}+h^2*g2_{ik}

and I want to calculate some expression by this metric, for example
Riemann tensor and obtain expressions in front of h^n. But I fail in
doing this even for Christoffel symbol

Code:

load(itensor)$
components(g([i,j],[]),g0([i,j],[])+h*g1([i,j],[])+h^2*g2([i,j],[]))$
imetric:g$
ishow(ichr2([i,j],[k]));
ishow(coeff(%,h,2));
ishow(icurvature([i,j,k],[l]));
ishow(coeff(%,h,2));

answer for coeff in ichr2 is 0. 
What is the problem? I should use another function instead coeff?
And why icurvature does not substitute Christoffel symbols?

Semyon.