There are some errors in the file 'MeshDef.m'

[yixincfd]

if norm(deltaXn-deltaX)<=tol && norm(deltaYn-deltaY)<=tol
    ddeltaX_dDeltaX(i,Nc{i}) = k/sum(k);
    ddeltaY_dDeltaY(i,Nc{i}) = k/sum(k);
    break;
else
    deltaX = deltaXn;
    deltaY = deltaYn;
end

It should be:

if norm(deltaXn-deltaX)<=tol && norm(deltaYn-deltaY)<=tol
    for i=1:nNodes
        if isempty(find(i==B,1))    
            k = 1./sqrt((nodes(i,1)-nodes(Nc{i},1)).^2+(nodes(i,2)-nodes(Nc{i},2)).^2);
            ddeltaX_dDeltaX(i,Nc{i}) = k/sum(k);
            ddeltaY_dDeltaY(i,Nc{i}) = k/sum(k);
        end
    end
    break;
else
    deltaX = deltaXn;
    deltaY = deltaYn;
end

[yixin-cfd]

ddeltaX_dDeltaX should be identity matrix