Hi Rick,

I tried to run the following code for making an integral but I got error msg, if you run it you will see that I got same error msg as you specified in this post. Can you shed some light? Thanks!

proc iml;

start grand(u) global(a1,a2,a3,a4,a5,a6);
grand=u##(a1-1)#(1-u)##(a4-a1-1)#(1-u#a5)##(-a2)#(1-u#a6)##(-a3);
return(grand);
finish;

start Appell(a1,a2,a3,a4,a5,a6);
call quad(z,"grand",{0 1}) peak=0.1234;
beta=fact(a4-1)/(fact(a1-1)#fact(a4-a1-1));
f=beta#z;
return(f);
finish;

test=Appell(186,464,-46,259,0.89,0.9879);
print test;
quit;