n=-700:1:-300;
x = pi * 10.^(n/100);
y = sin(x);
z = sqrt((1-10^(-7)*round(10^7*cos(x).^2)));
w = (1-cos(x).^2);
q = round(1e7*w)*1e-7;   
for i=1:length(q)
  fprintf('  %.20f %.20f\n', w(i), q(i));
end
