给定函数式和x值，怎么求一系列的y值的stata代码

如题，给定一个函数式Yₓ=31.62+0.83Yₓ₋₁-0.22Yₓ₋₂-0.05Yₓ₋₃-0.13Yₓ₋₄+εₓ-εₓ₋₁和x值，怎么求x=35到72所有的的y值的stata代码

clear all
set more off

// 定义初始值
local Y1 = 0
local Y2 = 0
local Y3 = 0
local Y4 = 0

// 定义函数
gen Y = .

forvalues x = 35/72 {
    replace Y = 31.62 + 0.83*Y[_n-1] - 0.22*Y[_n-2] - 0.05*Y[_n-3] - 0.13*Y[_n-4] + rnormal(0,1)
    if _n == 1 replace Y = Y1
    if _n == 2 replace Y = Y2
    if _n == 3 replace Y = Y3
    if _n == 4 replace Y = Y4

    // 更新前四个值
    if _n == 5 {
        local Y1 = Y[1]
        local Y2 = Y[2]
        local Y3 = Y[3]
        local Y4 = Y[4]
    }
}

// 显示结果
list Y in 35/72