import java.applet.*;
import java.awt.*;
/*
 applet code="Vortex" width=600 height=400>
  </applet>
*/

public class heat extends Applet 
implements Runnable {
  int xp = 0;
  int yp = 0;

  Thread t;
  Image buffer;
  Graphics bufferg;

  public void init() {

    // スレッドを開始する
    t = new Thread(this);
    t.start();

    // バッファを作成する
    Dimension d = getSize();
    buffer = createImage(d.width, d.height);
  }


  public void run() {
    try {
        while(true) {// 再描画を要求する
        repaint();

        // ウィンドウを更新する前に休止する
        Thread.sleep(1000);
      }
	}
    catch(Exception e) {
    }
  }

public void update(Graphics g) {
    paint(g);
  }

  public void paint(Graphics g) {

	int NI=21;
	double time;
	int i,k;
	double x[] = new double [NI];
	double temp[] = new double [NI];

    	if(bufferg == null)
      		bufferg = buffer.getGraphics();
	/*	INITIAL VALUES FOR TEMPERATURE		*/
        for(i=0;i<NI;i++)
	{
		temp[i]=1.0;
	}
	time=0.0;
		
	for(i=0;i<100;i++)
	{
		Dimension d = getSize();
		// フォントを設定する
    		g.setFont(new Font("Serif", Font.BOLD, 20));
		g.drawString("Time= "+(float)(time), 50, 50);
		heat(i,NI,time,x,temp);
		bufferg.drawLine((int)(x[0]*200.0+150),(int)(0.9*d.height-1.0*200.0),                            (int)(x[NI-1]*200.0+150), (int)(0.9*d.height-200.0));

		// バッファを描画する
 		bufferg.setColor(Color.white);
		bufferg.fillRect(0, 0, d.width, d.height);
    		bufferg.setColor(Color.red);
		for(k=1;k<NI;k++)
		{
		   bufferg.drawLine((int)(x[k-1]*200.0+150), (int)(0.9*d.height-temp[k-1]*200.0),                            (int)(x[k]*200.0+150), (int)(0.9*d.height-temp[k]*200.0));
		}
		g.drawImage(buffer, 0, 0, this);
	}
   }

/* サブルーチン：熱伝導計算
*************************************************************
     Program for heat conduction with implicit method
*************************************************************/
static void heat(int i, int NI,double time, double x[], double temp[])
 {
        int NIM1;
	double dx,xmax,r;
	double dtime,tmax;
	double temp_n[]    = new double [NI];
	NIM1=NI-1;
	xmax=1.0;
	r=0.2;
	dx=xmax/(double)NIM1;
	dtime=r*dx*dx;
		
	x[0]=0.0;
	for(i=1;i<NI;i++)
        {
		x[i]=x[i-1]+dx;
	}
	time += dtime;
        temp[0]=0.0;
        temp[NIM1]=0.0;
        temp_n[0]=0.0;
        temp_n[NIM1]=0.0;

	/*	CALCILATE TEMPERUTURE WITH MARCHING    */
     	for(i=1;i<NIM1;i++)
	{
		各自記入せよ;
	}
	for(i=0;i<NI;i++)
	{
        	各自記入せよ;
	}
	
}
}