clear all filename = 'signal'; r = 10000; %distance to source, in tens of kpc data = load(strcat(filename, '.dat')); t = 0.001*data(:,1); h = 8.8524e-24*(10/r)*data(:,2); clear data; plot(t,h); title('Gravitational Wave Burst'); xlabel('Time (s)'); ylabel('Strain');