Use MPI to parallelize the following code:

#include  #include  #include  // fonction qui calcule l'equation a resoudre // Au = [Id -k d2/dx2].u // par difference finies void equation(int N,double dx,double *k, double *u,double *Au) { int i; double d2udx2; for(i=1;imaxx) { maxx = fabs(u[i]); } } return maxx; } // schema iteratif de richardson // par le minimum du residu void richardson(int N, double dx,double *k, double *uj,double *smb,int itmax) { int i,j; double alpha, numer, denom,res; double rj[N],Arj[N]; // rj = erreur, Arj = L(rj) // calcul du residu initial equation(N,dx,k,uj,rj); //rj = erreur for(i=0;i1e-6 && j