MATLAB含参数的二重积分画图

推出式子：

当a一定，为0.5度即π/360rad；r1=40mm,r2=60mm；看e与τ、d的关系，其中τ从0°到180°变化，d从0mm-40mm变化,
大概得出的图像应该是

如何通过matlab编程序得到啊？
在计算的过程中我发现τ的取值会影响整个式子，能不能直接得出这样的图啊?

同学你好，把你积分的代码贴一贴，这里光看公式也不知道你的参数具体怎么处理的

想试一下这样固定τ，然后代入d看值发现也不行，我最后想出来的还是一个大图，不知道怎么算。