Řešení soustavy rovnic pomocí Newtonovy metody

U verze MAPLE 5.1 výpočet soustavy nelineárních rovnic pomocí příkazu fsolve často zkolabuje, protože navržený "solver" není dostatečně stabilní. Toto se děje především v případech, kdy soustava rovnic obsahuje polynomiální rovnice (stačí jedna), nebo jsou v rovnicích pohromadě velká a malá čísla. Jedním z řešení, jak vypočítat soustavu rovnic bez použití příkazu fsolve je užití procedury NEWTON.

Proceduru Newton načteme (obvykle na počátku výpočtu) pomocí příkazu read.

obecně:

read `<cesta>newton.mpl`:

Příklad při umístění procedury na disketě v mechanice A:

Pro řešení soustavy rovnic pomocí Newtonovy metody pak zadáme následující příkaz:

obecně:

Newton([Rovnice1, Rovnice2, ..., RovniceN],[Neznámá1=odhad, Neznámá2=odhad, ..., NeznámáN=odhad])

pro každou neznámou musíme zapsat počáteční (rozumný) odhad výsledku.

Příklad pro řešení tří rovnic:

newton.mpl