from PhyNetPy.BirthDeath import CBDP

        # a birth rate of .05, an extinction rate of .01, 
        # 6 extant taxa, and the default 100% sampling rate
        model = CBDP(.05, .01, 6) 
        
        model.generateTree() #to see the network, call .printGraph() on the result
        model.generateTrees(100) #to generate 100 trees

        #Use of the optional sampling rate parameter
        model2 = CBDP(.05, .01, 6, .5) 

        #to generate trees under the Yule model, set extinction rate to 0
        yuleModel = CBDP(.05, 0, 6)