Aide table loi normale Scilab [Résolu]

   // La fonction f

function y=f(x)                   // Création de la function à intégrer.
    y= 1/sqrt(2*%pi)*exp(-x^2/2) // La fonction f(x) de l'énoncé.
endfunction                     // Fin de la function.
x=linspace(-4,4,100)    // On définit les paramétres graphique.
clf                    // Réinitialise les valeurs par défault.
plot(x,f)             // Affichage de la courbe.

// Une fonction qui calcule une approximation de l'intégrale par la méthode de Simpson en prenant pour n l'entier calculé à la partie A. 

function ms=simpson(a,b,n) // Création de la function simpson.
    h=(b-a)/n             //Début de la formule de Simpson.
    ms=0
    for k=0:n-1          // Début de la boucle.
        x1=a+k*h
        x2=x1+h
        ms=ms+f(x1)+4*f((x1+x2)/2)+f(x2)
    end                 // Fin de la boucle.
    ms=ms*h/6
endfunction

format(7)
table=[]
table(1,1)="x"
for i=2:41
    table(i,1)=string((i-2)/10)
end
for j=2:11
    table(1,j)=string((j-2)/100)
end
for i=2:41
    for j=2:11
        pij=0.5+simpson(0,(i-2)/10+(j-2)/100,9)
        table(i,j)=string(pij)
    end
end
disp(table)