sousuo.m

数模中的钻井问题的matlab源码

function y=sousuo(x0,chuli,chidu)
a0=real(x0);
b0=imag(x0);
y=[];
for a=a0-chidu/2:.01:a0+chidu/2
    for b=b0-chidu/2:.01:b0+chidu/2
        [c d]=size(find(abs(complex(a,b)-chuli)<=chidu/2));          %ou juli
        y=[y d];
    end
end
y=max(y);