# TP 3 ANALYSE NUMERIQUE

# notations pour la quadrature :  
# f : fonction a integrer
# a,b : bornes d'integration
# n : nombre d'intervalles de decomposition de l'intervalle [a,b]
# 
# notations pour la resolution de l'equation differentielle y'=f(t,y) :
# a,b : bornes de l'intercalle sur lequel on veut calculer la solution y
# ya : valeur de y en a
# n : nombre de points du schema numerique
# yy : liste contenant les valeurs de y

from math import *
from sys import *

def trapeze_composite(f,a,b,n):
      h = (b-a)*1.0/n
      s = 0
      for i in range(0,n):
            s += (f(a+i*h)+f(a+(i+1)*h))*0.5
      return h*s

def simpson_composite(f,a,b,n):
      return 0

def euler_progressif(f,a,b,ya,n):
      h = (b-a)*1.0/n
      yy = [ya]
      for i in range(0,n):
            yy.append(yy[i]+h*f(a+i*h,yy[i]))
      return yy

def euler_modifie(f,a,b,ya,n):
      return 0