data=rand(10000,40,40);
result=reshape(data,[10000,40,40,1]);
Is there any sense? data (10000,40,40) might considering same for (10000,40,40,1). There might be no difference to implement 1 more dimemtion, which olready have in space vector,
e.g.
>> data=rand(3,3);
>> data(2,3,1,1)
ans =
0.7719

