Skip to content

Instantly share code, notes, and snippets.

@soh-i

soh-i/gist:5246952

Last active December 15, 2015 10:39
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save soh-i/5246952 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save soh-i/5246952 to your computer and use it in GitHub Desktop.
Download ZIP
問題9-8 (b)
Raw
gistfile1.r
R <- c(1:10) #[nm]
r <- 0.5 #[nm]
conc.mM <- c(1, 4, 25, 100)
color <- c('pink', 'lightblue', 'lightgreen', 'Orange')
lB <- c(0.757)
U <- c(0)
for (lambD in sqrt(1/(4*pi*0.757*conc.mM))) {
U <- 8*(R/r)^4*(lB*lambD/R*(R+lambD))
plot(U, ylim=c(0,2000), xlim=c(1,3), ylab="U(kBT)", xlab="R[nm]", type="l", col="blue", lwd=2)
par(new=T)
}
legends <- c('1mM', '4mM', '25mM', '100mM')
legend('topleft', legend=legends, pch=20, col='blue', lwd=2, bty='n')
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment