import java.applet.Applet;
import java.awt.*;

public class raytracer1 extends Applet{   
   int NSPHERE=5; /* 表示する球の個数 */
   int SIZE=150; /* 表示する画像の大きさ */
   int DEPTH=4; /* 光線追跡の深さ */
   Button button1;
   TextField tf1[],tf2[],tf3[],tf4[];
   TextField tf6,tf7,tf8,tf9,tf10,tf11,tf0;
   Choice paramchoice,paramchoice2;
   Label depthLabel,noLabel,paramLabel,viewLabel,ambLabel;
   Panel choicePanel,inputPanel,depthPanel,noPanel,butnPanel,canvasPanel,
         matelPanel,paramPanel0,viewPanel,ambPanel,dispPanel,allPanel;
   Panel paramPanel[];
   Phong phongCanvas;
   DispCanvas disp; 
   double r_view=10.,the_view=0.0,pi_view=0.0,
          amb1=0.1,amb2=0.1,amb3=0.8; /* 視点、環境光の初期値 */
   double ks[]={0.99,0.1,0.6,0.6,0.6};double kt[]={0.0,0.9,0.0,0.0,0.0};
   double ir[]={1.0,3.0,1.0,1.0,1.0};double deg[]={60.,60.,60.,60.,60.}; /* 球の属性の初期値 */

   public void init(){
      int nn=1;
      String A="0.0";
      System.out.println("test");
      phongCanvas = new Phong(SIZE);
      disp = new DispCanvas();      
      phongCanvas.setNSPHERE(NSPHERE);
      //phongCanvas.setSIZE(SIZE);
      phongCanvas.setDEPTH(SIZE);
      phongCanvas.setparam(ks,kt,ir,deg);
      System.out.println("test2");
      phongCanvas.setView(r_view,the_view,pi_view);
      phongCanvas.setAmb(amb1,amb2,amb3);
      disp.setNSPHERE(NSPHERE);
      disp.setparam(ks,kt,ir,deg);

      choicePanel = new Panel();allPanel= new Panel();                 
      noPanel = new Panel();inputPanel = new Panel();matelPanel  = new Panel();
      paramPanel0 = new Panel();viewPanel = new Panel();ambPanel = new Panel();
      butnPanel = new Panel();canvasPanel = new Panel();  
      dispPanel = new Panel();depthPanel = new Panel();
      
      button1 = new Button("trace");
      paramchoice = new Choice();paramchoice2 = new Choice();
   
      tf0 = new TextField(A.valueOf(DEPTH),1);      
      tf1 = new TextField[NSPHERE];tf2 = new TextField[NSPHERE];
      tf3 = new TextField[NSPHERE];tf4 = new TextField[NSPHERE];
      for(int i=0;i<NSPHERE;i++){
         tf1[i] = new TextField(A.valueOf(deg[i]),4);
         tf2[i] = new TextField(A.valueOf(ks[i]),4);
         tf3[i] = new TextField(A.valueOf(kt[i]),4);
         tf4[i] = new TextField(A.valueOf(ir[i]),4);
      }
      tf6 = new TextField(A.valueOf(r_view),4);
      tf7 = new TextField(A.valueOf(the_view),4);
      tf8 = new TextField(A.valueOf(pi_view),4);
      tf9 = new TextField(A.valueOf(amb1),4);
      tf10 = new TextField(A.valueOf(amb2),4);
      tf11 = new TextField(A.valueOf(amb3),4);
      
      depthLabel= new Label("depth :");
      noLabel =   new Label("sphere No. :");
      viewLabel = new Label("view point  r,th,pi :");
      ambLabel =  new Label("back color  r,g,b:");
      
      paramchoice.addItem("depth");
      paramchoice.addItem("material");
      paramchoice.addItem("view point");
      paramchoice.addItem("ambient color");
      choicePanel.setLayout(new FlowLayout(FlowLayout.LEFT));
      choicePanel.add(paramchoice);
            
      depthPanel.setLayout(new FlowLayout(FlowLayout.LEFT));
      depthPanel.add(depthLabel);
      depthPanel.add(tf0);      

      paramchoice2.addItem("0:red");paramchoice2.addItem("1:transparent");
      paramchoice2.addItem("2:green");paramchoice2.addItem("3:blue");
      paramchoice2.addItem("4:yellow");
      
      noPanel.setLayout(new FlowLayout(FlowLayout.LEFT));
      noPanel.add(noLabel);
      noPanel.add(paramchoice2);

      matelPanel.setLayout(new GridLayout(2,1));
      matelPanel.add(noPanel);      

      paramPanel = new Panel[NSPHERE];
      for(int i=0;i<NSPHERE;i++){
         paramPanel[i] = new Panel();
      }
      
      for (int i=0;i<NSPHERE;i++){                  
         paramLabel = new Label("n,W,kt,ir :");         
         paramPanel[i].add(paramLabel);
         paramPanel[i].add(tf1[i]);paramPanel[i].add(tf2[i]);
         paramPanel[i].add(tf3[i]);paramPanel[i].add(tf4[i]);        
      }
      paramPanel0.setLayout(new CardLayout());
      paramPanel0.add("0:red",paramPanel[0]);
      paramPanel0.add("1:transparent",paramPanel[1]);
      paramPanel0.add("2:green",paramPanel[2]);
      paramPanel0.add("3:blue",paramPanel[3]);
      paramPanel0.add("4:yellow",paramPanel[4]);;
      matelPanel.add(paramPanel0);
          
      viewPanel.setLayout(new FlowLayout(FlowLayout.LEFT));
      viewPanel.add(viewLabel);viewPanel.add(tf6);
      viewPanel.add(tf7);viewPanel.add(tf8);
      ambPanel.setLayout(new FlowLayout(FlowLayout.LEFT));
      ambPanel.add(ambLabel);ambPanel.add(tf9);
      ambPanel.add(tf10);ambPanel.add(tf11);
      butnPanel.setLayout(new FlowLayout(FlowLayout.RIGHT));  
      butnPanel.add(button1);
      
      inputPanel.setLayout(new CardLayout());
      inputPanel.add("depth",depthPanel);
	   inputPanel.add("material",matelPanel);
      inputPanel.add("view point",viewPanel);
	   inputPanel.add("ambient color",ambPanel);
	  
      //allPanel.setLayout(new GridLayout(3,1));
      //allPanel.add(choicePanel);inputPanel.add(inputPanel);
      //allPanel.add(butnPanel);
      
      canvasPanel.setLayout(new FlowLayout(FlowLayout.LEFT));
      phongCanvas.resize(SIZE,SIZE);
      phongCanvas.setBackground(Color.black);    
      canvasPanel.add(phongCanvas);
 
      disp.resize(130,100);
      disp.setBackground(Color.white);
      dispPanel.setLayout(new FlowLayout(FlowLayout.RIGHT));
      dispPanel.add(disp);     
      
      add(canvasPanel);
      add(dispPanel);
      //add(allPanel);
      add(choicePanel);
      add(inputPanel);
      add(butnPanel);
      resize(310,330);
      setBackground(Color.gray); 
      phongCanvas.start();
 }
 
