Rotation

> R := matrix( 2, 2, [ [ cos(Pi/6), -sin(Pi/6) ],
[ sin(Pi/6), cos(Pi/6) ] ] );

> Nr := evalm( R &* N );

> pointplot( Nr, opts, view=[-5..10,0..10] );

>