sp;size(model.Cov,2),&nbsp;
<br>
&nbsp;&nbsp;model.Cov&nbsp;=&nbsp;reshape(model.Cov,1,1,2);&nbsp;
<br>
<span class=keyword>end</span>
<br>

<br>
<span class=comment>%&nbsp;finds&nbsp;out&nbsp;decision&nbsp;function&nbsp;which&nbsp;is&nbsp;generaly&nbsp;quadratic
</span><br>
quad_model=bayesdf(model);
<br>
a&nbsp;=&nbsp;quad_model.A;
<br>
b&nbsp;=&nbsp;quad_model.B;
<br>
c&nbsp;=&nbsp;quad_model.C;
<br>

<br>
<span class=comment>%&nbsp;get&nbsp;parameters
</span><br>
p1=model.Prior(1);&nbsp;p2=model.Prior(2