   public boolean action(Event e,Object o) {
      int nn;
      String label = o.toString();
      Double A,B,C,D,Deg,Ks,Kt,Ir,Vr,Vt,Vp,Ambr,Ambg,Ambb;
      
      if(e.target instanceof Choice){
            ((CardLayout)inputPanel.getLayout()).show(inputPanel,label);
            ((CardLayout)paramPanel0.getLayout()).show(paramPanel0,label); 
      }
           
      if(e.target instanceof Button){
         System.out.println("trace");
         for(int i=0;i<NSPHERE;i++){
            A = new Double(0.);B = new Double(0.);C = new Double(0.);
            D = new Double(0.);Deg = new Double(0.);Ks = new Double(0.);             Kt = new Double(0.);Ir = new Double(0.);
            Deg = A.valueOf(tf1[i].getText());Ks = B.valueOf(tf2[i].getText());
            Kt = C.valueOf(tf3[i].getText());Ir = D.valueOf(tf4[i].getText());  
            ks[i]=Ks.doubleValue();kt[i]=Kt.doubleValue();
            ir[i]=Ir.doubleValue();deg[i]=Deg.doubleValue();
         }
         A = new Double(0.);B=new Double(0.);
         B = A.valueOf(tf0.getText());DEPTH=(int)B.doubleValue();
         if(DEPTH>0){
              phongCanvas.setDEPTH(DEPTH);
         }else{
              phongCanvas.setDEPTH(4);
         }       
         phongCanvas.setparam(ks,kt,ir,deg);
         disp.setparam(ks,kt,ir,deg);
         A = new Double(0.);B = new Double(0.);
         C = new Double(0.);         
         Vr = new Double(0.);Vt = new Double(0.);Vp = new Double(0.);
         Vr = A.valueOf(tf6.getText());Vt = B.valueOf(tf7.getText());
         Vp = C.valueOf(tf8.getText());
         phongCanvas.setView(Vr.doubleValue(),Vt.doubleValue(),Vp.doubleValue());
         A = new Double(0.);B = new Double(0.);
         C = new Double(0.);         
         Ambr = new Double(0.);Ambg = new Double(0.);Ambb = new Double(0.);
         Ambr = A.valueOf(tf9.getText());Ambg = B.valueOf(tf10.getText());
         Ambb = C.valueOf(tf11.getText());
         phongCanvas.setAmb(Ambr.doubleValue(),Ambg.doubleValue(),                 Ambb.doubleValue());
         
         disp.repaint();
         phongCanvas.start();
      }
      
   return true;
   }
}

class DispCanvas extends Canvas {  /* アプレットの右上にパラメータの情報を表示  */
   double ks[] = new double[10];double kt[] = new double[10];
   double ir[] = new double[10];double deg[] = new double[10];
   int NSPHERE;
   
   public void setNSPHERE(int i){ /* 球の個数を代入 */
      NSPHERE=i;
   }

   public void setparam(double a[],double bb[],double c[],double d[]){ /* 球の属性を表わす変数の代入 */
      for(int i =0;i<NSPHERE;i++){
         ks[i]=a[i];kt[i]=bb[i];ir[i]=c[i];deg[i]=d[i];
      }   
   }

   public void paint(Graphics g){ /* 各球の属性の値を表示するための文字列の描画 */
      String str_ks=" ",str_kt=" ",str_ir=" ",str_deg=" ",str_num=" ",str=" ";
      int i,y;
      
      g.setFont(new Font("TimerRoman",Font.PLAIN,12));
      g.drawString("sphere No.=n,W,kt,ir",0,10);    
      y=24;
      for(i=0;i<NSPHERE;i++){
         str_num=str.valueOf(i);
         str_ks=str.valueOf(ks[i]);str_kt=str.valueOf(kt[i]);
         str_ir=str.valueOf(ir[i]);str_deg=str.valueOf(deg[i]);
         g.drawString(str_num+"="+str_deg+","+str_ks+","+str_kt+","+str_ir,0,y);
         y+=14;
      }
   }

}
      
   

class Phong extends Canvas implements Runnable {  /* レイトレーシングを行うクラス */
   int AOV=35; /* 視野角 */

   double TOL=1e-7; /* 交差判定を行うときの許容誤差 */
   double u,b,tmin,xmax,xmin;

