Рис.18. Реализация метода Гаусса.

Метод Тейлора

Файл teilor.mac:

f(x) := -(x-3)^2 + 5;

teilor(a, b, e) := block( [x0:0, result:0, k:1], if a >= b then block(

print("Incorrect interval!!!"), break

),

if e <= 0.000001 then block( print("Incorrect presicion!"), break

),

x0 : (a+b)/2.0, result : f(x0)*(b-a),

while (abs(at(diff(f(x), x, k),[x = x0])/(k!)) > e) do(

result : result + at(diff(f(x),x,k),[x=x0])/(k+1)!*((b-x0)^(k+1) - (a-x0)^(k+1)), k : k + 1

),

return (result)

);

print("<Integral of function: -(x-3)^2+5>"); a:read("Enter lower limit: "); b:read("Enter upper limit: ");

e: read("Enter precision: "); if e <= 0 then block(

print("Incorrect e!"), break

);

print("<Method of teilor>::: ", teilor(a, b, e));

Рис. 19. Метод разложения подынтегральной функции в ряд Тейлора.