syms z
iv=g*K/(b*s0)-a*g/b;
eq1=(exp(g*z)-1)*exp(h*z.^2/2)-iv
solve(eval(eq1))


除了z 其他变量都有数值，但是需要迭代数值，所以不能直接写在solve里面，这样计算出来的结果Z是符号变量，
我要用它计算下面这个函数
  P=normcdf(z,0,1);
但是始终都有问题，已经试用了double 等常用函数了，都不可以，求大神指教呀？？？

比如说，结果是这样的
>> syms z
>> eq1 =exp(z^2/200000)*(exp((23*z)/100) - 1) + 299/490
eq1 =
exp(z^2/200000)*(exp((23*z)/100) - 1) + 299/490
>> z=solve(eval(eq1))
z =
z = -4.0956551315107424897773903470668
错误：
P=normcdf(z,0,1);
??? Error using ==> erfcore
Input must be single or double.

Error in ==> erfc at 23
y = erfcore(x,1);

Error in ==> normcdf at 68
p = 0.5 * erfc(-z ./ sqrt(2));

double(z)
??? Error using ==> sym.double at 29
DOUBLE cannot convert the input expression into a double array.
If the input expression contains a symbolic variable, use the VPA function instead.