//Example 13.3//

l=0.154//nm //length of a single bond
n=750;// number of bonds
L=l*sqrt(2*n)
mprintf("L = %f nm",L)
a=109.5;//degree
b=2;//given
Le=2*n*l*sind(a/b)
mprintf("\nLe = %i nm",Le)