[virtuexu]

明显是第二个除号前少了个点，改成点除就可以了，如下：
z=((log((xx+yy)/2))-((xx.*log(xx)-yy.*log(yy))./(xx-yy))+1);