interface(echo=4,screenwidth=72);
N := 9; # Change N to see a table of Eulerian numbers to order N.
en := array(1..N,0..N-1):

en[1,0] := 1:
lprint(en[1,0]);
for i to N-1 do
 en[i+1,0] := 1;
 en[i+1,i] := 1;
 for j to i-1 do
  en[i+1,j] := (j+1)*en[i,j] + (i+1-j)*en[i,j-1];
 od;
 lprint(seq(en[i+1,j],j=0..i));
od:

### For use by authors only
### Experimental editor for the Maple Form Interface
###