#!/usr/local/bin/tcc -run
#include <stdio.h>
#include <math.h>
#define NMAX 5

main(){
	FILE *fp;
	int i, kensu;
	double aveX, aveY, aveXX, aveYY, aveXY, vX, vY, vXY, v_XY, r, inpdata[NMAX][2];
	
	printf("データの件数:");
	fscanf(fp,"%d",&kensu);
	if(kensu>NMAX){
		printf("100以下の件数を入力してください")
	}
	for (i=0;i<kensu;i++){
		printf("i+1番目のデータ:");
		fscanf(fp,"%lf %lf",&inpdata[i][0],&inpdata[i][1]);
	}
    
    aveX=0;
    for(i=0;i<NMAX-1;i++){
    	aveX=aveX+inpdata[i][0];
    	aveY=aveY+inpdata[i][1];
    	aveXX=aveXX+inpdata[i][0]*inpdata[i][0];
    	aveYY=aveYY+inpdata[i][1]*inpdata[i][1];
    	aveXY=aveXY+inpdata[i][0]*inpdata[i][1];
    }
    
    aveX=aveX/kensu;
    aveY=aveY/kensu;
    aveXX=aveXX/kensu;
    aveYY=aveYY/kensu;
    aveXY=aveXY/kensu;
    vX=aveXX-aveX*aveX;
    vY=aveYY-aveY*aveY;
    vXY=aveXY-aveX*aveY;
    v_XY=sqrt(vX*vY);
    r=vXY/v_XY;
    
    printf("相関係数=%lf",r);
	
}