// SuSE Version 8.1 professional edition
//
// title image
//
// pro_81.pic
//
// A D_12-symmetric surface with 12*11=131 A_3-singularities.
// 

int width_fac=400; width=width_fac; height=width_fac; antialiasing=6; 

illumination = ambient_light + diffuse_light + reflected_light + transmitted_light;
ambient      = 60;
diffuse      = 60;
reflected    = 80;
transmitted  = 20;
smoothness   = 50;
transparence = 0;
thickness = 10;

background_red=255; background_green=255; background_blue=255;
surface_red=82; surface_green=193; surface_blue=71;
inside_red=112; inside_green=203; inside_blue=101;


// from front right up
light1_x=10;
light1_y=10;
light1_z=10;
light1_vol=60;

// from right in the front
light2_x=10;
light2_y=0;
light2_z=5;
light2_vol=10;
light2_red=255;
light2_green=255;
light2_blue=255;

// from top    
light3_x = 0;
    light3_y = 10;
    light3_z = -5;
    light3_vol = 60;
    light3_red   = 255;
    light3_green = 255;
    light3_blue  = 255;

// back left
  light4_x = -10;
    light4_y = 3;
    light4_z = -10;
    light4_vol = 60;
    light4_red   = 255;
    light4_green = 255;
    light4_blue  = 255;

// back right
  light5_x = 10;
    light5_y = 3;
    light5_z = -10;
    light5_vol = 60;
    light5_red   = 255;
    light5_green = 255;
    light5_blue  = 255;

light6_vol = 0;
light7_vol = 0;
light8_vol = 30;
light9_vol = 0;

double factor=0.4;
scale_x=factor*1.0;
scale_y=factor*1.0;
scale_z=factor*1.0;

rot_x=1.77;
rot_y=0.55;
rot_z=0.55;


radius=9.5;
double  Pi=2*arccos(0);
double  c=sqrt(3.38-sqrt(7))/2;
double  a=-(1.577)*(1+1/sqrt(7))*sqrt(7-sqrt(7));
double  r=(1-0.85*sqrt(7))/6;
double d=1;

a=-90;
//c=0.5;
r=-0.2;
d=1;

poly q=((x^2+y^2)+r*z^2-d)^6;

double k=12.0;
poly p1=x-z;
poly p2=cos(2*Pi/k)*x-sin(2*Pi/k)*y-z;
poly p3=cos(4*Pi/k)*x-sin(4*Pi/k)*y-z;
poly p4=cos(6*Pi/k)*x-sin(6*Pi/k)*y-z;
poly p5=cos(8*Pi/k)*x-sin(8*Pi/k)*y-z;
poly p6=cos(10*Pi/k)*x-sin(10*Pi/k)*y-z;
poly p7=cos(12*Pi/k)*x-sin(12*Pi/k)*y-z;
poly p8=cos(14*Pi/k)*x-sin(14*Pi/k)*y-z;
poly p9=cos(16*Pi/k)*x-sin(16*Pi/k)*y-z;
poly p10=cos(18*Pi/k)*x-sin(18*Pi/k)*y-z;
poly p11=cos(20*Pi/k)*x-sin(20*Pi/k)*y-z;
poly p12=cos(22*Pi/k)*x-sin(22*Pi/k)*y-z;

surface=(-q+a*p1*p2*p3*p4*p5*p6*p7*p8*p9*p10*p11*p12);
//surface2=c;
//transparence2=0;
//surface2=x^2+y^2+(z+3)^2-0.2;
draw_surface;

//color_file_format=sun; 
//filename="pro_81.ras";
color_file_format = jpeg;
filename = "pro_81.jpg";
save_color_image;