Problem in double integration

syms x;

syms y;

w0=0.46*10^-6;

a = 100*10^-6;

w1 = 0.401;

w2 = 1.611;

d0= 1*10^-6;

e0 = 8.854*10^(-12);

xmin = 0;

xmax = 200*10^-6;

ymin = 0;

ymax = 200*10^-6;

c = @(x,y) 1./(d0 - (w0 + w1.*(x.^2 + y.^2)/a.^2 + w2.*x.^2.*y.^2./a.^4) .* cos(pi.*x./(2.*a)).^2 .* cos(pi.*y./(2.*a)).^2)

quad2d(c,xmin,xmax,ymin,ymax)

I'm trying this code and getting errors like "Reached the maximum number of function evaluations (100000). The result fails the global error test.". the i tried the following options

'Abstol',1e-4,'reltol',1e-4,'Singular',true, 'FailurePlot',true,'MaxFunEvals', 100000

and got the answer in the order of 10^-5. But the value i have to get should be in the range of 10^-15. What might be the problem