function p = set_accel(p,xddot,n)
p.xddot(n,:) = xddot;
end