Cette thèse porte sur l'approximation des équations aux dérivées partielles, en particulier l'équation de diffusion en tomographie optique. Elle peut se présenter en deux parties essentielles. Dans la première partie on discute le problème direct alors que le problème inverse est abordé dans la seconde partie. Pour le problème direct, on suppose que les paramètres optiques et les fonctions sources sont donnés. On résout alors le problème de diffusion dans un domaine où la densité du flux lumineux est considérée comme une fonction inconnue à approcher numériquement. Le plus souvent, pour reconstruire le signal numérique dans ce genre de problème, une discrétisation dans le temps est nécessaire. Nous avons proposé d'utiliser la transformée de Fourier et son inverse afin d'éviter une telle discrétisation. Les techniques que nous avons utilisées sont la quadrature de Gauss-Hermite ainsi que la méthode de Galerkin basée sur les B-splines ou les B-splines tensorielles ainsi que sur les fonctions radiales. Les B-splines sont utilisées en dimension un alors que les B-splines tensorielles sont utilisées lorsque le domaine est rectangulaire avec un maillage uniforme. Lorsque le domaine n'est plus rectangulaire, nous avons proposé de remplacer la base des B-splines tensorielles par les fonctions à base radiale construites à partir d'un nuage de points dispersés dans le domaine. Du point de vue théorique, nous avons étudié l'existence, l'unicité et la régularité de la solution puis nous avons proposé quelques résultats sur l'estimation de l'erreur dans les espaces de type Sobolev ainsi que sur la convergence de la méthode. Dans la seconde partie de notre travail, nous nous sommes intéressés au problème inverse. Il s'agit d'un problème inverse non-linéaire dont la non-linéarité est liée aux paramètres optiques. On suppose qu'on dispose des mesures du flux lumineux aux bords du domaine étudié et des fonctions sources. On veut alors résoudre le problème inverse de façon à simuler numériquement l'indice de réfraction ainsi que les coefficients de diffusion et d'absorption. Du point de vue théorique, nous avons discuté certains résultats tels que la continuité et la dérivabilité, au sens de Fréchet, de l'opérateur mesurant le flux lumineux reçu aux bords. Nous avons établi les propriétés lipschitzienne de la dérivée de Fréchet en fonction des paramètres optiques. Du point de vue numérique nous nous somme intéressés au problème discret dans la base des B-splines et la base des fonctions radiales. En suite, nous avons abordé la résolution du problème inverse non-linéaire par la méthode de Newton-Gauss.