# Perception probability of visually estimated meteor of magnitude 3.0
# with a limiting magnitude of 5.6.
vmperception(5.6 - 3.0)
# plot
old_par <- par(mfrow = c(1,1))
plot(
vmperception,
-0.5, 8,
main = paste(
'perception probability of',
'visual meteor magnitudes'
),
col = "blue",
xlab = 'm',
ylab = 'p'
)
par(old_par)
Run the code above in your browser using DataLab