Plot the length-tension relationship of a sarcomere

sarcomere_length_tension_curve.R

# Plot length-tension curve of a sarcomere
length = c(1.2, 1.7, 2, 2.2, 3.65)  # µm
force = c(0, 85, 100, 100, 0)  # % of maximum force
xLabel = "Sarcomere Length (µm)"
yLabel = expression("Force (% "*P[0]*")")

svg("sarcomere_length_tension_curve.svg", width = 10, height = 5.5)
par(mar = c(4, 4, 0, 0.5) + 0.25)  # decrease bottom, top and right margin
plot(length, force, xlab = xLabel, ylab = yLabel, type = "l", col = "steelblue",
     bty = "n", xlim = c(0, 4), ylim = c(0, 105), lwd = 2, yaxs = "i", xaxs = "i")
dev.off()