📝 **Task:** Implement `fit_linear(xs, ys, lr=0.01, epochs=200)` returning `(w, b)` rounded to 3 decimals. 📋 Implement the function above. Tests run automatically. 💡 **Hint:** Re-read the theory if you get stuck.