Hi,
I tried to define the skew normal distribution function as you have suggested. I found the way you define the function is not appropriate, so I tried myself and rewrote the function y as following:
y=1/(w*PI)*exp(-(x-u)^2/2/w^2)*integral(myintegrand, -inf ,a*(x-u)/w ,a ,u ,w)
where you should have the following script as your definition for integrand “myintegrand”:
return exp(-t^2/2);
There is no “*dt” after the expression but should have “return” in front.
Can you try these steps and see whether the problem gets solved? If it persists, please do not hesitate to ask questions.
Thanks!
Sean
OriginLab Tech. Service
