OpenSCAD

module arm(r, marg) {
union() {
    for(i = [0:1:20])
    {
        translate([rands(-5,5,1)[0],rands(-5,5,1)[0],i*r*marg]) {
            sphere(r);
        }
    } 
    }
}

module branch() {
union() {
    arm(5,0.8); 
    arm(2,0.8); 
    arm(3,0.8); 
}
}

branch(); 
for(i=[0:1:10])
    {
    ang = rands(-45,45,1)[0]; 
    dis = rands(-30,30,1)[0];
    translate([0,0,dis]) {
        rotate([ang,ang,ang]) {
        branch(); 
        }
    }
}