エアリーディスクの中央ならばかんたんだろう。
2πr*cos(位相誤差(r))を0からRまでの積分をπR^2で割る。
回折像は開口をフーリエ変換すれば求められるので。