秩 = 0，tol = NaN

代码：
if(abs(duj)<10^-10 | abs(duj/delta_disp)<10^-5)
break
end
结果：
In Solution_of_SDOF_under_Earthquake (line 87)
警告: 秩不足，秩 = 0，tol = NaN。