part3

0001 u2 = u1 + L2 * cos( theta + gamma );
0002 v2 = v1 + L2 * sin( theta + gamma );
0003 
0004 u3 = u1 + L3 * cos( theta );
0005 v3 = v1 + L3 * sin( theta );
0006 
0007 plot( [u1 u2 u3 u1], [v1 v2 v3 v1], 'r' ); hold on
0008 plot( [ 0 x1 x2], [ 0  0 y2], 'bo' ); hold on
0009 plot( [ 0 u1 ], [ 0 v1 ], 'b-', [ x2 u2 ], [ y2  v2 ], 'b-', [ x1 u3 ], [ 0 v3 ], 'b-' );