function fern
%ifs code for fern. first 6 columns = M 7th =p'
fern=[ 0 0 0 0.16 0 0 .01;
.85 .04  -0.04   0.85 0 1.6 .85;
.2  -.26 0.23 .22 0 1.6 .07;
 -.15 .28 .26 .24 0  .44 .07];
