!m=50
!n=20
wfcreate mcarlo_results_1 u 1 !n
!theta=5
!gamma=2
rndseed 1213859
vector(!m) weibull
vector(!m) psi_m
vector(!m) qd
smpl 1 !n
for !i = 1 to !m
series y=@rweib(!theta,!gamma)

logl qnmb
qnmb.append @logl log_wei
qnmb.append @param c(1)
qnmb.append log_wei=ln((!gamma*y^(!gamma-1))/(c(1)^!gamma)))-(y/c(1))^!gamma
weibull(!i)=c(1)
psi_m(!i)=c(1)^!gamma
next