clear all
d=1e-50;
t=[-1e-48:d:1e-48];
c=length(t);
for i=1:c
y(i) = log_pnorm(t(i),0, 1e-100);
end

for i=1:c
    z(i) = exp(y(i));
end

sum(z)*d
plot(t,z)