MU = 70;
SIGMA = 8.2;

% Define random variable x (speed of cars on highway)

x=50:1:90;

% Define the function of the random variable x (PDF function)

fx = 1/(SIGMA * sqrt(2*pi)) * exp(-1/2 * ((x-MU)/SIGMA).^2 );

% Make a plot
plot(x,fx,'o--r')
xlabel('Random Variable')
ylabel('Value of F(x)')
grid

% calculate the area under the curve

area_Under_Curve= trapz(x,fx);

  disp(' ')
  disp([blanks(5),'Area under the Normal Distribution : '])
  disp(' ')
  disp([blanks(5),'     P(x>75)=',num2str(area_Under_Curve)])
  disp(' ')