## Not run:
# ## Loading from the file system
# fn <- system.file("extdata", "1AJC_chainA.pdb", package="Rknots")
# protein <- loadProtein(fn)
# protein <- newKnot(protein$A)
# protein.cp <- closeAndProject(protein)
#
# # Plot the results
# par(mfrow = c(1,2))
# plot(protein, lwd = 2)
# plot(protein.cp, lwd = 2)
# ## End(Not run)
Run the code above in your browser using DataLab