# TP 3 ANALYSE NUMERIQUE

# notations pour la quadrature :  
# f : fonction à intégrer
# a,b : bornes d'intégration
# n : nombre d'intervalles de décomposition de l'intervalle [a,b]
# 
# notations pour la résolution de l'équation différentielle y'=f(t,y) :
# a,b : bornes de l'intervalle sur lequel on veut calculer la solution y
# ya : valeur de y en a
# n : nombre de points du schéma numérique
# yy : liste contenant les valeurs de y

import sys
from math import *
from pylab 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