T=1; // an
S_0=100; 
r=0.05; // par an
sigma=0.3; // par racine d'annee sigma^2 * T  est sans dimension
K=100;

// Question 1
N=1000;
W_T=sqrt(T)*rand(1,N,"gauss");
S_T=...

clf;
histplot(100,W_T);

clf;
histplot(100,S_T);