# MATLAB: Why use the equationsToMatrix can get the precise result

MATLAB

Why do you use equationsToMatrix to get the root and bring back the equation, which is not exactly 0?
codes are as this
clc;clear all;close all;syms g h d a omega m0 m1 m2 m3 m4 m5 m6 m7 m8 m9;%************************************************************%************************************************************g=9.8;h=20;a=2.25;d=19.17187672;%************************************************************%************************************************************omega=1.000000;m0=0.105131;m1=0.122316;m2=0.297646;m3=0.460332;m4=0.620165;m5=0.778885;m6=0.937054;m7=1.094911;m8=1.252573;m9=1.410105;m10=1.567546;%************************************************************%************************************************************fprintf('omega=%d',omega)syms a_0k a_1k a_2k a_3k a_4k a_5k a_6k a_7k a_8k a_9k;a_2k=1;a_3k=1;a_4k=1;a_5k=1 ;a_6k=1;a_7k=1;a_8k=1;a_9k=1;equation0you =   (2*sinh(d*m0)*(((pi*a_1k*d*m0*sinh(d*m0)*besseli(1, (pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 + d^2*m0^2)) - (2*pi*a_2k*d*m0*sinh(d*m0)*besseli(1, (2*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 + d^2*m0^2)) + (3*pi*a_3k*d*m0*sinh(d*m0)*besseli(1, (3*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 + d^2*m0^2)) - (4*pi*a_4k*d*m0*sinh(d*m0)*besseli(1, (4*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 + d^2*m0^2)) + (5*pi*a_5k*d*m0*sinh(d*m0)*besseli(1, (5*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 + d^2*m0^2)) - (6*pi*a_6k*d*m0*sinh(d*m0)*besseli(1, (6*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 + d^2*m0^2)) + (7*pi*a_7k*d*m0*sinh(d*m0)*besseli(1, (7*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 + d^2*m0^2)) - (8*pi*a_7k*d*m0*sinh(d*m0)*besseli(1, (8*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 + d^2*m0^2)) + (9*pi*a_9k*d*m0*sinh(d*m0)*besseli(1, (9*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 + d^2*m0^2)))/(h*m0*(besselj(1, a*m0) + bessely(1, a*m0)*1i)) + (a*sinh(d*m0)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(2*d^2*h*m0^2*(besselj(1, a*m0) + bessely(1, a*m0)*1i)))*(besselj(0, a*m0) + bessely(0, a*m0)*1i)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(d*m0) - (2*(d/6 - a^2/(4*d)))/d - (2*sin(d*m1)*besselk(0, a*m1)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2)*(((pi*a_1k*d*m1*sin(d*m1)*besseli(1, (pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m1^2)) - (2*pi*a_2k*d*m1*sin(d*m1)*besseli(1, (2*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m1^2)) + (3*pi*a_3k*d*m1*sin(d*m1)*besseli(1, (3*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m1^2)) - (4*pi*a_4k*d*m1*sin(d*m1)*besseli(1, (4*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m1^2)) + (5*pi*a_5k*d*m1*sin(d*m1)*besseli(1, (5*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m1^2)) - (6*pi*a_6k*d*m1*sin(d*m1)*besseli(1, (6*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m1^2)) + (7*pi*a_7k*d*m1*sin(d*m1)*besseli(1, (7*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m1^2)) - (8*pi*a_7k*d*m1*sin(d*m1)*besseli(1, (8*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m1^2)) + (9*pi*a_9k*d*m1*sin(d*m1)*besseli(1, (9*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m1^2)))/(h*m1*besselk(1, a*m1)) - (a*sin(d*m1)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(2*d^2*h*m1^2*besselk(1, a*m1))))/(d*m1) - (2*sin(d*m2)*besselk(0, a*m2)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2)*(((pi*a_1k*d*m2*sin(d*m2)*besseli(1, (pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m2^2)) - (2*pi*a_2k*d*m2*sin(d*m2)*besseli(1, (2*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m2^2)) + (3*pi*a_3k*d*m2*sin(d*m2)*besseli(1, (3*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m2^2)) - (4*pi*a_4k*d*m2*sin(d*m2)*besseli(1, (4*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m2^2)) + (5*pi*a_5k*d*m2*sin(d*m2)*besseli(1, (5*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m2^2)) - (6*pi*a_6k*d*m2*sin(d*m2)*besseli(1, (6*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m2^2)) + (7*pi*a_7k*d*m2*sin(d*m2)*besseli(1, (7*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m2^2)) - (8*pi*a_7k*d*m2*sin(d*m2)*besseli(1, (8*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m2^2)) + (9*pi*a_9k*d*m2*sin(d*m2)*besseli(1, (9*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m2^2)))/(h*m2*besselk(1, a*m2)) - (a*sin(d*m2)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(2*d^2*h*m2^2*besselk(1, a*m2))))/(d*m2) - (2*sin(d*m3)*besselk(0, a*m3)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2)*(((pi*a_1k*d*m3*sin(d*m3)*besseli(1, (pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m3^2)) - (2*pi*a_2k*d*m3*sin(d*m3)*besseli(1, (2*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m3^2)) + (3*pi*a_3k*d*m3*sin(d*m3)*besseli(1, (3*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m3^2)) - (4*pi*a_4k*d*m3*sin(d*m3)*besseli(1, (4*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m3^2)) + (5*pi*a_5k*d*m3*sin(d*m3)*besseli(1, (5*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m3^2)) - (6*pi*a_6k*d*m3*sin(d*m3)*besseli(1, (6*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m3^2)) + (7*pi*a_7k*d*m3*sin(d*m3)*besseli(1, (7*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m3^2)) - (8*pi*a_7k*d*m3*sin(d*m3)*besseli(1, (8*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m3^2)) + (9*pi*a_9k*d*m3*sin(d*m3)*besseli(1, (9*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m3^2)))/(h*m3*besselk(1, a*m3)) - (a*sin(d*m3)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(2*d^2*h*m3^2*besselk(1, a*m3))))/(d*m3) - (2*sin(d*m4)*besselk(0, a*m4)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2)*(((pi*a_1k*d*m4*sin(d*m4)*besseli(1, (pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m4^2)) - (2*pi*a_2k*d*m4*sin(d*m4)*besseli(1, (2*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m4^2)) + (3*pi*a_3k*d*m4*sin(d*m4)*besseli(1, (3*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m4^2)) - (4*pi*a_4k*d*m4*sin(d*m4)*besseli(1, (4*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m4^2)) + (5*pi*a_5k*d*m4*sin(d*m4)*besseli(1, (5*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m4^2)) - (6*pi*a_6k*d*m4*sin(d*m4)*besseli(1, (6*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m4^2)) + (7*pi*a_7k*d*m4*sin(d*m4)*besseli(1, (7*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m4^2)) - (8*pi*a_7k*d*m4*sin(d*m4)*besseli(1, (8*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m4^2)) + (9*pi*a_9k*d*m4*sin(d*m4)*besseli(1, (9*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m4^2)))/(h*m4*besselk(1, a*m4)) - (a*sin(d*m4)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(2*d^2*h*m4^2*besselk(1, a*m4))))/(d*m4) - (2*sin(d*m5)*besselk(0, a*m5)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2)*(((pi*a_1k*d*m5*sin(d*m5)*besseli(1, (pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m5^2)) - (2*pi*a_2k*d*m5*sin(d*m5)*besseli(1, (2*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m5^2)) + (3*pi*a_3k*d*m5*sin(d*m5)*besseli(1, (3*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m5^2)) - (4*pi*a_4k*d*m5*sin(d*m5)*besseli(1, (4*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m5^2)) + (5*pi*a_5k*d*m5*sin(d*m5)*besseli(1, (5*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m5^2)) - (6*pi*a_6k*d*m5*sin(d*m5)*besseli(1, (6*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m5^2)) + (7*pi*a_7k*d*m5*sin(d*m5)*besseli(1, (7*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m5^2)) - (8*pi*a_7k*d*m5*sin(d*m5)*besseli(1, (8*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m5^2)) + (9*pi*a_9k*d*m5*sin(d*m5)*besseli(1, (9*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m5^2)))/(h*m5*besselk(1, a*m5)) - (a*sin(d*m5)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(2*d^2*h*m5^2*besselk(1, a*m5))))/(d*m5) - (2*sin(d*m6)*besselk(0, a*m6)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2)*(((pi*a_1k*d*m6*sin(d*m6)*besseli(1, (pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m6^2)) - (2*pi*a_2k*d*m6*sin(d*m6)*besseli(1, (2*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m6^2)) + (3*pi*a_3k*d*m6*sin(d*m6)*besseli(1, (3*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m6^2)) - (4*pi*a_4k*d*m6*sin(d*m6)*besseli(1, (4*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m6^2)) + (5*pi*a_5k*d*m6*sin(d*m6)*besseli(1, (5*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m6^2)) - (6*pi*a_6k*d*m6*sin(d*m6)*besseli(1, (6*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m6^2)) + (7*pi*a_7k*d*m6*sin(d*m6)*besseli(1, (7*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m6^2)) - (8*pi*a_7k*d*m6*sin(d*m6)*besseli(1, (8*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m6^2)) + (9*pi*a_9k*d*m6*sin(d*m6)*besseli(1, (9*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m6^2)))/(h*m6*besselk(1, a*m6)) - (a*sin(d*m6)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(2*d^2*h*m6^2*besselk(1, a*m6))))/(d*m6) - (2*sin(d*m7)*besselk(0, a*m7)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2)*(((pi*a_1k*d*m7*sin(d*m7)*besseli(1, (pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m7^2)) - (2*pi*a_2k*d*m7*sin(d*m7)*besseli(1, (2*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m7^2)) + (3*pi*a_3k*d*m7*sin(d*m7)*besseli(1, (3*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m7^2)) - (4*pi*a_4k*d*m7*sin(d*m7)*besseli(1, (4*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m7^2)) + (5*pi*a_5k*d*m7*sin(d*m7)*besseli(1, (5*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m7^2)) - (6*pi*a_6k*d*m7*sin(d*m7)*besseli(1, (6*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m7^2)) + (7*pi*a_7k*d*m7*sin(d*m7)*besseli(1, (7*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m7^2)) - (8*pi*a_7k*d*m7*sin(d*m7)*besseli(1, (8*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m7^2)) + (9*pi*a_9k*d*m7*sin(d*m7)*besseli(1, (9*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m7^2)))/(h*m7*besselk(1, a*m7)) - (a*sin(d*m7)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(2*d^2*h*m7^2*besselk(1, a*m7))))/(d*m7) - (2*sin(d*m8)*besselk(0, a*m8)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2)*(((pi*a_1k*d*m8*sin(d*m8)*besseli(1, (pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m8^2)) - (2*pi*a_2k*d*m8*sin(d*m8)*besseli(1, (2*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m8^2)) + (3*pi*a_3k*d*m8*sin(d*m8)*besseli(1, (3*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m8^2)) - (4*pi*a_4k*d*m8*sin(d*m8)*besseli(1, (4*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m8^2)) + (5*pi*a_5k*d*m8*sin(d*m8)*besseli(1, (5*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m8^2)) - (6*pi*a_6k*d*m8*sin(d*m8)*besseli(1, (6*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m8^2)) + (7*pi*a_7k*d*m8*sin(d*m8)*besseli(1, (7*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m8^2)) - (8*pi*a_7k*d*m8*sin(d*m8)*besseli(1, (8*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m8^2)) + (9*pi*a_9k*d*m8*sin(d*m8)*besseli(1, (9*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m8^2)))/(h*m8*besselk(1, a*m8)) - (a*sin(d*m8)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(2*d^2*h*m8^2*besselk(1, a*m8))))/(d*m8) - (2*sin(d*m9)*besselk(0, a*m9)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2)*(((pi*a_1k*d*m9*sin(d*m9)*besseli(1, (pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m9^2)) - (2*pi*a_2k*d*m9*sin(d*m9)*besseli(1, (2*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m9^2)) + (3*pi*a_3k*d*m9*sin(d*m9)*besseli(1, (3*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m9^2)) - (4*pi*a_4k*d*m9*sin(d*m9)*besseli(1, (4*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m9^2)) + (5*pi*a_5k*d*m9*sin(d*m9)*besseli(1, (5*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m9^2)) - (6*pi*a_6k*d*m9*sin(d*m9)*besseli(1, (6*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m9^2)) + (7*pi*a_7k*d*m9*sin(d*m9)*besseli(1, (7*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m9^2)) - (8*pi*a_7k*d*m9*sin(d*m9)*besseli(1, (8*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m9^2)) + (9*pi*a_9k*d*m9*sin(d*m9)*besseli(1, (9*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m9^2)))/(h*m9*besselk(1, a*m9)) - (a*sin(d*m9)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(2*d^2*h*m9^2*besselk(1, a*m9))))/(d*m9)   ;                             equation1you =   2/pi^2 - (2*d*m0*sinh(d*m0)*(((pi*a_1k*d*m0*sinh(d*m0)*besseli(1, (pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 + d^2*m0^2)) - (2*pi*a_2k*d*m0*sinh(d*m0)*besseli(1, (2*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 + d^2*m0^2)) + (3*pi*a_3k*d*m0*sinh(d*m0)*besseli(1, (3*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 + d^2*m0^2)) - (4*pi*a_4k*d*m0*sinh(d*m0)*besseli(1, (4*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 + d^2*m0^2)) + (5*pi*a_5k*d*m0*sinh(d*m0)*besseli(1, (5*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 + d^2*m0^2)) - (6*pi*a_6k*d*m0*sinh(d*m0)*besseli(1, (6*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 + d^2*m0^2)) + (7*pi*a_7k*d*m0*sinh(d*m0)*besseli(1, (7*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 + d^2*m0^2)) - (8*pi*a_7k*d*m0*sinh(d*m0)*besseli(1, (8*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 + d^2*m0^2)) + (9*pi*a_9k*d*m0*sinh(d*m0)*besseli(1, (9*pi*a)/d)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 + d^2*m0^2)))/(h*m0*(besselj(1, a*m0) + bessely(1, a*m0)*1i)) + (a*sinh(d*m0)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(2*d^2*h*m0^2*(besselj(1, a*m0) + bessely(1, a*m0)*1i)))*(besselj(0, a*m0) + bessely(0, a*m0)*1i)*(sinh(2*h*m0)/(4*h*m0) + 1/2)^(1/2))/(pi^2 + d^2*m0^2) - (2*d*m1*sin(d*m1)*besselk(0, a*m1)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2)*(((pi*a_1k*d*m1*sin(d*m1)*besseli(1, (pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m1^2)) - (2*pi*a_2k*d*m1*sin(d*m1)*besseli(1, (2*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m1^2)) + (3*pi*a_3k*d*m1*sin(d*m1)*besseli(1, (3*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m1^2)) - (4*pi*a_4k*d*m1*sin(d*m1)*besseli(1, (4*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m1^2)) + (5*pi*a_5k*d*m1*sin(d*m1)*besseli(1, (5*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m1^2)) - (6*pi*a_6k*d*m1*sin(d*m1)*besseli(1, (6*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m1^2)) + (7*pi*a_7k*d*m1*sin(d*m1)*besseli(1, (7*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m1^2)) - (8*pi*a_7k*d*m1*sin(d*m1)*besseli(1, (8*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m1^2)) + (9*pi*a_9k*d*m1*sin(d*m1)*besseli(1, (9*pi*a)/d)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m1^2)))/(h*m1*besselk(1, a*m1)) - (a*sin(d*m1)*(sin(2*h*m1)/(4*h*m1) + 1/2)^(1/2))/(2*d^2*h*m1^2*besselk(1, a*m1))))/(pi^2 - d^2*m1^2) - (2*d*m2*sin(d*m2)*besselk(0, a*m2)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2)*(((pi*a_1k*d*m2*sin(d*m2)*besseli(1, (pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m2^2)) - (2*pi*a_2k*d*m2*sin(d*m2)*besseli(1, (2*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m2^2)) + (3*pi*a_3k*d*m2*sin(d*m2)*besseli(1, (3*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m2^2)) - (4*pi*a_4k*d*m2*sin(d*m2)*besseli(1, (4*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m2^2)) + (5*pi*a_5k*d*m2*sin(d*m2)*besseli(1, (5*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m2^2)) - (6*pi*a_6k*d*m2*sin(d*m2)*besseli(1, (6*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m2^2)) + (7*pi*a_7k*d*m2*sin(d*m2)*besseli(1, (7*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m2^2)) - (8*pi*a_7k*d*m2*sin(d*m2)*besseli(1, (8*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m2^2)) + (9*pi*a_9k*d*m2*sin(d*m2)*besseli(1, (9*pi*a)/d)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m2^2)))/(h*m2*besselk(1, a*m2)) - (a*sin(d*m2)*(sin(2*h*m2)/(4*h*m2) + 1/2)^(1/2))/(2*d^2*h*m2^2*besselk(1, a*m2))))/(pi^2 - d^2*m2^2) - (2*d*m3*sin(d*m3)*besselk(0, a*m3)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2)*(((pi*a_1k*d*m3*sin(d*m3)*besseli(1, (pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m3^2)) - (2*pi*a_2k*d*m3*sin(d*m3)*besseli(1, (2*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m3^2)) + (3*pi*a_3k*d*m3*sin(d*m3)*besseli(1, (3*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m3^2)) - (4*pi*a_4k*d*m3*sin(d*m3)*besseli(1, (4*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m3^2)) + (5*pi*a_5k*d*m3*sin(d*m3)*besseli(1, (5*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m3^2)) - (6*pi*a_6k*d*m3*sin(d*m3)*besseli(1, (6*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m3^2)) + (7*pi*a_7k*d*m3*sin(d*m3)*besseli(1, (7*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m3^2)) - (8*pi*a_7k*d*m3*sin(d*m3)*besseli(1, (8*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m3^2)) + (9*pi*a_9k*d*m3*sin(d*m3)*besseli(1, (9*pi*a)/d)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m3^2)))/(h*m3*besselk(1, a*m3)) - (a*sin(d*m3)*(sin(2*h*m3)/(4*h*m3) + 1/2)^(1/2))/(2*d^2*h*m3^2*besselk(1, a*m3))))/(pi^2 - d^2*m3^2) - (2*d*m4*sin(d*m4)*besselk(0, a*m4)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2)*(((pi*a_1k*d*m4*sin(d*m4)*besseli(1, (pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m4^2)) - (2*pi*a_2k*d*m4*sin(d*m4)*besseli(1, (2*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m4^2)) + (3*pi*a_3k*d*m4*sin(d*m4)*besseli(1, (3*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m4^2)) - (4*pi*a_4k*d*m4*sin(d*m4)*besseli(1, (4*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m4^2)) + (5*pi*a_5k*d*m4*sin(d*m4)*besseli(1, (5*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m4^2)) - (6*pi*a_6k*d*m4*sin(d*m4)*besseli(1, (6*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m4^2)) + (7*pi*a_7k*d*m4*sin(d*m4)*besseli(1, (7*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m4^2)) - (8*pi*a_7k*d*m4*sin(d*m4)*besseli(1, (8*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m4^2)) + (9*pi*a_9k*d*m4*sin(d*m4)*besseli(1, (9*pi*a)/d)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m4^2)))/(h*m4*besselk(1, a*m4)) - (a*sin(d*m4)*(sin(2*h*m4)/(4*h*m4) + 1/2)^(1/2))/(2*d^2*h*m4^2*besselk(1, a*m4))))/(pi^2 - d^2*m4^2) - (2*d*m5*sin(d*m5)*besselk(0, a*m5)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2)*(((pi*a_1k*d*m5*sin(d*m5)*besseli(1, (pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m5^2)) - (2*pi*a_2k*d*m5*sin(d*m5)*besseli(1, (2*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m5^2)) + (3*pi*a_3k*d*m5*sin(d*m5)*besseli(1, (3*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m5^2)) - (4*pi*a_4k*d*m5*sin(d*m5)*besseli(1, (4*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m5^2)) + (5*pi*a_5k*d*m5*sin(d*m5)*besseli(1, (5*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m5^2)) - (6*pi*a_6k*d*m5*sin(d*m5)*besseli(1, (6*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m5^2)) + (7*pi*a_7k*d*m5*sin(d*m5)*besseli(1, (7*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m5^2)) - (8*pi*a_7k*d*m5*sin(d*m5)*besseli(1, (8*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m5^2)) + (9*pi*a_9k*d*m5*sin(d*m5)*besseli(1, (9*pi*a)/d)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m5^2)))/(h*m5*besselk(1, a*m5)) - (a*sin(d*m5)*(sin(2*h*m5)/(4*h*m5) + 1/2)^(1/2))/(2*d^2*h*m5^2*besselk(1, a*m5))))/(pi^2 - d^2*m5^2) - (2*d*m6*sin(d*m6)*besselk(0, a*m6)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2)*(((pi*a_1k*d*m6*sin(d*m6)*besseli(1, (pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m6^2)) - (2*pi*a_2k*d*m6*sin(d*m6)*besseli(1, (2*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m6^2)) + (3*pi*a_3k*d*m6*sin(d*m6)*besseli(1, (3*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m6^2)) - (4*pi*a_4k*d*m6*sin(d*m6)*besseli(1, (4*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m6^2)) + (5*pi*a_5k*d*m6*sin(d*m6)*besseli(1, (5*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m6^2)) - (6*pi*a_6k*d*m6*sin(d*m6)*besseli(1, (6*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m6^2)) + (7*pi*a_7k*d*m6*sin(d*m6)*besseli(1, (7*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m6^2)) - (8*pi*a_7k*d*m6*sin(d*m6)*besseli(1, (8*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m6^2)) + (9*pi*a_9k*d*m6*sin(d*m6)*besseli(1, (9*pi*a)/d)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m6^2)))/(h*m6*besselk(1, a*m6)) - (a*sin(d*m6)*(sin(2*h*m6)/(4*h*m6) + 1/2)^(1/2))/(2*d^2*h*m6^2*besselk(1, a*m6))))/(pi^2 - d^2*m6^2) - (2*d*m7*sin(d*m7)*besselk(0, a*m7)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2)*(((pi*a_1k*d*m7*sin(d*m7)*besseli(1, (pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m7^2)) - (2*pi*a_2k*d*m7*sin(d*m7)*besseli(1, (2*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m7^2)) + (3*pi*a_3k*d*m7*sin(d*m7)*besseli(1, (3*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m7^2)) - (4*pi*a_4k*d*m7*sin(d*m7)*besseli(1, (4*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m7^2)) + (5*pi*a_5k*d*m7*sin(d*m7)*besseli(1, (5*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m7^2)) - (6*pi*a_6k*d*m7*sin(d*m7)*besseli(1, (6*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m7^2)) + (7*pi*a_7k*d*m7*sin(d*m7)*besseli(1, (7*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m7^2)) - (8*pi*a_7k*d*m7*sin(d*m7)*besseli(1, (8*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m7^2)) + (9*pi*a_9k*d*m7*sin(d*m7)*besseli(1, (9*pi*a)/d)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m7^2)))/(h*m7*besselk(1, a*m7)) - (a*sin(d*m7)*(sin(2*h*m7)/(4*h*m7) + 1/2)^(1/2))/(2*d^2*h*m7^2*besselk(1, a*m7))))/(pi^2 - d^2*m7^2) - (2*d*m8*sin(d*m8)*besselk(0, a*m8)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2)*(((pi*a_1k*d*m8*sin(d*m8)*besseli(1, (pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m8^2)) - (2*pi*a_2k*d*m8*sin(d*m8)*besseli(1, (2*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m8^2)) + (3*pi*a_3k*d*m8*sin(d*m8)*besseli(1, (3*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m8^2)) - (4*pi*a_4k*d*m8*sin(d*m8)*besseli(1, (4*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m8^2)) + (5*pi*a_5k*d*m8*sin(d*m8)*besseli(1, (5*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m8^2)) - (6*pi*a_6k*d*m8*sin(d*m8)*besseli(1, (6*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m8^2)) + (7*pi*a_7k*d*m8*sin(d*m8)*besseli(1, (7*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m8^2)) - (8*pi*a_7k*d*m8*sin(d*m8)*besseli(1, (8*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m8^2)) + (9*pi*a_9k*d*m8*sin(d*m8)*besseli(1, (9*pi*a)/d)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m8^2)))/(h*m8*besselk(1, a*m8)) - (a*sin(d*m8)*(sin(2*h*m8)/(4*h*m8) + 1/2)^(1/2))/(2*d^2*h*m8^2*besselk(1, a*m8))))/(pi^2 - d^2*m8^2) - (2*d*m9*sin(d*m9)*besselk(0, a*m9)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2)*(((pi*a_1k*d*m9*sin(d*m9)*besseli(1, (pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (pi*a)/d)*(pi^2 - d^2*m9^2)) - (2*pi*a_2k*d*m9*sin(d*m9)*besseli(1, (2*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (2*pi*a)/d)*(4*pi^2 - d^2*m9^2)) + (3*pi*a_3k*d*m9*sin(d*m9)*besseli(1, (3*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (3*pi*a)/d)*(9*pi^2 - d^2*m9^2)) - (4*pi*a_4k*d*m9*sin(d*m9)*besseli(1, (4*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (4*pi*a)/d)*(16*pi^2 - d^2*m9^2)) + (5*pi*a_5k*d*m9*sin(d*m9)*besseli(1, (5*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (5*pi*a)/d)*(25*pi^2 - d^2*m9^2)) - (6*pi*a_6k*d*m9*sin(d*m9)*besseli(1, (6*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (6*pi*a)/d)*(36*pi^2 - d^2*m9^2)) + (7*pi*a_7k*d*m9*sin(d*m9)*besseli(1, (7*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (7*pi*a)/d)*(49*pi^2 - d^2*m9^2)) - (8*pi*a_7k*d*m9*sin(d*m9)*besseli(1, (8*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (8*pi*a)/d)*(64*pi^2 - d^2*m9^2)) + (9*pi*a_9k*d*m9*sin(d*m9)*besseli(1, (9*pi*a)/d)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(besseli(0, (9*pi*a)/d)*(81*pi^2 - d^2*m9^2)))/(h*m9*besselk(1, a*m9)) - (a*sin(d*m9)*(sin(2*h*m9)/(4*h*m9) + 1/2)^(1/2))/(2*d^2*h*m9^2*besselk(1, a*m9))))/(pi^2 - d^2*m9^2)    ;    equation0=a_0k-equation0you;equation1=a_1k-equation1you;[A, b] = equationsToMatrix([equation0, equation1], [a_0k ,a_1k]);x = A\b;format longsolution=@(n)(vpa(x(n+1),10));solution(0)solution(1)

>> simplify(A*x-b) ans =  0 0 >> simplify(A*vpa(x)-b) ans =  3.275394916534423828125e-34 + 1.535350322723388671875e-34i - 2.019959735870361328125e-34 - 2.6699886322021484375e-35i