% FisherTea.m
D=[3 1;1 3]
[H,P,STATS] = fishertest(D)
[p,cs,expected] = contincy(D, 'chi2')