矩阵必须是正定的

Question

我使用的代码最初用于流过半圆柱（对称）。我需要更改几何体，以便将半圆替换为半三角形。当我通过修改几何坐标进行此更改时，我得到如下错误：

使用==＆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/ka4p6f/8

我改变坐标的几何代码如下所示：

http://tinypic.com/r/14nzd3m/8

我如何使其工作或我需要做出哪些修改？