function figbug_example()
t = 0:0.01:200;
q0a = [1, 1, 1, 1];
qa = lsode("dpddpde",q0a,t);
x1a = sin(qa(:,1));
y1a = -cos(qa(:,1));
x2a = (sin(qa(:,1)) + sin(qa(:,2)));
y2a = -(cos(qa(:,1)) + cos(qa(:,2)));
plot(x1a,y1a, "r;m_1;", x2a,y2a, "b;m_2;");
print("-dfig", "-color", "-textspecial", "doublependulumposition_a_q1_1.fig");
endfunction
