fc06d37Bead-and-spring percussion simulator: numbl-scripted setup, Verlet/Yoshida solver with exact damping, linear modes, browser appJeremy Magland 1function [pos, E] = bp_chain(n, len)
2% BP_CHAIN A straight line of n beads along x, from 0 to len.
3% [pos, E] = bp_chain(n, len) returns positions (n-by-3) and the n-1
4% neighbor pairs E. The end beads are 1 and n.
5x = linspace(0, len, n)';
6pos = [x, zeros(n, 2)];
7E = [(1:n-1)', (2:n)'];
8end