[暖阳三杯]

这是函数中的命令：
function PlotStress(iStress)
%显示应力云图
%iStress----应力分量指示，可以为：
%1----x方向正应力
%2----y方向正应力
%3----剪应力
%4----第一主应力
%5----第二主应力
%6----最大剪应力
global gNode gElement gNodeStress

switch iStress
case 1
title='x方向正应力';
case 2
title='y方向正应力';
case 3
title='剪应力';
case 4
title='最大主应力';
case 5
title='最小主应力';
case 6
title='最大剪应力';
end
figure;
axis equal;
axis off;
set(gcf,'NumberTitle','off');
set(gcf,'Name',title);

Minstress=gNodeStress(:,iStress)
Maxstress=max( gNodeStress(:,iStress));
caxis([Minstress,Maxstress]);
colormap('gray');

当运行PlotStress(4)
则出现下列错误提示：
??? Attempted to access gNodeStress(:,4); index out of bounds because size(gNodeStress)=[0,0].

Error in ==> PlotStress at 32
Minstress=gNodeStress(:,iStress)
不明白哪儿出问题了，请高手帮忙看看，谢谢啊！！