How to compute and map day-night terminator using MATLAB?

Hello all,

Could someone offer suggestions for matlab code to compute day/night terminator position?

The projections will be Mercator and polar stereographic; presumably, it will suffice to have Mercator-based terminator map only.

Thanks for any suggestions.

-- Liming

[Kristina Collins]

https://xkcd.com/979/

[Kristina Collins]

https://www.mathworks.com/matlabcentral/mlc-downloads/downloads/submissions/54875/versions/9/previews/Geostationary%20Satellites%20Tracking%20-%20Tony%20C%20Huang%20V1.11/private/plotdaynightterminator.m/index.html

[Denis Anikiev]

I came up with the same request, I believe it is from this project:

https://www.mathworks.com/matlabcentral/fileexchange/54875-geostationary-satellites-tracking

https://huangproject.wordpress.com/