R script (22/23)

R script