# NOT RUN {
## Fetch stucture
pdb <- read.pdb( system.file("examples/1hel.pdb", package="bio3d") )
## Calculate (vibrational) normal modes
modes <- nma(pdb)
## Fluctuations
f <- fluct.nma(modes)
## Fluctuations of first non-trivial mode
f <- fluct.nma(modes, mode.inds=c(7,8))
# }
Run the code above in your browser using DataLab