   Vec U = new Vec();
   Vec black = new Vec(0.,0.,0.);
   Vec black2 = new Vec(0.,0.,0.);
   Vec amb = new Vec(0.5,0.5,0.9);
   Vec light = new Vec(1.0,0.0,1.0);
   double ambient=0.15;
   double r_view=15.0,the_view=0.0,pi_view=0.0; 
   double xv=0.0,yv=0.0,zv=0.0;
   double ct,st,cf,sf;
   int SIZE,flug,NSPHERE,DEPTH;
   Sphere sph[];
   double ks[]=new double[10];double kt[]=new double[10];
   double ir[]=new double[10];double deg[]=new double[10];

   Image doubleBuffer;
   Graphics myGC;
   Thread kick=null;   

   Phong(int a){
      SIZE = a;
   }

   public void setNSPHERE(int i){ /* 球の個数を代入　*/
      NSPHERE=i;
   }

   public void setDEPTH(int i){ /* 追跡の深さを代入 */
      DEPTH=i;
   }
      
   //public void setSIZE(int w){  /* 画像のサイズを代入 */
     // SIZE=w;
   //}

   public void setparam(double a[],double bb[],double c[],double d[]){ /* 球の属性を代入 */
      for(int i =0;i<NSPHERE;i++){
         ks[i]=a[i];kt[i]=bb[i];ir[i]=c[i];deg[i]=d[i];
      }   
   }

   public void setView(double r,double t,double p){ /*　視点に関する情報を代入 */
      r_view=r;the_view=t;pi_view=p;
   }
   
   public void setAmb(double r,double g,double bb){  /*　環境光の色を代入 */
      amb.setval(r,g,bb);
      System.out.println("amb = "+amb.x+" "+amb.y+" "+amb.z);
   }
   
   public void start(){
      if(kick == null){
         doubleBuffer = createImage(SIZE,SIZE);
         myGC = doubleBuffer.getGraphics();
         kick = new Thread(this);
         kick.start();
      }    
   }
   
   public void stop(){
      if(kick != null){
         kick.stop();
         kick = null;
      }
   }
   
   public void update(Graphics g){
       g.drawImage(doubleBuffer,0,0,this);
   }
   
   public void paint(Graphics g){
       update(g);
   }

