Здравствуйте, piero_, Вы писали:

_>int TCENS::Init()
_>{
_> ...
_> V4 Q2(g1, V3(g1, V2(g2, V1(g2,0))));

Здесь ты проинитил локальную переменную Q2.

_>int TCENS::DoOneStep()
_>{
_> ...
_> Q2[i1][i2][i3][i4] += hr; // здесь происходит ошибка при первом проходе т.е. i1=i2=i3=i4=0

А здесь ты обращаешься к непроиниченному (нулевого размера) мемберу Q2.