`fi=demR[Last/@kr] gamma[Take[k0,{2,-2}],Take[k0,{3,-1}]];
qf=Fq[[Mod[j,4]+1]]=floF[Most@k0,Rest@(Fk[[1]]),Total@Fq,Join[Total@Fin,{0.}],Join[fi dx,{0.}]];
qin=Fin[[Mod[j,4]+1]]=Subtract[fi,b Most@qf] dx;
Fk=RotateLeft[Fk];
k0=Fk[[-1]]=updateF[Join[{0.},Divide[Subtract[Most@qf,Rest@qf]+qin,Vf],{0.}],k0];
qsum=Total@Rq;
qr=Transpose[Join[Transpose[floR[Most/@kr,Rest/@Rk[[1]],Most/@qsum]],{fi}]];
qr[[All,RML]]=compiledMin[RMori,floR[kr\[Transpose][[RML]],Rk[[1]][[All,RML+1]],qsum\[Transpose][[RML]]]];