%  Insolate TD
% AP 15/09/2021

% trouver date ordinale JJ de votre naissance  24/10 (attention a ne pas
% utiliser le terme Julian Day)
% jj=datenum('24.10.1966','dd.mm.yyyy')-datenum('01.01.1966','dd.mm.yyyy')
jj=296;
% figure de la solar radiation du jour de votre naissance avec les heures
% en abscisse
plot((1:1440)/60,sji(jj,:))
title('Solar radiation on Oct 24')