J
jadaayu
Guest
Kompletní eliptický integrál prvního druhu ve tvaru Legendrovy K (K2), 0 <k2 <1, / 2 0 2 2 2 1 sin () () dt K kt k nelze hodnotit z hlediska elementárních funkcí. Následující algoritmus, díky CF Gauss, generuje posloupnost aritmetických {} a posloupnost geometrické prostředky {} mld. Kč, kde 43 a0 = 1, b0 = 1 k 2 = (-1 + BN- 1) / 2, BN = 1 1 nnabn = 1, 2, .... Je známo, že obě sekvence mají společný limit g a že mld. Kč, pro všechny n. Navíc, K (K2) = 2g zápisu MATLAB funkce K = compK (K2), která implementuje tento algoritmus. Vstupní parametr K2 kandiduje na K2. Pomocí smyčky a zároveň vytvářet po sobě jdoucí členy obou posloupností, ale nelze uložit všechna čísla získané v průběhu výpočtů. Pokračovat v provádění smyčky while tak dlouho, dokud - BN eps, kde eps je stroj epsilon