# -*- coding: utf-8 -*-
"""
Created on Sun Jun 15 14:57:13 2025

@author: franc
"""

import matplotlib.pyplot as plt
import math
import scipy.integrate as integr


def phi(t):
    return math.exp(-1/t**2)
    
lx=[-3.001+6/100*k for k in range(101)]
ly=[phi(x) for x in lx]
plt.plot(lx,ly)

# Ci-dessous, une fonction retournant f(0).
# Pour retourner b, il suffit de l'appliquer  à  a/2.
def Q2d(a):
    I=0
    depart=0
    while I<a:
        increm=integr.quad(phi,depart,depart+0.1)[0]
        depart+=0.1
        I+=increm
    return depart