Please help me with the listed revisions for the code below. Thank you. from math import pi, sin def m_trapz(f, a, b, n=100): res = (f(a)+f(b))/2 h = (b-a)/n for k in range(1, n): res = res+f(a+h*k) da = (f(a+h)-f(a))/h db = (f(b)-f(b-h))/h return res*h - (db-da)*(h**2)/12 def f(x): return sin(x) a = 0 b = pi/2 exact = 1 print('n\tI\tError') for n in [10, 100, 1000, 10000]: I = m_trapz(f, a, b, n) print('%d\t%.4f\t%.2e' % (n, I, I-exact))

C++ for Engineers and Scientists
4th Edition
ISBN:9781133187844
Author:Bronson, Gary J.
Publisher:Bronson, Gary J.
Chapter7: Arrays
Section: Chapter Questions
Problem 7PP: (Numerical) Using the srand() and rand() C++ library functions, fill an array of 1000 floating-point...
icon
Related questions
Question
100%
Please help me with the listed revisions for the code below. Thank you. from math import pi, sin def m_trapz(f, a, b, n=100): res = (f(a)+f(b))/2 h = (b-a)/n for k in range(1, n): res = res+f(a+h*k) da = (f(a+h)-f(a))/h db = (f(b)-f(b-h))/h return res*h - (db-da)*(h**2)/12 def f(x): return sin(x) a = 0 b = pi/2 exact = 1 print('n\tI\tError') for n in [10, 100, 1000, 10000]: I = m_trapz(f, a, b, n) print('%d\t%.4f\t%.2e' % (n, I, I-exact))
Expert Solution
steps

Step by step

Solved in 2 steps

Blurred answer
Knowledge Booster
Array
Learn more about
Need a deep-dive on the concept behind this application? Look no further. Learn more about this topic, computer-science and related others by exploring similar questions and additional content below.
Similar questions
  • SEE MORE QUESTIONS
Recommended textbooks for you
C++ for Engineers and Scientists
C++ for Engineers and Scientists
Computer Science
ISBN:
9781133187844
Author:
Bronson, Gary J.
Publisher:
Course Technology Ptr