Matlab代码错误请教

ds是dx*dy/夹角的余弦，算得结果是1/9pi，照理说这个思路应该是可以运行的，有没有人能指点一下是哪一步错了