YW function

Yule-Walker solver