% Part (b):
z = (2000-2200)/230; 
p = 1 - normcdf( z );

Ps = binopdf( 2, 3, p ) + binopdf( 3, 3, p )