matlab解贝塞尔函数

matlab解贝塞尔函数
能帮我解这个线性贝塞尔方程吗 可化为矩阵{H}*Ai=0 Ai不为零 求{H}=0的解 贝塞尔函数可以全取为零阶 即n=0 希望得出圆频率w与a,b,v1,v2,h1,h2,E1,E2,p1,p2的关系式

syms a k b v1 v2 D1 D11 D12 k1
>> H11 =besselj(0,k*a);
>> H12 =bessely(0,k*a);
>> H13 =besseli(0,k*a);
>> H14 =besselk(0,k*a);
>> H15=0;H16=0;
>> H21=diff(H11);
>> H22=diff(H12);
>> H23=diff(H13);
>> H24=diff(H14);
>> H25=0;
>> H26=0;
>> H31 =besselj(0,k*b);
>> H32 =bessely(0,k*b);
>> H33 =besseli(0,k*b);
>> H34 =besselk(0,k*b);
>> H35 =-besselj(0,k1*b);
>> H36 =-besseli(0,0.k1*b);
>> H41 =diff(H31);
>> H42 =diff(H32);
>> H43 =diff(H33);
>> H44 =diff(H34);
>> H45 =diff(H35);
>> H46 =diff(H36);
>> H51=-D1*diff((diff(H31,b,2)+1/b*diff(H31,b)),b);
>> H52=-D1*diff((diff(H32,b,2)+1/b*diff(H32,b)),b);
>> H53=-D1*diff((diff(H33,b,2)+1/b*diff(H33,b)),b);
>> H54=-D1*diff((diff(H34,b,2)+1/b*diff(H34,b)),b);
>> H55=(D11+D12)*diff((diff(-H35,b,2)+1/b*diff(-H35,b)),b);
>> H56=(D11+D12)*diff((diff(-H36,b,2)+1/b*diff(-H36,b)),b);
>> H61=-D1*(V1/b*diff(H31,b)+diff(H31,b,2));
>> H62=-D1*(V1/b*diff(H32,b)+diff(H32,b,2));
>> H63=-D1*(V1/b*diff(H33,b)+diff(H33,b,2));
>> H64=-D1*(V1/b*diff(H34,b)+diff(H34,b,2));
>> H65=(D11+D12)*diff(-H35,b,2)+(D11*V1+D12*V2)/b*diff(-H35,b);
>> H66=(D11+D12)*diff(-H36,b,2)+(D11*V1+D12*V2)/b*diff(-H36,b);
>> A=[H11 H12 H13 H14 H15 H16; H21 H22 H23 H24 H25 H26; H31 H32 H33 H34 H35 H36; H41 H42 H43 H44 H45 H46; H51 H52 H53 H54 H55 H56; H61 H62 H63 H64 H65 H66]
solve(det(A));