#include <stdio.h>
#include <math.h>

int main(int argc, char **argv){
int Narray;
int k;
Narray = 2250;
double zeta[Narray];
double dzeta = 0.225/Narray;
double tau[Narray];
zeta[0]=0.9;
tau[0]=0;

for (k=1; k<Narray-1 ; k++){
    zeta[k]=zeta[k-1]+dzeta;
    tau[k]=tau[k-1]+dzeta*zeta[k]/sqrt(0.987654*(zeta[k]-0.9)*(9/8-zeta[k]));
    printf("zeta = %lf tau = %lf n", zeta[k], tau[k]);
    }
return 0;
}