Discussion :

[Andromedae]

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):

[marco056]

Il y a pléthore de solutions grâce à moteur de recherche basique.

[Andromedae]

J'ai déjà tapé Euler sur ce forum, 5 sujets, aucuns ne m'aident réellement :/

[marco056]

J'ai déjà tapé Euler sur ce forum, 5 sujets, aucuns ne m'aident réellement :/

Pas très perspicace :
https://www.google.fr/search?client=...D9Ss8wft_4SQBA
Il suffit de taper "méthode euler python" et il tombe des dizaines de réponses pertinentes (les solutions sont même données !).

[Andromedae]

Tu crois que je ne l'ai pas fait ?
J'attendais juste d'avoir quelques pistes en rapport avec mes équations ... qui ne sont pas y' = y ...

[marco056]

C'est sûr : si tu attends que quelqu'un fasse tes devoirs, tu risques d'attendre longtemps.
Propose un bout de code, montre un peu de bonne volonté et tu auras sans doute de l'aide ici.