diff --git a/source/functions/simulink/model/irregExcF.m b/source/functions/simulink/model/irregExcF.m
index d8c41bc28..90b60fd4c 100644
--- a/source/functions/simulink/model/irregExcF.m
+++ b/source/functions/simulink/model/irregExcF.m
@@ -8,7 +8,7 @@
     B1= sin(bsxfun(@plus,A1,phaseRand(:,ii)));
     B11 = sin(bsxfun(@plus,w*time,phaseRand(:,ii)));
     C0 = bsxfun(@times,A*spread(ii),dw);
-    C1 = sqrt(bsxfun(@times,A*spread(ii),dw));
+    C1 = sqrt(bsxfun(@times,A*spread(ii).^2,dw));
     D0 =bsxfun(@times,squeeze(fExtMD(ii,:,:)),C0);
     D1 =bsxfun(@times,squeeze(fExtRE(ii,:,:)),C1);
     D11 = bsxfun(@times,squeeze(fExtIM(ii,:,:)),C1);
diff --git a/source/functions/simulink/model/irregnLYaw.m b/source/functions/simulink/model/irregnLYaw.m
index 2f86276c2..42073c399 100644
--- a/source/functions/simulink/model/irregnLYaw.m
+++ b/source/functions/simulink/model/irregnLYaw.m
@@ -75,7 +75,7 @@
     B1= sin(bsxfun(@plus,A1,phaseRandint));
     B11 = sin(bsxfun(@plus,w*time,phaseRandint));
     C0 = bsxfun(@times,A*WaveSpreadint,dw);
-    C1 = sqrt(bsxfun(@times,A*WaveSpreadint,dw));
+    C1 = sqrt(bsxfun(@times,A*WaveSpreadint.^2,dw));
     D0 =bsxfun(@times,fExtMDint,C0);
     D1 =bsxfun(@times,fExtREint,C1);
     D11 = bsxfun(@times,fExtIMint,C1);