diff --git a/VBA_bDCM_lesion.m b/VBA_bDCM_lesion.m
index 03bf0651..905c650a 100755
--- a/VBA_bDCM_lesion.m
+++ b/VBA_bDCM_lesion.m
@@ -1,4 +1,4 @@
-function results=VBA_lesionAnalysis(posterior,out)
+function results=VBA_bDCM_lesion(posterior,out)
 
 %% Get normal profile
 
@@ -22,7 +22,9 @@
         posterior.muTheta,...
         posterior.muPhi,...
         out.u,...
-        Inf,Inf,out.options,...
+        Inf, ...
+        inf(size(posterior.a_sigma)), ...
+        out.options,...
         posterior.muX0);
    
         rnormal.y=yp-er;
@@ -41,7 +43,9 @@
         posterior_l.muTheta,...
         posterior.muPhi,...
         out.u,...
-        Inf,Inf,out.options,...
+        Inf, ...
+        inf(size(posterior.a_sigma)), ...
+        out.options,...
         posterior.muX0);
    
         rlesion(iRegion).y=yp-er;