rmeddis@0: % these data are generated by the testBM program rmeddis@0: rmeddis@0: rmeddis@0: f=[250 500 1000 2000 4000 8000]; rmeddis@0: levels=-10:10:90; rmeddis@0: rmeddis@0: with=[ rmeddis@0: 3.00036e-010 2.9476e-010 3.67369e-010 2.69424e-010 2.61262e-010 7.18365e-011 rmeddis@0: 9.48796e-010 9.32111e-010 1.16172e-009 8.51994e-010 8.26184e-010 2.27167e-010 rmeddis@0: 3.00036e-009 2.9476e-009 3.67369e-009 2.69424e-009 2.61262e-009 7.18365e-010 rmeddis@0: 9.48796e-009 9.32126e-009 1.16172e-008 8.51994e-009 8.26184e-009 2.27167e-009 rmeddis@0: 2.98964e-008 3.30801e-008 4.3799e-008 2.85983e-008 2.12715e-008 7.18365e-009 rmeddis@0: 8.08142e-008 1.0005e-007 1.0854e-007 8.23884e-008 6.41103e-008 1.57936e-008 rmeddis@0: 1.4489e-007 1.70294e-007 1.67989e-007 1.37031e-007 9.77912e-008 3.74565e-008 rmeddis@0: 1.9871e-007 2.56245e-007 2.50429e-007 1.82704e-007 1.2266e-007 5.73054e-008 rmeddis@0: 2.80266e-007 4.15464e-007 4.42073e-007 3.65116e-007 3.24109e-007 6.7619e-008 rmeddis@0: 4.49308e-007 8.26643e-007 1.0501e-006 9.22625e-007 9.46576e-007 2.13766e-007 rmeddis@0: 8.73737e-007 2.1498e-006 2.9006e-006 2.61691e-006 2.90337e-006 6.36048e-007 rmeddis@0: ]; rmeddis@0: rmeddis@0: without=[ rmeddis@0: 1.85179e-009 1.84842e-009 2.29365e-009 1.68885e-009 1.64802e-009 4.53194e-010 rmeddis@0: 5.85587e-009 5.84521e-009 7.25317e-009 5.34061e-009 5.2115e-009 1.43313e-009 rmeddis@0: 1.85179e-008 1.84842e-008 2.39525e-008 1.67348e-008 1.41015e-008 4.53194e-009 rmeddis@0: 5.70059e-008 6.97074e-008 8.049e-008 5.93827e-008 4.75289e-008 1.43313e-008 rmeddis@0: 1.21363e-007 1.42557e-007 1.42791e-007 1.14796e-007 8.59727e-008 2.31062e-008 rmeddis@0: 1.7631e-007 2.1607e-007 2.05558e-007 1.61584e-007 1.09552e-007 5.05732e-008 rmeddis@0: 2.41242e-007 3.44485e-007 3.25688e-007 2.57527e-007 2.17035e-007 6.31655e-008 rmeddis@0: 3.64591e-007 5.90911e-007 7.23077e-007 6.17822e-007 6.09472e-007 1.38217e-007 rmeddis@0: 6.52753e-007 1.44434e-006 1.89854e-006 1.69908e-006 1.84634e-006 4.09033e-007 rmeddis@0: 1.43691e-006 4.03513e-006 5.54045e-006 5.04855e-006 5.73398e-006 1.24307e-006 rmeddis@0: 3.75649e-006 1.20961e-005 1.69372e-005 1.55783e-005 1.7989e-005 3.87293e-006 rmeddis@0: ]; rmeddis@0: rmeddis@0: rmeddis@0: figure(7), clf rmeddis@0: criterion = 1e-7; rmeddis@0: c1=contour(f, levels, with,[criterion criterion]) rmeddis@0: rmeddis@0: hold on rmeddis@0: c2=contour(f, levels, without,[criterion criterion]) rmeddis@0: title(['criterion= ' num2str(criterion)]) rmeddis@0: set(gca,'Xscale','log')