%
% example stiff ODE on page 340 in the text.
%
Intv = [0;1];
alpha=-1;
N=5;
M=10;
tol=1e-6;
[w,t,flg] = TrapNewton(@ty4,@ty4dy,Intv, alpha, N, M, tol);