我使用的代码最初用于流过半圆柱（对称）。我需要更改几何体，以便将半圆替换为半三角形。当我通过修改几何坐标进行此更改时，我得到如下错误：
使用==＆gt;时出错CHOL
矩阵必须是肯定的。
==＆gt;中的错误ellt1 at 35
R = chol（A）; Y = R＆＃39; \ B; Z = R \ Y;
==＆gt;中的错误demo2 at 67
ZN = ellipt1(p,t,RDZ,[],W);
代码的一部分，这个＆＃39; chol＆＃39;使用函数，几何代码如下图所示：
我改变坐标的几何代码如下所示：
http://tinypic.com/r/14nzd3m/8
我如何使其工作或我需要做出哪些修改？