IRM matlab

. Salut les jeunes, je suis confronté à un problème matlab
Voici le sujet: Avec un certain type de séquence, on peut montrer que le signal délivré par une IRM est proportionnel à p[-exp(-TR/T1)]exp(-TE/T2), où p est la densité de proton, T1 le temps de relaxation spin-réseau et T2 le temps de relaxation spin-spin. Ces trois quantité sont des fonctions spatiales dépendant des tissus dans lesquels se trouvent l'eau. TE et TR sont deux paramètre ficés par l'opérateur. Ecrire un programme permettant de simuler la dépendance du contraste en fonction des différents paramètres (images pondérées en densité, en T1 ou en T2)
Ci des passionnés pouvaient me réaliser ce progamme, j'en serais très reconnaissant.