Szukanie przybliżonego pierwiastka funkcji metodą siecznych w Octave / MatLab

Wpis gościnny hakerzy, znowu dzisiaj poruszający tematykę Octave i MatLaba. Dzisiaj szukamy przybliżonego pierwiastka funkcji metodą siecznych. Algorytm został napisany w Octavie (3.6.1), prawdopodobnie i tym razem działa w pełni w MatLabie.  Poniższy skrypt metodą siecznych znajduje przybliżenie z dokładnością 10^-4 w przedziale [-3, -2], pierwiastka funkcji f(x)=x^{3}+3x^{2}-1. Z tego względu warto spojrzeć na 5 linijkę kodu, tam została zdefiniowana powyższa funkcja. Program posiada kilka komentarzy w kodzie źródłowym, więc zrozumienie jego nie powinno być trudne. Czytaj więcej