options=optimset('disp','iter','LargeScale','off',...
    'TolFun',0.000000000000000000000001,'MaxIter',100000,'MaxFunEvals',1000000,...
    'Algorithm','trust-region','TolX',0.0000000000000000000000000000001, ...
    'TypicalX', [1e-15, 1e9]);
x0 = [440e-15,1e9];
norm(root2d(x0))
[x, dev] = fsolve(@root2d,[440e-15,2*pi*801e6],options)