matlab运行了很长时间,很简单的一个matlab程序运行时间过长，得不到结果

主程序：

clc

clear all

global kesai omega0 m rou D a1 a3 K C %%主系统基本参数

global epsilong_T s_33 C_pzt C_pzt1 K_pzt K_pzt1 k0 d_33 s_0 l_0 theta L R %压电材料参数

kesai=0.0013;m=0.44;omega0=62.83;rou=1.2;D=0.015;

a1=2.3;a3=-18;

K=1736.947916;

C=0.07187752;

s_0=0.005*0.005*pi; %cross section area of the PZT rod

l_0=0.01; %length of the PZT rod

epsilong_T=3400*8.854e-12; % permittivity under constant stress 常应力下的介电常数

s_33=20.7e-12; %compliance of the piezoceramic under constant electric 恒定电场下的弹性柔顺系数

d_33=593e-12; %压电常数

C_pzt=epsilong_T*s_0/l_0; %capacitance between the electrodes of the PZT rod with no external force无外力下的电容

K_pzt=(1/s_33)*s_0/l_0; %stiffness of the short-circuited PZT rod 短路刚度

k0=d_33*(1/(sqrt(epsilong_T*s_33))); % electromechanical c