quad1.m

有兴趣学系matlab的朋友

function S=quad1(x)
s=0;
k=0;
a=rand(100);
for i=1:100
    for j=1:100
        if j/100<sin(i/100*x)
            s=s+a(i,j);
        end    
    end
end
for l=1:100
    for m=1:100
        k=k+a(l,m);
    end
end
s/k*pi;
end