1/* a few sample plots */ 2 3 4numer:true; 5 6/* some nice figures. fig1.ps from fig1.mac 7 shows the units of a quadratic number field 8 as lying on the yx= +- 1 9 */ 10 11viewps("fig1.ps"); 12 13/* REal part of z ^ 1/3 */ 14 15viewps("maxout-3.ps"); 16 17plot2d(sin(x),[-%pi,%pi]); 18plot2d(3*sin(x),[-%pi,%pi]); 19 20block([ps_scale:[40,20],ps_translate:[5,15]], 21 psaxes(4,10), 22 plot2d(3*x^2*sin(x),[-%pi,%pi])); 23 24viewps("fig21.ps"); 25 26/* REal part of z ^ 1/3 */ 27 28block([ps_scale:[200,200], 29 ps_translate:[1.5,1.5], 30 colour_z:true,transform_xy:polar_to_xy], 31 closeps(), 32 plot3d(r^.3333*cos(th/3),[1,1,1.4],[0,1,0,6*%pi],[12,81]), 33 closeps(),viewps()); 34 35 36 37/* REal part of z ^ 1/6 */ 38block([ps_scale:[200,200], 39 ps_translate:[1.5,1.5], 40 colour_z:true,transform_xy:polar_to_xy], 41 plot3d(r^(1/6.0)*cos(th/6),[1,1,1.4],[0,1,0,2*6*%pi],[12,121]), 42 closeps(),viewps() 43 ); 44 45 46 47 48display2d:false; 49set_up_dot_simplifications([z.z+y.y+y.x+x.y,-a*z.z+y.x+x.y+x.x,a*z.y,z.x-a*y.z-x.z],4); 50centrals:fast_central_elements([x,y,z],3); 51centrals4: fast_central_elements([x,y,z],4); 52set_up_dot_simplifications(append(centrals,centrals4,[z.z+y.y+y.x+x.y,-a*z.z+y.x+x.y+x.x,a*z.y,z.x-a*y.z-x.z]),8); 53monomial_dimensions(8); 54