圆形矢量airy beam

已知圆形矢量涡旋airy beam有自动聚焦效应
他在源平面上的光场表达式：

a.源平面光强图c.焦点处光强图
这里的r0=1mm，b =0.1，w=0.5mm，a=0.05仅供参考
最后得到的图a，c分别是

源平面光强图本人已经画出，主要问题是代码传输到焦点处的时候无法聚焦
求一个复现图a，c的代码