   public void run(){  /* 描画メソッド */
      int x,y,i,ysn;
      double fa,fb,fc,fd,scale,txw;
      sph = new Sphere[NSPHERE];
      sph[0] = new Sphere(  0.,  0.,-1.0, 1.,0.1, .1,1.5,
               0.8,ks[0],kt[0],0.,ir[0],deg[0]);
      sph[1] = new Sphere( 4.5,-0.5,-1.0, .0, .1, .3,1.0,
               0.8,ks[1],kt[1],0.,ir[1],deg[1]); 
      sph[2] = new Sphere( 2.5, 2.5,-0.5, .0, .8, .0,1.0,
               0.8,ks[2],kt[2],0.,ir[2],deg[2]);       
      sph[3] = new Sphere(-2.5,-2.5, 0.4, .4, .9, .8,1.0,
               0.8,ks[3],kt[3],0.,ir[3],deg[3]);
      sph[4] = new Sphere(-2.5, 2.5, 1.5, .8, .8, .0,1.0,
               0.8,ks[4],kt[4],0.,ir[4],deg[4]);
      Vec COL = new Vec();
      Vec AA = new Vec();
      Vec BB = new Vec(); 
      Vec CC = new Vec();

      myGC.clearRect(0,0,SIZE-1,SIZE-1);
      light.vunit(light);
      System.out.println("light="+light.x+" "+light.y+" "+light.z);
      transfrm();
      scale=(double)SIZE/2/Math.tan(AOV/114.5915590261)/r_view;
      txw=r_view*Math.tan(AOV/114.5915590261); 

      for(y=0;y<=SIZE;y++){
         ysn=yscan(((float)SIZE/2-y)/scale,scale,txw); 
         //System.out.println(" ");
         //System.out.println("ysn="+ysn+" "+xmin+" "+xmax);
         //ysn=1;xmin=0;xmax=SIZE;
         if(ysn!=0){ 
            for(x=0;x<(int)xmin;x++){
               COL.x=255.*amb.x;COL.y=255.*amb.y;COL.z=255.*amb.z;
               myGC.setColor(new Color((int)COL.x,(int)COL.y,(int)COL.z));
               myGC.drawLine(x,y,x,y);
               repaint();
            }
            for(x=(int)xmin;x<=(int)xmax;x++){
               flug = 0;
               fa=(x-SIZE/2)/scale;
               fb=(SIZE/2-y)/scale;
               fc=r_view;
               U.setval(fa,fb,fc);
               AA.vunit(U);
               BB=trace(DEPTH,black,AA);
               COL.x=255.*BB.x;COL.y=255.*BB.y;COL.z=255.*BB.z;
               if(COL.x<0) COL.x=0.;if(COL.x>255) COL.x=255.;
               if(COL.y<0) COL.y=0.;if(COL.y>255) COL.y=255.;
               if(COL.z<0) COL.z=0.;if(COL.z>255) COL.z=255.;
               myGC.setColor(new Color((int)COL.x,(int)COL.y,(int)COL.z));
               myGC.drawLine(x,y,x,y);
               repaint();
            }
            for(x=(int)(xmax+1);x<SIZE;x++){
               COL.x=255.*amb.x;COL.y=255.*amb.y;COL.z=255.*amb.z;
               myGC.setColor(new Color((int)COL.x,(int)COL.y,(int)COL.z));
               myGC.drawLine(x,y,x,y);
               repaint();
            }
         }else{
            for(x=0;x<SIZE;x++){
               COL.x=255.*amb.x;COL.y=255.*amb.y;COL.z=255.*amb.z;
               myGC.setColor(new Color((int)COL.x,(int)COL.y,(int)COL.z));
               myGC.drawLine(x,y,x,y);
               repaint();
            }
         }
      }
      kick=null;
   }
   
   int yscan(double ys,double scale,double txw){  /* 描画を高速化するため、１スキャンラインごとに球の存在するｘの範囲を求める */
      int i,num;
      double ddy2,ddy,day,ya,za,s1,s2;
      double d1,dd, d2, x1, x2, xz,r2,zz,zd;

      ddy2 = r_view*r_view+ys*ys;
      ddy = Math.sqrt(ddy2);
      ya = r_view/ddy; za = ys/ddy;
      xmin = 1000.; xmax = -1000.;
      //System.out.println("txw="+txw+" ys="+ys+" "+ddy+" "+ya+" "+za);
      num=0;
      for (i=0; i<NSPHERE; i++) {
         day = sph[i].cen.y*ya - sph[i].cen.z*za;
         //System.out.println("Math.abs(day)"+Math.abs(day)+" "+sph[i].rad);
         if(Math.abs(day) > sph[i].rad) {
         s1 = 0.;
         }else {
         s1 = 1.;
         num++;
         r2 = sph[i].rad*sph[i].rad; 
         d1 = sph[i].cen.z*sph[i].cen.z - r2;
         dd = sph[i].cen.x*sph[i].cen.x+ d1;
         xz = sph[i].cen.x*sph[i].cen.z;
         d2 = r2-day*day;
         zz = sph[i].cen.z + day*za; zd = zz*zz - d2;
         dd = Math.sqrt((sph[i].cen.x*sph[i].cen.x + zd)*d2);
         x1 = (xz - dd)/zd*r_view; x2 = (xz + dd)/zd*r_view;
         s1 = (x1 + txw)*scale - 1.;
         s2 = (x2 + txw)*scale - 1.;
         //System.out.println("d1="+d1+" "+dd+" "+xz+" "+d2);
         //System.out.println("x1="+x1+" "+x2+" "+s1+" "+s2);
         if(s1 < xmin) xmin = s1;
         if(s2 > xmax) xmax = s2; /* ysでのスキャンライン上には、xminからxmaxまでの間に球が存在している */
         }
      }
      return num;
   }

