#beta known and alpha unknown   DST
rm(list = ls())
ltmode1=c(300,400,500)
utmode1=c(400,500,650)
pa=array(3)
ba=array(3)
lR=array(3)
uR=array(3)
mf=array(3)
mw=array(3)
mwf=array(3)
lalphaaggY=c()
ualphaaggY=c()
uTmin=c()
lTmin=c()
lmineta=c()
umineta=c()
#T=1950
B=2
c=0.3
#T=c(50,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,300,310,320,345,346,347,348,349,350,360,370,380,390,400,410,420,421,422,423,424,425,426,427,428,429,430,431,432,433,434,435,436,437,438,439,440,450,460,461,462,463,464,465,466,470,472,475,478,480,485,488,490,495,496,497,
498,499,500,510,520,530,540,550,560,570,580,590,591,592,593,594,595,596,597,598,599,600,610,620,630,640,650,660,670,680,690,691,692,693,694,695,696,697,698,699,700,710,720,730,731,800,810,820,830,840,841,842,843,844,845,846,847,848,849,900,1000,1100,1120,1130,1140,1150,1200,1250,1300,1350,1400,1600,1700,1900,2000)
for(jj in 1:length(T)){

for(j in 1:3){
f=function(a) {((B-1)/B)^(1/B)*a-ltmode1[j]}
pa[j]=uniroot(f,c(1,10000))$root

f1=function(a) {((B-1)/B)^(1/B)*a-utmode1[j]}
ba[j]=uniroot(f1,c(1,10000))$root[1]

lR[j]=exp(-(T[jj]/pa[j])^2)
uR[j]=exp(-(T[jj]/ba[j])^2)

mw[j]=lR[j]
mf[j]=1-uR[j]
mwf[j]=1-mw[j]-mf[j]

}

myw=prod(mw)+(mw[1]*mwf[2]*mw[3])+(mw[2]*mwf[1]*mw[3])+(mw[1]*mwf[3]*mw[2])+(mw[1]*mwf[2]*mwf[3])+(mwf[1]*mw[2]*mwf[3])+(mwf[1]*mwf[2]*mw[3])
myf=prod(mf)+(mf[1]*mwf[2]*mf[3])+(mf[2]*mwf[1]*mf[3])+(mf[1]*mwf[3]*mf[2])+(mf[1]*mwf[2]*mwf[3])+(mwf[1]*mf[2]*mwf[3])+(mwf[1]*mwf[2]*mf[3])
k=((mf[1]*mf[2]*mw[3])+(mf[1]*mw[2]*mf[3])+(mw[1]*mf[2]*mf[3])+(mw[1]*mw[2]*mf[3])+(mw[1]*mf[2]*mw[3])+(mf[1]*mw[2]*mw[3])+(mf[1]*mwf[2]*mw[3])
+(mf[1]*mw[2]*mwf[3])+(mw[1]*mwf[2]*mf[3])+(mw[1]*mf[2]*mwf[3])+(mwf[1]*mw[2]*mf[3])+(mwf[1]*mf[2]*mw[3]))
mywf=(mwf[1]*mwf[2]*mwf[3])+k


mdw=myw/(1-k)
mdw
mdf=myf/(1-k)
mdf
mdwf=(mwf[1]*mwf[2]*mwf[3])/(1-k)
mdwf
mdw+mdf+mdwf
lARY=mdw
lARY
uARY=1-mdf
uARY
lalphaaggY[jj]=T[jj]/(sqrt(-log(lARY)))
ualphaaggY[jj]=T[jj]/(sqrt(-log(uARY)))

uTmin[jj]=ualphaaggY[jj]*(c/(B-1))^(1/B)
lTmin[jj]=lalphaaggY[jj]*(c/(B-1))^(1/B)
lmineta[jj]=10^5*( uTmin[jj]^(B-1)/(ualphaaggY[jj])^(B)+c/uTmin[jj])
umineta[jj]=10^5*( lTmin[jj]^(B-1)/(lalphaaggY[jj])^(B)+c/lTmin[jj])
}
lTmin
uTmin
uTmin-lTmin
lmineta
umineta
umineta-lmineta
sort(umineta-lmineta)[1]
orgeta=which.min(umineta-lmineta)
lmineta[orgeta]
umineta[orgeta]
lTmin[orgeta]
uTmin[orgeta]
orgeta1=which.min(uTmin-lTmin)
lmineta[orgeta1]
umineta[orgeta1]
lTmin[orgeta1]
uTmin[orgeta1]

<pre>An Error occurred while handling another error:
yii\web\HeadersAlreadySentException: Headers already sent in  on line 0. in /var/www/html/prof-homepages/vendor/yiisoft/yii2/web/Response.php:366
Stack trace:
#0 /var/www/html/prof-homepages/vendor/yiisoft/yii2/web/Response.php(339): yii\web\Response-&gt;sendHeaders()
#1 /var/www/html/prof-homepages/vendor/yiisoft/yii2/web/ErrorHandler.php(136): yii\web\Response-&gt;send()
#2 /var/www/html/prof-homepages/vendor/yiisoft/yii2/base/ErrorHandler.php(135): yii\web\ErrorHandler-&gt;renderException()
#3 [internal function]: yii\base\ErrorHandler-&gt;handleException()
#4 {main}
Previous exception:
yii\web\HeadersAlreadySentException: Headers already sent in  on line 0. in /var/www/html/prof-homepages/vendor/yiisoft/yii2/web/Response.php:366
Stack trace:
#0 /var/www/html/prof-homepages/vendor/yiisoft/yii2/web/Response.php(339): yii\web\Response-&gt;sendHeaders()
#1 /var/www/html/prof-homepages/vendor/yiisoft/yii2/base/Application.php(656): yii\web\Response-&gt;send()
#2 /var/www/html/prof-homepages/vendor/faravaghi/yii2-filemanager/models/Files.php(696): yii\base\Application-&gt;end()
#3 /var/www/html/prof-homepages/vendor/faravaghi/yii2-filemanager/controllers/FilesController.php(484): faravaghi\filemanager\models\Files-&gt;getFile()
#4 [internal function]: faravaghi\filemanager\controllers\FilesController-&gt;actionGetFile()
#5 /var/www/html/prof-homepages/vendor/yiisoft/yii2/base/InlineAction.php(57): call_user_func_array()
#6 /var/www/html/prof-homepages/vendor/yiisoft/yii2/base/Controller.php(180): yii\base\InlineAction-&gt;runWithParams()
#7 /var/www/html/prof-homepages/vendor/yiisoft/yii2/base/Module.php(528): yii\base\Controller-&gt;runAction()
#8 /var/www/html/prof-homepages/vendor/yiisoft/yii2/web/Application.php(103): yii\base\Module-&gt;runAction()
#9 /var/www/html/prof-homepages/vendor/yiisoft/yii2/base/Application.php(386): yii\web\Application-&gt;handleRequest()
#10 /var/www/html/prof-homepages/frontend/web/index.php(18): yii\base\Application-&gt;run()
#11 {main}</pre>