& 2698*invtk*is**1.5 + 1776*invtk*is2 +
_RL is2
is2=is*is
& (-12.1652 _d 0*invtk + 0.07871 _d 0) * is2 +
& 2698. _d 0*invtk*is**1.5 _d 0 + 1776. _d 0*invtk*is2 +
& (-12.1652*invtk + 0.07871) * is2 +