   void transfrm(){  /* 視線がｚ軸上になるように座標変換する　*/
      double fef,tef,xy,px,py,pz;
      int s;
      
      fef=3.1415/180.*pi_view;tef=3.1415/180.*the_view;
      ct=Math.cos(tef);st=Math.sin(tef);
      cf=Math.cos(fef);sf=Math.sin(fef);
      System.out.println("fef="+fef+" "+tef+" "+ct+" "+st);
      for(s=0;s<NSPHERE;s++){
         xy=ct*sph[s].cen.x+st*sph[s].cen.y;
         pz=r_view-cf*xy-sf*sph[s].cen.z;
         px=ct*sph[s].cen.y-st*sph[s].cen.x;
         py=cf*sph[s].cen.z-sf*xy;
         sph[s].cen.x=px;sph[s].cen.y=py;sph[s].cen.z=pz;
         System.out.println("px,py,pz="+px+" "+py+" "+pz);
      }
   }
   
   double vdot(Vec A,Vec B){ /* ベクトルＡとベクトルＢとの内積関数 */
      return A.x*B.x+A.y*B.y+A.z*B.z;
   }

   int intersect(Vec P,Vec D){ /* 球と視線との交差判定 */

      int best=-1;
      tmin=1e30;
      //System.out.println("NSPHERE= "+NSPHERE);
      for(int i=0;i<NSPHERE;i++){
         //b=vdot(D,U=vcomb(-1,P,sph[i].cen));
         U.vcomb(-1,P,sph[i].cen);
         b=vdot(D,U); 
         u=b*b-vdot(U,U)+sph[i].rad*sph[i].rad;
         // System.out.println("tmin,b,u"+tmin+" "+b+" "+u+" "+TOL);

         u = u>0 ? Math.sqrt(u) : 1e31;
         u = b-u>TOL ? b-u : b+u;
         if((u>=TOL)&&(u<tmin)){
            best = i;
            tmin = u;
         }else{
            tmin=tmin;
         }
      }
      return best;
   }

