import math import numpy print("x: f*: f:") x = 0.2 for i in range(11): y = 3 * x + 0.6 z = 1 - x**2 g = 0.9 * x + 1 s1 = 0 s2 = 0 s3 = 0 n = 0 while y**(2*n+1)/math.factorial(2*n+1) > 10**(-6)/2: n += 1 for j in range(n + 1): s1 += (-1)**j*y**(2*j+1)/math.factorial(2*j+1) n = 0 w = 1.0 w2 = (1 + (g))/2.0 while abs(w2-w)*g > 10**(-6)/1.92: w = w2 w2 = (w+ g/w) / 2.0 s2 = w2/z while s2**-(2*n+1)/(2*n+1) > 10**(-6)/2: n += 1 for j in range(n + 1): s3 += (-1)**n*s2**(-2*n-1)/(2*n+1) s_3 = (math.pi/2)*numpy.sign(s2)-s3 print(round(x,2), " ", s1+s_3 ," ", math.atan((0.9*x+1)**0.5/(1-x**2))+math.sin(3*x+0.6)) x += 0.01