program fib args n clear qui set obs `n' qui gen a=1 qui replace a=a[_n-1]+a[_n-2] in 3/l end