    Vec trace(int level,Vec P,Vec D){ /* 視線を追跡する　levelが０になるまで再帰的に交点を追跡 */
      int s;
      Vec R ;
      Vec color2;Vec color;Vec color3;
      double Px,Py,Pz,Dx,Dy,Dz,Nx,Ny,Nz,gamma,cos_alfa,Hx,Hy,Hz,cx,cy,cz;
      double fd,specular,diffuse,lightx,lighty,lightz,Rx,Ry,Rz;
      
      if(level -- <= 0) return black;
      s = intersect(P,D);
      if((level == DEPTH-1)&&(s == -1)) return amb;
      if((flug == 1)&&(s == -1)) return amb;
      if((flug == 1)&&(s != -1)) flug = 0;
      if(s == -1) return black;
      
      Px=P.x;Py=P.y;Pz=P.z;Dx=D.x;Dy=D.y;Dz=D.z;
      lightx=light.x;lighty=light.y;lightz=light.z;
      Px=tmin*Dx+Px;Py=tmin*Dy+Py;Pz=tmin*Dz+Pz;
      Nx=-Px+sph[s].cen.x;Ny=-Py+sph[s].cen.y;Nz=-Pz+sph[s].cen.z;
      fd=Math.sqrt(Nx*Nx+Ny*Ny+Nz*Nz);
      Nx=Nx/fd;Ny=Ny/fd;Nz=Nz/fd;
      cos_alfa = lightx*Nx + lighty*Ny+lightz*Nz;
      specular=0.;
      if(cos_alfa>0.){
         diffuse = cos_alfa*sph[s].kd + ambient;
         Hx=lightx+Dx;Hy=lighty+Dy;Hz=lightz+Dz;
         gamma = (Hx*Nx + Hy*Ny + Hz*Nz)/Math.sqrt(Hx*Hx + Hy*Hy + Hz*Hz); 
         if(gamma>0.) specular = sph[s].ks*Math.pow(gamma,(int)sph[s].deg); 
      }else{
         diffuse=ambient;
      }  
      fd=Dx*Nx+Dy*Ny+Dz*Nz;
      fd=(double)Math.abs(fd);
      cx=Dx/fd;cy=Dy/fd;cz=Dy/fd;
      Rx=-(cx+2*Nx);Ry=-(cy+2*Ny);Rz=-(cz+2*Nz);
      R=new Vec(Rx,Ry,Rz);R.vunit(R);
      P=new Vec(Px,Py,Pz);
      color2=new Vec();
      color2=trace(level,P,R);
      color3=new Vec();
      color3=black;
      
      if(sph[s].kt>0.0){
          flug = 1;
          Px = cx*cx+cy*cy+cz*cz;
          Py = (cx+Nx)*(cx+Nx)+(cy+Ny)*(cy+Ny)+(cz+Nz)*(cz+Nz);
          Pz =sph[s].ir*sph[s].ir*Px-Py;
          if(Pz > 0.0){          
             Rx=1./Math.sqrt(Pz)*fd*(Nx+cx)-Nx;Ry=fd*(Ny+cy)-Ny;Rz=fd*(Nz+cz)-Nz;
             R=new Vec(Rx,Ry,Rz);R.vunit(R);
          }else{
             R=new Vec();R=D;

          }
          color3=trace(level,P,R);
      }                

      cx=sph[s].color.x*diffuse + specular + sph[s].ks*color2.x + 
         sph[s].kt*color3.x;
      cy=sph[s].color.y*diffuse + specular + sph[s].ks*color2.y +
         sph[s].kt*color3.y;
      cz=sph[s].color.z*diffuse + specular + sph[s].ks*color2.z +
         sph[s].kt*color3.z;

      return (new Vec(cx,cy,cz));
  } 


}


class Vec{ /* ベクトルクラス */
   public double x,y,z;

   Vec(){
   }

   Vec(double nwx,double nwy,double nwz){
      x=nwx;
      y=nwy;
      z=nwz;
   }

   void setval(double nwx,double nwy,double nwz){ /* ベクトルの値を代入 */
      x=nwx;
      y=nwy;
      z=nwz;
   }  

   void vcomb(double a,Vec A,Vec B){ /* ベクトルＡとベクトルＢとの１次結合関数 */
      x=B.x + a*A.x;
      y=B.y + a*A.y;
      z=B.z + a*A.z;
   }

   void vunit(Vec A){ /* ベクトルＡの正規化　*/
      double a,v;

      v=A.x*A.x+A.y*A.y+A.z*A.z;
      a=1/Math.sqrt(v);
      x=a*A.x;
      y=a*A.y;
      z=a*A.z;      
   }


}

class Sphere{ /* 球のクラス */
   public Vec cen,color;
   public double x,y,z,r,g,b,rad,kd,ks,kt,kl,ir,deg;

   Sphere(){
   }

   Sphere(double ocenx,double oceny,double ocenz,double ocolr,
      double ocolg,double ocolb,double orad,double okd,
      double oks,double okt,double okl,double oir,double odeg){
      cen = new Vec(ocenx,oceny,ocenz);
      color = new Vec(ocolr,ocolg,ocolb);
      rad=orad;kd=okd;ks=oks;kt=okt;kl=okl;ir=oir;deg=odeg;
   }
   
   void setval(double ocenx,double oceny,double ocenz,double ocolr,
      double ocolg,double ocolb,double orad,double okd,
      double oks,double okt,double okl,double oir,double odeg){ /* 球の変数を代入 */
      cen = new Vec(ocenx,oceny,ocenz);
      color = new Vec(ocolr,ocolg,ocolb);
      rad=orad;kd=okd;ks=oks;kt=okt;kl=okl;ir=oir;deg=odeg;
   }

}