Using R language:

Generate 5000 random samples of x.

Use below function to generate a set of y values given the model parameter a and the generated input data x. Plot y against x.

output = function(a,x){ y=-1 +a*x

return(y) }