####DST3 be ravesh dr asadi
rm(list = ls())
ltmode1=c(300,400,500)
utmode1=c(400,500,650)
lE=c(700,750,800)
uE=c(750,800,850)
pa=array(3)
ba=array(3)
pb=array(3)
bb=array(3)
a1=array(3)
a2=array(3)
a3=array(3)
a4=array(3)
b1=array(3)
b2=array(3)
b3=array(3)
b4=array(3)
lR=array(3)
uR=array(3)
mf=array(3)
mw=array(3)
mwf=array(3)
lRt=array(3)
uRt=array(3)
mf1=array(3)
mw1=array(3)
mwf1=array(3)
uTmin=c()
lTmin=c()
lmineta=c()
umineta=c()
lbetaaggY1=c()
ubetaaggY1=c()
ualphaaggY=c()
lalphaaggY=c()
lbetaaggY=c()
ubetaaggY=c()
lbeta=c()
c=0.3
#T=150
#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(b) {(uE[j]/ltmode1[j])^b*((b-1)/b)+log(0.01)}
b1[j]=uniroot(f,c(1,1000))$root
a1[j]=ltmode1[j]*(b1[j]/(b1[j]-1))^(1/b1[j])
f1=function(b) {(lE[j]/utmode1[j])^b*((b-1)/b)+log(0.01)}
b2[j]=uniroot(f1,c(1,1000))$root[1]
a2[j]=utmode1[j]*(b2[j]/(b2[j]-1))^(1/b2[j])
f=function(b) {(uE[j]/utmode1[j])^b*((b-1)/b)+log(0.01)}
b3[j]=uniroot(f,c(1,1000))$root
a3[j]=ltmode1[j]*(b3[j]/(b3[j]-1))^(1/b3[j])
f1=function(b) {(lE[j]/ltmode1[j])^b*((b-1)/b)+log(0.01)}
b4[j]=uniroot(f1,c(1,1000))$root[1]
a4[j]=utmode1[j]*(b4[j]/(b4[j]-1))^(1/b4[j])
pa[j]=min(a1[j],a2[j],a3[j],a4[j])
ba[j]=max(a1[j],a2[j],a3[j],a4[j])
pb[j]=min(b1[j],b2[j],b3[j],b4[j])
bb[j]=max(b1[j],b2[j],b3[j],b4[j])
t=70
if (T[jj] > pa[j]){
lR[j]=exp(-(T[jj]/pa[j])^(bb[j]))}
else {
lR[j]=exp(-(T[jj]/pa[j])^(pb[j]))}
if (T[jj] > ba[j]){
uR[j]=exp(-(T[jj]/ba[j])^(pb[j]))}
else {
uR[j]=exp(-(T[jj]/ba[j])^(bb[j]))}


if (T[jj] > pa[j]){
lRt[j]=exp(-((t+T[jj])/pa[j])^bb[j])/exp(-(t/pa[j])^bb[j])}
else {
lRt[j]=exp(-((t+T[jj])/pa[j])^pb[j])/exp(-(t/pa[j])^pb[j])}
if (T[jj] > ba[j]){
uRt[j]=exp(-((t+T[jj])/ba[j])^pb[j])/exp(-(t/ba[j])^pb[j])}
else {
uRt[j]=exp(-((t+T[jj])/ba[j])^bb[j])/exp(-(t/ba[j])^bb[j])}


mw[j]=lR[j]
mf[j]=1-uR[j]
mwf[j]=1-mw[j]-mf[j]
mw+mf+mwf
mw1[j]=lRt[j]
mf1[j]=1-uRt[j]
mwf1[j]=1-mw1[j]-mf1[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


myw1=prod(mw1)+(mw1[1]*mwf1[2]*mw1[3])+(mw1[2]*mwf1[1]*mw1[3])+(mw1[1]*mwf1[3]*mw1[2])+(mw1[1]*mwf1[2]*mwf1[3])+(mwf1[1]*mw1[2]*mwf1[3])+(mwf1[1]*mwf1[2]*mw1[3])
myf1=prod(mf1)+(mf1[1]*mwf1[2]*mf1[3])+(mf1[2]*mwf1[1]*mf1[3])+(mf1[1]*mwf1[3]*mf1[2])+(mf1[1]*mwf1[2]*mwf1[3])+(mwf1[1]*mf1[2]*mwf1[3])+(mwf1[1]*mwf1[2]*mf1[3])
k=((mf1[1]*mf1[2]*mw1[3])+(mf1[1]*mw1[2]*mf1[3])+(mw1[1]*mf1[2]*mf1[3])+(mw1[1]*mw1[2]*mf1[3])+(mw1[1]*mf1[2]*mw1[3])+(mf1[1]*mw1[2]*mw1[3])+(mf1[1]*mwf1[2]*mw1[3])
+(mf1[1]*mw1[2]*mwf1[3])+(mw1[1]*mwf1[2]*mf1[3])+(mw1[1]*mf1[2]*mwf1[3])+(mwf1[1]*mw1[2]*mf1[3])+(mwf1[1]*mf1[2]*mw1[3]))
mywf1=(mwf1[1]*mwf1[2]*mwf1[3])+k



mdw1=myw1/(1-k)
mdw1
mdf1=myf1/(1-k)
mdf1
mdwf1=(mwf1[1]*mwf1[2]*mwf1[3])/(1-k)
mdwf1


lARtY1=mdw1
lARtY1
uARtY1=1-mdf1
uARtY1

ff=function(b){(((t+T[jj])^(b)-t^b)*(-log(uARY)))/(T[jj]^b)+log(uARtY1)}
ubetaaggY[jj]=uniroot(ff,c(2,100))$root
ff1=function(b){(((t+T[jj])^(b)-t^b)*(-log(lARY)))/(T[jj]^b)+log(lARtY1)}
lbetaaggY[jj]=uniroot(ff1,c(2,100))$root[1]

if (-log(lARY) < 1 ){
lalphaaggY[jj]=T[jj]/(-log(lARY))^(1/ubetaaggY[jj])}
else {
lalphaaggY[jj]=T[jj]/(-log(lARY))^(1/lbetaaggY[jj])}

if (-log(uARY) < 1 ){
ualphaaggY[jj]=T[jj]/(-log(uARY))^(1/lbetaaggY[jj])}
else {
ualphaaggY[jj]=T[jj]/(-log(uARY))^(1/ubetaaggY[jj])}



lTmin[jj]=lalphaaggY[jj]*(c/(lbetaaggY[jj]-1))^(1/lbetaaggY[jj])
uTmin[jj]=ualphaaggY[jj]*(c/(ubetaaggY[jj]-1))^(1/ubetaaggY[jj])
lmineta[jj]=10^5*( uTmin[jj]^(ubetaaggY[jj]-1)/(ualphaaggY[jj])^(ubetaaggY[jj])+c/uTmin[jj])
umineta[jj]=10^5*( lTmin[jj]^(lbetaaggY[jj]-1)/(lalphaaggY[jj])^(lbetaaggY[jj])+c/lTmin[jj])

}


lbetaaggY
ubetaaggY
ubetaaggY-lbetaaggY
lalphaaggY
ualphaaggY
ualphaaggY-lalphaaggY
lTmin
uTmin
uTmin-lTmin
lmineta
umineta
umineta-lmineta<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>