matlab隧道塑性区

求问利用摩尔库伦准测绘制 隧道开挖诱发周围地层塑性区范围的代码 ！

已经求出各个点的大小主应力，怎么画图啊？大小主应力有正有负需要都变成正数吗（用绝对值函数abs)?