元のプログラム（１)使用環境の変化により修正を繰り返しています。
（２）VBA（Visual Basic Application）から呼び出すDLL（ダイナミックリンクライブラリ） ４章：光線追跡計算プログラムのソース公開と修正をベースにしました。
光線追跡計算プログラムのダウンロード光線追跡計算プログラムは以下からダウンロードできます。
ダウンロード後は解凍してから使用してください。
光線追跡計算プログラム]をダウンロードする。解凍すると「光線追跡」フォルダーがあります。
「光線追跡」フォルダー内に
（１）光線追跡実行ファイルフォルダー
（２）光線追跡.sln
（３）光線追跡フォルダー
（４）.vsフォルダー
があります。
注（１）光線追跡.slnを起動して、ビルドすると「光線追跡.VC.db」「光線追跡.VC.VC.opendb」「ipchフォルダー 」「Debugフォルダー 」「Releaseフォルダー 」が生成されますが容量が大きいので削除しました。
注（2）ソースコードは「光線追跡」フォルダーの「光線追跡.cpp」ファイルのみに記載されています。
実行ファイル群「光線追跡実行ファイル 」フォルダー内に以下のファイルがあります。
（１）光線追跡.exe：コンソールアプリ
（２）VBA_C_W10光線追跡.xls：VBAアプリ
（３）flag.txt：コントロールフラグファイル
（４）IN_FM.txt：計算条件ファイル
（４）IN_FM.txt2：計算条件ファイル
（５）OUT_M.xls：計算結果ファイル
注（1） 「光線追跡実行ファイル 」フォルダーの保存場所の制限は無く、任意の場所に保存して実行できます。
操作方法（１）近接露光.exeをダブルクリックします。
（２）コンソールウインドウが表示されます。（文字等の表示はありません。）
（３）VBA_C_W10光線追跡.xlsをダブルクリックします。
注（１） VBA_C_W10光線追跡の操作はかなり複雑であり、使いこなすのは試行錯誤が必要となります。
注（２）条件設定は、光の不思議と応用の26章〜38章を参照する必要があります。
注（３）「IN_FM」と「IN_FM2」は「 IN_FM (演習1）」と「 IN_FM2 (演習1） 」の値がデホルトで設定されています。これは、光の不思議と応用32章の図32-2 タイプA等倍光学系の配置となっています。
（４）「スポットダイアグラム計算実行」ボタンを押します。
（５）スポット要約シートが表示されます。
（６）完全等倍でテレセントリックであることが理解できます。
（７）光線図シートをクリックします。
（８）光線が通過する範囲が理解できます。
（９）シート「スポット図」をクリックします。
（１０）シート「球面レンズ作図（連動用）」をクリックします。
（１１）「球面レンズ作図実行」ボタンを押します。
（１２）シート「作図」が表示されます。
（１３）シート「操作」の「波面収差計算実行」ボタンを押します。
（１４）シート「OUT_M」に波面収差計算結果が表示されます。
*これから最大収差は波長656.3nmの円周方向に発生することがわかります。
| V[0] | Y[NUMB-1] | V[NUMB-1] | AB |
| -5.00E-02 | 8.06E-03 | 5.00E-02 | 3.50E-01 |
| -4.38E-02 | 7.48E-03 | 4.38E-02 | 2.76E-01 |
| -3.75E-02 | 6.75E-03 | 3.75E-02 | 2.08E-01 |
| -3.13E-02 | 5.89E-03 | 3.13E-02 | 1.48E-01 |
| -2.50E-02 | 4.89E-03 | 2.50E-02 | 9.66E-02 |
| -1.88E-02 | 3.78E-03 | 1.88E-02 | 5.52E-02 |
| -1.25E-02 | 2.58E-03 | 1.25E-02 | 2.48E-02 |
| -6.25E-03 | 1.30E-03 | 6.25E-03 | 6.24E-03 |
| 0.00E+00 | 0.00E+00 | 0.00E+00 | 8.66E-11 |
| 6.25E-03 | -1.30E-03 | -6.25E-03 | 6.24E-03 |
| 1.25E-02 | -2.58E-03 | -1.25E-02 | 2.48E-02 |
| 1.88E-02 | -3.78E-03 | -1.88E-02 | 5.52E-02 |
| 2.50E-02 | -4.89E-03 | -2.50E-02 | 9.66E-02 |
| 3.13E-02 | -5.89E-03 | -3.13E-02 | 1.48E-01 |
| 3.75E-02 | -6.75E-03 | -3.75E-02 | 2.08E-01 |
| 4.38E-02 | -7.48E-03 | -4.38E-02 | 2.76E-01 |
| 5.00E-02 | -8.06E-03 | -5.00E-02 | 3.50E-01 |
| 0次 | -4.90E-07 | ||
| 1次 | -1.65E-16 | ||
| 2次 | 4.00E-01 | ||
| 3次 | 5.10E-16 | ||
| 4次 | -5.56E-02 | ||
| 5次 | -3.31E-16 | ||
| 6次 | 6.06E-03 |
*上記表は波長656.3nmの円周方向の計算結果です。
*最大収差はNA0.05の位置で0.35λです。
*収差の分析結果は２次の焦点収差で0.4λです。
（１５）シート「レンズ間隔最適化計算実行」ボタンを押します。
（１６）シート「OUT_M」にレンズ間隔最適化計算結果が表示されます。
N次元Gauss-Newton法によるレンズ間隔(D)最適化計算結果
| No | E | Xa | Xs | Ya | Ys | D(0) |
| 1 | 0.003778 | 8.98482E-07 | 0.00382 | 0 | 0 | 96.247 |
| 2 | 0.003570 | 5.4499E-07 | 0.00361 | 0 | 0 | 96.259 |
| 3 | 0.003547 | 4.27119E-07 | 0.00359 | 0 | 0 | 96.263 |
| 4 | 0.003544 | 3.87847E-07 | 0.00359 | 0 | 0 | 96.264 |
| 5 | 0.003544 | 3.74764E-07 | 0.00359 | 0 | 0 | 96.265 |
| 6 | 0.003544 | 3.70416E-07 | 0.00359 | 0 | 0 | 96.265 |
| 7 | 0.003544 | 3.68953E-07 | 0.00359 | 0 | 0 | 96.265 |
| 8 | 0.003544 | 3.68492E-07 | 0.00359 | 0 | 0 | 96.265 |
| 9 | 0.003544 | 3.68358E-07 | 0.00359 | 0 | 0 | 96.265 |
| 10 | 0.003544 | 3.68282E-07 | 0.00359 | 0 | 0 | 96.265 |
*軸間距離D(0)＝96.265mmが最適であることがわかります。
（１７）シート「IN_FM」のD18セルの値を96.265に変更します。
（１８）シート「操作」 の「スポットダイアグラム計算実行」ボタンを押します。
（１９）シート「スポット図」をクリックします。
*スポット範囲が小さくなっているのが確認できます。
（２０）シート「操作」の「波面収差計算実行」ボタンを押します。
（２１）シート「OUT_M」に波面収差計算結果が表示されます。
*波長656.3nmの円周方向の収差は以下に変化します。
| V[0] | Y[NUMB-1] | V[NUMB-1] | AB |
| -5.00E-02 | 6.26E-03 | 5.00E-02 | 2.82E-01 |
| -4.38E-02 | 5.90E-03 | 4.38E-02 | 2.24E-01 |
| -3.75E-02 | 5.40E-03 | 3.75E-02 | 1.70E-01 |
| -3.13E-02 | 4.76E-03 | 3.13E-02 | 1.21E-01 |
| -2.50E-02 | 3.99E-03 | 2.50E-02 | 7.94E-02 |
| -1.88E-02 | 3.11E-03 | 1.88E-02 | 4.55E-02 |
| -1.25E-02 | 2.13E-03 | 1.25E-02 | 2.05E-02 |
| -6.25E-03 | 1.08E-03 | 6.25E-03 | 5.16E-03 |
| 0.00E+00 | 0.00E+00 | 0.00E+00 | 8.66E-11 |
| 6.25E-03 | -1.08E-03 | -6.25E-03 | 5.16E-03 |
| 1.25E-02 | -2.13E-03 | -1.25E-02 | 2.05E-02 |
| 1.88E-02 | -3.11E-03 | -1.88E-02 | 4.55E-02 |
| 2.50E-02 | -3.99E-03 | -2.50E-02 | 7.94E-02 |
| 3.13E-02 | -4.76E-03 | -3.13E-02 | 1.21E-01 |
| 3.75E-02 | -5.40E-03 | -3.75E-02 | 1.70E-01 |
| 4.38E-02 | -5.90E-03 | -4.38E-02 | 2.24E-01 |
| 5.00E-02 | -6.26E-03 | -5.00E-02 | 2.82E-01 |
| 0次 | -4.88E-07 | ||
| 1次 | -6.03E-16 | ||
| 2次 | 3.31E-01 | ||
| 3次 | 2.40E-15 | ||
| 4次 | -5.55E-02 | ||
| 5次 | -1.90E-15 | ||
| 6次 | 6.04E-03 |
*上記表は波長656.3nmの円周方向の計算結果です。
*最大収差はNA0.05の位置で0.282λです。
*収差の分析結果は２次の焦点収差で0.331λです。
（２２）「コンソール終了」ボタンを押すとflag.txt（コントロールフラグファイル）に文字コード‘Z’が書き込まれます。
（２３）文字コード‘Z’を認識して光線追跡.exe（コンソールアプリ）が終了します。
（２４）最後にVBA_C_W10光線追跡.xls（VBAアプリ）を終了します。
光線追跡プログラムまとめ（１） Windows10アップグレード後から、VBA（Visual Basic Application）から呼び出すDLL（ダイナミックリンクライブラリ）が動作しなくなりました。
（２）対応策として、 DLL（ダイナミックリンクライブラリ）をコンソールアプリケーションに変更しました。
（３） 「近接露光実行ファイル 」フォルダーの保存場所の制限は無く、任意の場所に保存して実行できます。
（４）コンソールアプリケーションの起動はダブルクリックが必要となります。
（５）コンソールアプリケーションの終了はflag.txt（コントロールフラグファイル）で行うことができます。
（６）光線追跡の計算条件の解説は光の不思議と応用の26章〜38章を参照願います。
*光線追跡プログラムの作成は大変苦労しました。使用方法も難しく、簡単にはうまく説明できません。
*光学条件の最適化は特に難しく、多くのノウハウが必要です。
４章：照度分布計算プログラムに行く。