why use the equationsToMatrix can get the precise result?

2 visualizzazioni (ultimi 30 giorni)
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 long
solution=@(n)(vpa(x(n+1),10));
solution(0)
solution(1)

Risposta accettata

Walter Roberson
Walter Roberson il 11 Mag 2019
Your equations involve the irrational number Pi. No finite truncation of the results in x can possibly exactly match the infinitely long Pi.
>> 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
  2 Commenti
dcydhb dcydhb
dcydhb dcydhb il 11 Mag 2019
in fact my question has 9 equations,and i can't paste it as it is too long,in fact the way
simplify(A*x-b)
simplify(A*vpa(x)-b)
are great and in this way i can get the 9 equations Checking result
about e-34
but when certify it in the convensional way
a_0k=...
a_1k =...
equation0you = ...
equation1you = ...
equation0=a_0k-equation0you
equation1=a_1k-equation1you
and in this traditional way of 9 equations,
the certified result is as this
equation0you =
1.351151124950071e-01 + 2.424118832106468e-01i
equation1you =
6.809596339796679e-02 - 6.967954118599076e-02i
equation2you =
-7.642085494819821e-03 + 2.164658528943209e-02i
equation3you =
2.369729764183772e-03 - 9.895950999635576e-03i
equation4you =
-1.126122185829983e-03 + 5.576256528491903e-03i
equation5you =
6.661879693648914e-04 - 3.565578348703271e-03i
equation6you =
-4.495053393553626e-04 + 2.482327854624580e-03i
equation7you =
3.326572160341770e-04 - 1.840839544059650e-03i
equation8you =
-3.714702983427300e-04 + 2.010614845948146e-03i
equation9you =
2.528750271046693e-04 - 1.238731056260830e-03i
>>
completely different from the
about e-34
why?

Accedi per commentare.

Più risposte (0)

Categorie

Scopri di più su Mathematics in Help Center e File Exchange

Tag

Prodotti


Release

R2014a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by