c=Sqrt[(1/2)^2+1/2*((1)^2-(1/2)^2)*(1+(2/Pi)*ArcTan[(x+(8*Tan[(Pi/((1)^2-(1/2)^2))*((3/4)^2-1/2*((1/2)^2+1^2))]))/(8)])];

v0=3/4;

x1=10+1/2;
x3=x1;
x2=-x1;
x4=x2;
x1Start=-800-3/5;
x1End=800+2/5;
x2Start=-800-2/5;
x2End=800+3/5;

xstar[y_]:=NIntegrate[c/(c^2-v0^2), {x,  x3,y}, WorkingPrecision->30];

T=0;

tout[y_]:=T-NIntegrate[v0/(c^2-v0^2), {x,  x1,y}, WorkingPrecision->30];

XLargeStep=Abs[( x1Start-x1End)/(1001)];
XLargeStepNumber1=Floor[Abs[(x1Start-x1)]/XLargeStep];
XLargeStepNumber2=Floor[Abs[(x2Start -x1)]/XLargeStep];
XLargeStart1=( x1Start+XLargeStep*(XLargeStepNumber1+1));
XLargeStart2=( x2Start+XLargeStep*(XLargeStepNumber2+1));

Plot[{xstar[x], -tout[x]}, {x, 1, 400}, PlotLegends->{xstar, t}]

datsmall=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 100-3/5, 800+2/5,100 }];
C:\Users\dudlra13\Dropbox\BH\BEC\Massive\DensityDensity\InputData\XandT

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4Small.txt",datsmall,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4Small.txt"

datmed=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 20-3/5, 800+2/5,10 }];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4Med.txt",datmed,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4Med.txt"

datlarge=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 12-3/5, 800+2/5,4 }];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4Large.txt",datlarge,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4Large.txt"

datxlarge=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, XLargeStart1, 800+2/5,XLargeStep}];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4XLarge.txt",datxlarge,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4XLarge.txt"

dat2single=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 602+1/3, 603+1/3,100 }];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeSingle.txt",dat2single,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeSingle.txt"

dat2small=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 100-2/5, 800+3/5,100 }];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeSmall.txt",dat2small,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeSmall.txt"

dat2med=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 20-2/5, 800+3/5,10 }];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeMed.txt",dat2med,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeMed.txt"

dat2large=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, 12-2/5, 800+3/5,4 }];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeLarge.txt",dat2large,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeLarge.txt"

dat2xlarge=Table[{N[xstar[x], 30],N[tout[x], 30], N[x, 30]},{x, XLargeStart2, 800+3/5,XLargeStep}];

Export["C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeXLarge.txt",dat2xlarge,"Table","FieldSeparators"->" "]

"C://Users//dudlra13//Dropbox//BH//BEC//Massive//DensityDensity//InputData//XandT//X4PrimeXLarge.txt"