Math

k = 1
pi = 0
 
for i in range(1000000): 
    if i % 2 == 0: 
        pi += 4/k 
    else: 
        pi -= 4/k 
    k += 2
 
print(pi)
 

# integration numerique par la methode des rectangles avec alpha = a
 
import numpy as np
import matplotlib.pyplot as plt
 
xmin = 0
xmax = 3*np.pi/2
nbx = 20
nbi = nbx - 1 # nombre d'intervalles
 
x = np.linspace(xmin, xmax, nbx)
y = np.cos(x)
plt.plot(x,y,"bo-")
 
integrale = 0
for i in range(nbi):
    integrale = integrale + y[i]*(x[i+1]-x[i])
    # dessin du rectangle
    x_rect = [x[i], x[i], x[i+1], x[i+1], x[i]] # abscisses des sommets
    y_rect = [0   , y[i], y[i]  , 0     , 0   ] # ordonnees des sommets
    plt.plot(x_rect, y_rect,"r")
print("integrale =", integrale)
 
plt.show()
 

from scipy import interpolate
import numpy as np
import matplotlib.pyplot as plt
x = np.arange(-15.01, 15.01, 1.00)
y = np.arange(-15.01, 15.01, 1.00)
xx, yy = np.meshgrid(x, y)
z = np.cos(xx**2+yy**2)
f = interpolate.interp2d(x, y, z, kind='quintic')
xnew = np.arange(-15.01, 15.01, 1e-2)
ynew = np.arange(-15.01, 15.01, 1e-2)
znew = f(xnew, ynew)
plt.plot(x, z[0, :], 'ro-', xnew, znew[0, :], 'b-')
plt.show()
 

# -*- coding: utf-8 -*-
"""
Created on Wed Aug 24 16:23:30 2022
 
@author: Ricore
"""
fibonacciList = [0, 1]
# finding 10 terms of the series starting from 3rd term
N = 25
term = 3
while term < N + 1:
    value = fibonacciList[term - 2] + fibonacciList[term - 3]
    fibonacciList.append(value)
    term = term + 1
 
print(fibonacciList)