set deltah = 0.43010E-04 set n=1 data trans.dat read kh 1 read trans 2 set del2 = deltah**2 * (kh*2997.9)**(3+n)*trans**2 set lkh = lg(kh) set ldel2 = lg(del2) limits -1.5 -0.5 -1 1 connect lkh ldel2 ticksize -1 10 -1 10 box rel -3 0 draw 3 0 rel -0.86646 -3 draw -0.86646 3