Bonjour à toutes et à tous,
Voilà j'ai un projet à rendre dans 2 semaines exactement sur la résolution de l'équation de Schrödinger stationnaire à 1D.
Nous avons une méthode qui nous est donnée, sur laquelle nous devons nous baser pour programmer cela pour plusieurs types de potentiels (potentiel plat, Lennard-Jones etc ...).
Néanmoins, j'ai appris la programmation seulement depuis 2 mois et je ne sais mais alors pas du tout comment débuter mon projet ...
Je vous mets le début de mon code :
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 #include <cstdlib> #include <vector> #include "matprod.h" #include <fstream> #define h 6.62606957*pow(10,-34) // Définition de la constante de Planck #define Pi 3.1415926 // Définition de la valeur de Pi // -------------------------------------------------- // // ------ Définition du potentiel carré simple ------ // // -------------------------------------------------- // double V (double x, double xmin, double xmax, double V0) { if(x<xmin || x>xmax) return 0; else return V0; } int main() { cout<<"###########################################################################"<<endl; cout<<"###########################################################################"<<endl; cout<<"########################## Projet Informatique ############################"<<endl; cout<<"############# Résolution de l'équation de Schroedinger à 1D ###############"<<endl; cout<<"###########################################################################"<<endl; cout<<"###########################################################################"<<endl; return 0 ; }
Et ci-joint les quelques documents qui m'ont été donné :
Le pdf : projet M1 ==> Mon sujet correspond au premier
Le pdf : Notice ==> Explication d'une méthode de résolution
projets-M1.pdf
Notice 2014.pdf
Je ne souhaite pas forcément que l'on me donne des codes déjà fait, car je veux coder plus ou moins moi-même afin de comprendre ce que je fait.
En espérant avoir un petit coup de pouce
Merci d'avance !
-----------------------------------------------------
Je viens d'ajouter j'espère de façon correcte au niveau de la syntaxe la fonction potentiel J(x) :
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
38
39
40
41
42
43
44
45
46
47
48 #include <iostream> #include <cmath> using namespace std ; #include <cstdlib> #include <vector> #include "matprod.h" #include <fstream> #define h 6.62606957*pow(10,-34) // Définition de la constante de Planck #define Pi 3.1415926 // Définition de la valeur de Pi // -------------------------------------------------- // // ------ Définition du potentiel carré simple ------ // // -------------------------------------------------- // double V (double x, double xmin, double xmax, double V0) { if(x<xmin || x>xmax) return 0; else return V0; } double J (double V (double x)) { double J = ((2*m)/(pow(h,2)))*(E - V (double x)) ; return J ; } int main() { cout<<"###########################################################################"<<endl; cout<<"###########################################################################"<<endl; cout<<"########################## Projet Informatique ############################"<<endl; cout<<"############# Résolution de l'équation de Schroedinger à 1D ###############"<<endl; cout<<"###########################################################################"<<endl; cout<<"###########################################################################"<<endl; cout << V << endl; return 0 ; }
Partager