Bonjour tout le monde;
Voilà, j'ai un petit projet à faire avec l'intégration numérique par plusieurs méthodes.
Je vois bien le principe je dirais analytique (crayon et papier), mais quand à programmer cela j'ai bien plus de mal.
Voilà l'équation, la méthode sur papier que j'ai trouvé et un début de code dont j'aurais besoin de la plus grande aide ^^
L'équation correspond à la première des 3.
Merci d'avance !
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37 # usr/bin/python from pylab import * import numpy as np ####################################### # ORBITES AUTOUR D'UN TROU NOIR DE SCHWARZSCHILD # ####################################### # -------------------------------- # # Methode d'Euler explicite # # -------------------------------- # Rs = 2 # Rayon de Swharzschild m = 1 # Masse quelconque d'une particule c = 1 # Vitesse lumiere r = 0.8 # rayon r E = 1 # Energie particule F = ( (1-(Rs/r))**(-1) ) * E/(m*c**2) def Euler(F, y0, t):
Partager