我需要做一个牛顿法迭代，即
[at1,bt1]=[at;bt]-([laa,lab;lab,lbb]^(-1))*[da;db];

laa，lab，lbb是二阶导数
但是运行时会提示出现奇异阵：Matrix is singular, close to singular or badly scaled. Results
may be inaccurate. RCOND = NaN.

laa,lab,lbb,da,db都是固定的方程，无法修改，所以不知道怎么办了，求高手帮助下解决。