请教一下哪位高手Matlab中如果已知一个指数形式的概率密度函数p（x）,怎么求得其概率分布函数直方图?

clear;
clc;
a=20;b=2;c=150;
%概率密度函数%
for x=0:255
if x>c
p(x+1)=exp((x-c)/b)/(a+b);
else
p(x+1)=exp((x-c)/a)/(a+b);
end
end
p=p/sum(p); %归一化
total=0;
%概率分布函数%
for x=0:255
F(x+1)=total+p(x+1);
total=F(x+1);
end
bar(0:255,F)axis([0,255,0,1])