#!/bin/bash
math <<EOF
gamma=0.05;
epsilon=0.3;
Delta[w_]:=gamma (I+1/Pi Log[(1-w)/(1+w)]);
A[w_]:=-1/Pi Im[1/(w-epsilon+Delta[w])];
tab=Table[{w,A[w]},{w,-1+0.001/2,1,0.001}];
Export["reference.dat",tab];
EOF
