function deriv
% MATAB Plot Code "deriv" for Newton's Quotient versus the reciprocal exponent -log(h):
clc
m1 = 200;  m2 =10;
iv = 1:m1;
for i=1:m1, 
    mloghv(i) = 1+(iv(i)-1)/m2;
    hv(i) = 10^(-mloghv(i)); 
    dv(i) = (2*(1+hv(i))-2)/(hv(i));
end
%
hplot=plot(mloghv,dv,'r-');
axis tight
htitle=title('Newton''s Quotient for Double Precision');
hylabel=ylabel('Q_N');
hxlabel=xlabel('-log(h)');
% Advanced Plot Formatting for Classroom Presenta
set(gca,'Fontsize',28,'FontName','Helvetica','FontWeight','Bold');
set(hplot,'linewidth',4);
set(htitle,'Fontsize',36,'FontName','Helvetica','FontWeight','Bold');
set(hylabel,'Fontsize',32,'FontName','Helvetica','FontWeight','Bold');
set(hxlabel,'Fontsize',32,'FontName','Helvetica','FontWeight','Bold');