smax = 0.7838531572;
steps=smax/1000;
n=1000;
x=zeros(0,n);
y=zeros(0,n);
c=0;
p=0;
counter=0;

for b=0:steps:smax
    counter=counter+1;
    p = @(s) 0.06598252621*exp(((-0.06222716681)./(s.^3)+0.1292039801))./((((0.01250766053)./(s.^3))-0.02597).^0.4.*s.^4);
    c= integral(p,0,b);
    x(counter)=b;
    y(counter)=c;
end

plot(x,y)