杜哈梅积分的matlab程序

cｌccleaｒ％输入数据,开始可以改参数了。ａa=1０;%输入时间长度bb=０.0１;％输入精度%%%％％ｔ＝bb:bb:aa;t1=t;%不用改％%%％％ｔheta=1;%输入荷载频率w=2;%输