用Matlab计算静刚度

[柏拉图的世界]

现在在做毕业设计，牵扯到用Matlab计算静刚度。很多地方不熟悉，各位大牛救个急啊~
比如说：
Dw=5
ri=7
cos_tau=(1/ri)/((4/Dw)-(1/ri))
syms k phi
E=sym(int(sqrt(1-(k*sin(phi))^2),phi,0,pi/2))
K=sym(int(1/(sqrt(1-(k*sin(phi))^2)),phi,0,pi/2))
k=solve('((2-k^2)*E-2*(1-k^2)*K)/(k^2*E)-cos_tau')
k_ans=subs(k)

其中，E和K是只与k有关的椭圆积分，k、E、K又满足((2-k^2)*E-2*(1-k^2)*K)/(k^2*E)=cos_tau。所以想先用k表达E、K，然后代入那个等式，解得k的数值。但是我试了半天最后按上面的程序只能得到以k的椭圆积分表达k的结果，而不能得到k的数值。怎么办啊？
很着急啊，这些算不出来毕设就很难进行下去了。高手一定帮帮忙啊!~
PS：我连数值计算和符号计算的指令都弄不太清楚。