for i=1:1:5724 D(i)= sqrt((contour(:,2)-579).^2 + (contour(:,1)-603).^2); end

这段代码中存在语法错误，应该修改为：

for i = 1:1:5724
    D(i) = sqrt((contour(i,2)-579).^2 + (contour(i,1)-603).^2);
end

在计算每个点到坐标 $(603, 579)$ 的距离时，应该使用当前迭代变量 $i$ 所代表的点的坐标，即 `contour(i,1)` 和 `contour(i,2)`。此外，计算平方和时应该使用加号进行运算。