% hw020.m
% John Burkardt
% Count derangements.
%
d1 = 0;
d2 = 1;
d3 = 0;

for n = 2 : 10

  d1 = d2;
  d2 = d3;
  d3 = ( n - 1 ) * ( d2 + d1 );

  fprintf ( '  D(%d) = %d\n', n, d3 );

end
