## Revision 1486

25 25
```    \n\n The implementation of this algorithm was inspired by
```
26 26
```    "Numerical Recipes in C", 2nd edition,
```
27 27
```    Press, Teukolsky, Vetterling, Flannery, chapter 2
```
28

29
```    Edited by Steve Martell on Dec 30, 2013 to address a bug
```
30
```    where the lower bound of the square dvar_matrix was greater than 1.
```
28 31
```*/
```
29 32
```dvariable det(const dvar_matrix& aa)
```
30 33
```{
```
......
129 132
```    }
```
130 133
```  }
```
131 134
```  double det=d;
```
132
```  part_prod(1)=d*bb(1,1);
```
135

136
```  // SM Bug 129, issue appears to be at this line
```
137
```  // part_prod is declared above as dvector(lb,ub)
```
138
```  // cout<<"Bug 129 ="<<part_prod(lb)<<endl;
```
139
```  // part_prod(1)=d*bb(1,1);  // replaced this line with:
```
140
```  part_prod(lb) = d*bb(lb,lb);
```
141
```  // cout<<"Ok got this far; det = "<<det<<endl;
```
133 142
```  for (j=lb+1;j<=ub;j++)
```
134 143
```  {
```
135 144
```    part_prod(j)=part_prod(j-1)*bb(j,j);
```
trunk/contrib/statslib/statsLib.h (revision 1486)
92 92
```dvar_vector plogis( const dvector& x, const prevariable& location, const prevariable& scale );
```
93 93
```dvar_vector plogis( const dvar_vector& x, const prevariable& location, const prevariable& scale );
```
94 94
```// template <typename T1, typename T2>
```
95
```// T1 plogis( const T1& x, const T2& location, const T2& scale );
```
95
```// T1 plogis( const T1& x, const T2& locatiogn, const T2& scale );
```
96 96

97 97
```// Exponential logistic distribution
```
98 98
```dvar_vector eplogis(const dvar_vector& x, const dvariable& alpha, const dvariable& beta, const dvariable& gamma);
```

Also available in: Unified diff