Use MPI to parallelize the following code:

#include #include #include /* Fonction calculant le second membre des equations du mvt ici d2u/dt2 = -(EI/(rhoS))*d4u/dx4 + NL avec NL = (1/2)*(E/rho) . (du/dx)^2 . d2u/dx2 */ double secondMembre(int N,double dt,double dx,double E, double rho, double S, double I, double F, double *u, double *SM) { int i; double du,d2u,d3u,d4u,NL,lin; double clin, cnl, cm,cmax=0.,maxu=0.; // application des CL // bord encastre u = 0, du/dx = 0 u[0] = 0.; u[1] = 0.; SM[0]= 0.; SM[1]= 0.; // bord encastre u[N-2] = 0.; u[N-1] = 0.; SM[N-2]= 0.; SM[N-1]= 0.; // on ne resoud qu'a l'interieur du domaine, les autres points sont les CL for(i=2;i