title "Trapezoidal Rule Area Estimate for P(&LOW<Z<&HIGH)";
title2 "(based on &NPTS equally spaced points)";
var low high incr area_est;
run;
%end;
data trapper2;
set trapper;
array x_value(&npts) x1-x&npts;
array f_value(&npts) y1-y&npts;
do ii=1 to &npts;
xout = x_value[ii];
yout = f_value[ii];
output;
end;
run;
%if %upcase(&print_pts)=TRUE %then %do;
proc print data=trapper2;
title "Interpolation Points for Trapezoidal Rule";
var ii low high incr area_est xout yout;
run;
%end;
%if %upcase(&display ...
Become an O’Reilly member and get unlimited access to this title plus top books and audiobooks from O’Reilly and nearly 200 top publishers, thousands of courses curated by job role, 150+ live events each month, and much more.