1st: you're using incorrect statistics when you change the output of RAND the way you do (setting 0 values to 1 and dividing by the max value). Just use the output of RAND directly in your comparisons, that is what RAND is designed for.
2nd: A function that asks to enter a choice every time it is run is inconsistent with all the other functions in MATLAB (and annoying), the choice should be an input parameter to the function. Also, you could provide default values for the ND and D parameters.
3rd: There is no need to output the result to a figure window. Let the user decide whether they want that or not.
Comment only