Diagramme en bâton - scilab

[bonjourrr]

Bonjour, je dois écrire un algorithme traçant un diagramme en bâtons représentant la loi de Bernoulli de paramètre 1/2.
J'ai essayé avec

function x=loidebernoulli(p)
[x]=grand(1,100,'bin',1,p)
plot2d3([100;1],[x])
endfunction

A vrai dire, je suis un peu perdu. J'ai demandé conseil et on m'a dit qu'il fallait utilisé "rand" plutôt que "grand"... Pouvez-vous m'aider ?