x = linspace(0.01, 20,50000);

y = 0.003./sin(1.09*x).^2;
loglog(x,y);
title('8.05 PS#7, P2(d): Parameters consistent with results');
xlabel('\Delta m^2c^4 (in eV^2)');
ylabel('sin^2(2\theta)');
xlim([0.01 20]);
ylim([0.001 2]);

z = ones(1,50000);
hold on;
plot(x,z,'r');