diff --git a/Virtual Sensing/README.md b/Virtual Sensing/README.md
index 14c1e5b432b6f64a42d15f6d81989c1aa42bb6ff..116957cf83a44495fdb2a2cb88dae8f2db479f32 100644
--- a/Virtual Sensing/README.md	
+++ b/Virtual Sensing/README.md	
@@ -39,7 +39,8 @@ Pm = ptSrcField(vPos, mPos, 7.5e2); % Primary-to-monitoring transfer function (f
 Pe = ptSrcField(vPos, ePos, 7.5e2); % Primary-to-virtual transfer function (f = 750 Hz)
 
 % Perform estimation and acquire the normalised square estimation error
-[~, ~, ~, ~, ~, ~, normSqrErr] = obsFilt(Pe, Pm, [], 1e2); % Use regularisation factor 100
+Oopt = obsFilt(Pe, Pm, [], 1e2); % Use regularisation factor 100 to calculate the observation filter
+[~, ~, ~, normSqrErr] = obsFileEst(Pm, Oopt, Pe, []); % Perform estimation
 ```
 
 #### Acknowledgments