* This is an illustrative example.  You may modify this code to adopt to your model.

cal 1967 1 1
all 2001:1

* Please use a proper folder name.
open data ukenergy.txt
data(format=prn, org=obs) / year Ew Xw GDPcurr GDP95 Tave Tjan POP

print / year ew xw GDPcurr GDP95 Tave Tjan POP

set GDPdef =GDPcurr / GDP95
set NPenergy = Xw/Ew
set RPenergy = Npenergy / GDPdef
set Lenergy = log(Ew)
set LenergyX = log(Xw)
set LGDP95 = log(GDP95)
set LRPenergy = log(RPenergy)
set TREND = t

smpl 1967:1 1997:1
print / year

* LEnergy = a + b1*LGDP95 + b2*LRPenergy + b3*Tave + b4*trend + .. + e
LIN LEnergyX
# constant LGDP95 LRPenergy Tave trend

LIN(define=eq1) LEnergyX
# constant LEnergyX{1} LGDP95 LRPenergy Tave trend

forecast(print) 1 4 1998:1
# eq1 y_for

print 1967:1 2001:1 LEnergyX y_for

smpl 1967:1 2001:1
graph(style=line) 2
# LEnergyX
# y_for