A new algorithm is derived for the standard fixed-interval linear smoothing problem in which the signal is generated by a state model. The structure of this new algorithm allows the smoothed estimate to be easily updated in response to a change in the initial state covariance matrix 

 , since, unlike in existing algorifiuns, the relevant Riccati equation is entirely independent of 

 . The derivation of the algorithm is based on properties of complementary models.