#!/bin/sh
math <<EOF

lu = Import["Delta-u.dat"];
ld = Import["Delta-d.dat"];
fu = Interpolation[lu];
fd = Interpolation[ld];
\[Omega]s = Table[1/1.01^n, {n, 0, 1500}];
eps = 10^-8;
\[Omega]s = Join[\[Omega]s, {eps}];
\[Omega]s = N[Sort[Join[\[Omega]s, -\[Omega]s]]];

r = RotationMatrix[Pi/4];

l = Map[{#, Transpose[r] . {{fu[#], 0}, {0, fd[#]}} . r} &, \[Omega]s];
Flatten2[{x_, {{a_, b_}, {c_, d_}}}] := {x, Re[a], Re[b], Re[c], Re[d]};
lfl = Map[Flatten2, l];
Export["Delta.dat", lfl];

EOF