[cooleggyh]

因为对那个ellipse_tri不太了解，所以不是很明白。
不过有一点注意到，就是一般做cell array的时候都是
numrow=length(x);
而不是size(x);
因为如果用你的算法numrow会是一个2-1 matrix，所以那个会多出很多行，因为for i=1:numrow会从全是1的矩阵运行到numrow