求解有关自己写BS代码的问题

代码：如何改动可以算出C，delta,vega 的值，用自己输入的数值

1. clear all;
   
2. syms S0 k r T sigma
   
3. d1=(log(S0/k)+(r+0.5*sigma^2)*T)/(sigma*sqrt(T));
   
4. d2=(log(S0/k)+(r-0.5*sigma^2)*T)/(sigma*sqrt(T));
   
5. C=S0*normcdf(d1)-k*exp(-r*T)*normcdf(d2);
   
6. dCdS0=diff(C,S0);
   
7. dCdsigma=diff(C,sigma);
   
8. S0=input('请输入S0=');
   
9. k=input('请输入k=');
   
10. r=input('请输入r=');
    
11. T=input('请输入T=');
    
12. sigma=input('请输入sigma=');
    
13. d1=subs(d1);
    
14. d2=subs(d2);
    
15. C=subs(C)
    
16. delta=subs(dCdS0)
    
17. vega=subs(dCdsigma)

复制代码

这段代码运行出错结果：

1. 错误使用 symfun>validateArgNames (line 211)
   
2. Second input must be a scalar or vector of unique symbolic variables.
   
3. 
   
4. 出错 symfun (line 45)
   
5.             y.vars = validateArgNames(inputs);
   
6. 
   
7. 出错 sym/subsasgn (line 762)
   
8.                 C = symfun(B,[inds{:}]);
   
9. 
   
10. 出错 normcdf>localnormcdf (line 100)
    
11.     p(sigma==0 & x<mu) = 0;
    
12. 
    
13. 出错 normcdf (line 46)
    
14. [varargout{1:max(1,nargout)}] = localnormcdf(uflag,x,varargin{:});
    
15. 
    
16. 出错 Untitled6 (line 5)
    
17. C=S0*normcdf(d1)-k*exp(-r*T)*normcdf(d2);

复制代码