// *  <applet code="toukousen3" width=500 height=500></applet>

import java.applet.*;
import java.awt.*;
import java.awt.event.*;
import Mitsui.*;   //Mitsui_Itoの入っているパッケージMitsuiを呼び出す。

public class toukousen3 extends Applet {
    Mitsui_and_Ito m = new Mitsui_and_Ito();

    public void init(){
        Graphics g=getGraphics();
        m.setGraphics(g);
        m.setScreenSize(getSize().width,getSize().height);
        m.setColor(Color.red);
        //[a,b]x[c,d]
        m.setArea(-1.2,1.2,-1.2,1.2);
        
    }
    
    public void paint(Graphics g) {
        //0, 1, -1, x, yを表示
        g.setColor(Color.red);
        int nx=100,ny=100;
        double xmin=-1.0,xmax=1.0,ymin=-1.0,ymax=1.0;
        double x,y;
        double [][]u=new double[nx+1][ny+1];
        
        double dx = (xmax-xmin)/nx;
        double dy = (ymax-ymin)/ny;
        for (int i=0;i<=nx;i++ ){
            x = xmin + i * dx;
            for (int j=0;j<=ny;j++ ) {
                y = ymin + j * dy;
                u[i][j] = x*x-y*y;
            }
        }
        // 等高線描画
        for (double h=-1.0;h<=1.0;h+=0.1) {
            if (h < 0.0)
                m.setColor(Color.blue);
            else
                m.setColor(Color.red);
            m.contln(xmin, xmax, ymin, ymax, u, nx, ny, h); 
/**     
        m.contln(xmin, xmax, ymin, ymax, u, nx, ny, h);は 
　　   Mitsui_and_Itoから下のcontlnを呼び出している。([xmin,xmax]×[ymin,ymax])

         public void contln(double xmin, double xmax, double ymin, double ymax,
          double [][]u, int nx, int ny, double h) {
          int i, j;
            double dx = (xmax-xmin)/nx, dy=(ymax-ymin)/ny;
            double x[] = new double[nx+1];
            double y[] = new double[ny+1];
            for (i=0;i<=nx;i++)
            x[i] = xmin + i * dx;
            for (j=0;j<=ny;j++)
            y[j] = ymin + j * dy;
            for (i=0;i<nx;i++){
            for (j=0;j<ny;j++){
            celldraw(x[i],x[i+1],x[i+1],
                             y[j],y[j],y[j+1],
                             u[i][j],u[i+1][j],u[i+1][j+1],h);
                             celldraw(x[i],x[i],x[i+1],
                             y[j],y[j+1],y[j+1],
                             u[i][j],u[i][j+1],u[i+1][j+1],h);
                             }
                             }
        }
        // 等高線描画ルーチン
        public void celldraw(double x1, double x2, double x3,
        double y1, double y2, double y3,
                             double u1, double u2, double u3,
                             double h){
            
            // 以下、オリジナルでは new double[3] だったが、
            // 明らかなバグであると思われる。
            double px[] = new double[4];
            double py[] = new double[4];
            int n;
            double r;
            
            n = 0;
            
            if ((u1-h)*(u2-h)<=0){
                n++;
                if (u1==h){
                    px[n] = x1;
                    py[n] = y1;
                } else if(u2==h) {
                    px[n] = x2;
                    py[n] = y2;
                } else {
                    r = (u1 - h) / (u1 - u2);
                    px[n] = (1 - r) * x1 + r * x2;
                    py[n] = (1 - r) * y1 + r * y2;
                }
            }
            if ((u2-h)*(u3-h)<=0){
                n++;
                if (u2==h){
                    px[n] = x2;
                    py[n] = y2;
                } else if(u3==h) {
                    px[n] = x3;
                    py[n] = y3;
                } else {
                    r = (u2 - h) / (u2 - u3);
                    px[n] = (1 - r) * x2 + r * x3;
                    py[n] = (1 - r) * y2 + r * y3;
                }
            }
            if ((u3-h)*(u1-h)<=0){
                n++;
                if (u3==h){
                    px[n] = x3;
                    py[n] = y3;
                } else if(u1==h) {
                    px[n] = x1;
                    py[n] = y1;
                } else {
                    r = (u3 - h) / (u3 - u1);
                    px[n] = (1 - r) * x3 + r * x1;
                    py[n] = (1 - r) * y3 + r * y1;
                }
            }
            
            if (n==2){
                move(px[1], py[1]);
                draw(px[2], py[2]);
            } else if (n==3) {
                move(px[1], py[1]);
                draw(px[2], py[2]);
                move(px[2], py[2]);
                draw(px[3], py[3]);
                move(px[3], py[3]);
                draw(px[1], py[1]);
                }
            }
*/
                
        }
    }
}