matlab求杜哈梅积分