Implement tree error boundaries
Ignore .RData
Add arguments H01 and H02 for subpop
Add multiple nmax 300-500
Delete H0 in argument
Minor modification
Delete H0 in arguments
Add maxiter=200 in uniroot
Add return(bdry)
Modify boundary c when at c_5 stage not leaf
TO DO: consider final analyses more carefully
Compute type 12 errors with different nmax
Add functionality to compute type I and II error