14c14
<       real fact2,w1cm,betafix
---
>       real fact2,w1cm
18c18
<       common /one/ bestW1,bestW2,bestbeta,bestE0,minchi2,timbeta
---
>       common /one/ bestW1,bestW2,bestbeta,bestE0,minchi2
43c43
<        read (5,*)timbeta
---
>        read (5,*)beta
102c102
<        write (6,*)'Initial beta is ',timbeta
---
>        write (6,*)'Initial beta is ',beta
131c131
< C       betac = (betac / 100) * beta
---
>        betac = (betac / 100) * beta
135c135
< C       betac = 0
---
>        betac = beta + betac
144c144
<        do i=1,3
---
>        do i=1,5
147,148c147,148
< C          vertix(i,3) = beta
<           vertix(i,3) = E0
---
>           vertix(i,3) = beta
>           vertix(i,4) = E0
153,154c153,154
< C       vertix(4,3) = betac
<        vertix(4,3) = E0c
---
>        vertix(4,3) = betac
>        vertix(5,4) = E0c
156,157c156,157
<        do i=1,4
<           do k=1,3
---
>        do i=1,5
>           do k=1,4
167c167
<        call amoeba(vertix,theresult,4,3,3,ftol,phi2,iter)
---
>        call amoeba(vertix,theresult,5,4,4,ftol,phi2,iter)
169c169
<        do i=1,4
---
>        do i=1,5
172,173c172,173
<           vertix(i,3) = bestE0
< C          vertix(i,4) = 
---
>           vertix(i,3) = bestbeta
>           vertix(i,4) = bestE0
178,179c178,179
< C       vertix(4,3) = betac - (betac / 15)
<        vertix(4,3) = E0c - (E0c /500)
---
>        vertix(4,3) = betac - (betac / 15)
>        vertix(5,4) = E0c - (E0c /500)
184,186c184,186
<   
<        do i=1,4
<           do k=1,3
---
> 
>        do i=1,5
>           do k=1,4
253c253
< 	real phi,m1,m2,corner(4),calcenergy1(1000),calcprob(1000)
---
> 	real phi,m1,m2,corner(4),calcenergy1(1000),calcprob(1000),beta
257c257
<         common /one/ bestW1,bestW2,bestbeta,bestE0,minchi2,timbeta
---
>         common /one/ bestW1,bestW2,bestbeta,bestE0,minchi2
272,273c272
< C        timbeta=-1
< C        write(6,*)timbeta
---
> 
276c275
< 	  calcprob(iphi)=(1.0/(2*3.14159))*(1.0+timbeta*0.5*
---
> 	  calcprob(iphi)=(1.0/(2*3.14159))*(1.0+corner(3)*0.5*
282c281
<      &     (corner(3)/m2+(2.0*(sqrt(E0*
---
>      &     (corner(4)/m2+(2.0*(sqrt(E0*
