Je suis ATOUI HAMZA, voila le code source sous MATLAB 6.5 des temps de la prières
% Calcul des temps de prières
%============================
clear all,clc
% Paramètres d'entrée
%====================
% longitude
long = 7.49;
%long = 10.2167;
%long = 51.6;
%long = 31.2500;
% Altitude
Lat = 36.50;
%Lat = 36.8333;
%Lat = 25.25;
%Lat = 30.0500;
% Décalage ou la zone par rapport au grinitsh
zone = 1;
% date en Day/Month/Year
Year = 2008;
Month = 8;
Day = 22;
% La partie calcul
% Calcul du Jour Julien
D = (367*Year)-fix((7/4)*(Year+fix((Month+9)/12)))+fix(275*(Month/9))+Day-730531.5
% Longeur de Soleil
L = 280.461+0.9856474*D
if L > 360
L = (L/360 - fix(L/360))*360
end
% partie de Soleil moyenne
M = 357.528+0.9856003*D;
if M > 360
M = (M/360 - fix(M/360))*360;
end
% Longeur de Soleil toural
Lambda = L+(1.915*sin(M*pi/180))+(0.02*sin(2*M*pi/180))
% Obliquity
Obliquity = 23.439-0.0000004*D
% Alpha
Alpha = atan(cos(Obliquity*pi/180)*tan(Lambda*pi/180))*(180/pi)
if (Alpha > 360)
Alpha = (Alpha/360 - fix(Alpha/360))*360
end
if (Lambda>0 && Lambda<90)
Alpha = Alpha
elseif (Lambda>90 && Lambda<180)
Alpha = Alpha+180
elseif (Lambda>180 && Lambda<360)
Alpha = Alpha+360
end
% if (Alpha > 360)
% Alpha = (Alpha/360 - fix(Alpha/360))*360
% end
Alpha = Alpha+(90*(fix(Lambda/90)-fix(Alpha/90)))
% Temps de Soleil
ST = 100.46+0.985647352*D
if (ST > 360)
ST = (ST/360 - fix(ST/360))*360
end
% Declinaison
Dec = asin(sin(Obliquity*pi/180)*sin(Lambda*pi/180))*(180/pi);
%Dec = abs(Dec)
% Noon
Noon = Alpha-ST
if Noon > 360
Noon = (Noon/360 - fix(Noon/360))*360
end
if Noon < 0
Noon = Noon+360
end
%UTnoon
UTnoon = Noon-long
% Douhr
Localnoon = UTnoon/15+zone
% Altitude du centre de soleil "Shafii"
AsrAlt = atan(1+tan(abs(Lat-Dec)*pi/180))*(180/pi)
% Altitude du centre de soleil "Hanafi"
%AsrAlt = atan(2+tan(abs(Lat-Dec)*pi/180))*(180/pi);
% Arc de Asr
AsrArc = acos((sin((90-AsrAlt)*pi/180)-sin(Dec*pi/180)*sin(Lat*pi/180))/(cos(Dec*pi/180)*cos(Lat*pi/180)))*(180/pi)
% Asr
AsrTime = Localnoon+(AsrArc/15)
% Durinal Arc
DurinalArc = acos((sin(-0.8333*pi/180)-sin(Dec*pi/180)*sin(Lat*pi/180))/(cos(Dec*pi/180)*cos(Lat*pi/180)))*(180/pi)
% Sun Set
SunSet = Localnoon+(DurinalArc/15)
% Isha Arc
IshaArc = acos((sin(-18*pi/180)-sin(Dec*pi/180)*sin(Lat*pi/180))/(cos(Dec*pi/180)*cos(Lat*pi/180)))*(180/pi)
% Isha Time
IshaTime = Localnoon+(IshaArc/15)
% Fajr Arc
FajrArc = IshaArc
% Fajr Time
FajrTime = Localnoon-(FajrArc/15)
% Affichage des temps de prières
HR = fix(FajrTime);
MN = fix((FajrTime-fix(FajrTime))*60);
msg1 = sprintf('Fajr Time %2.2d:%2.2d',HR,MN);
disp(msg1);
HR = fix(Localnoon);
MN = fix((Localnoon-fix(Localnoon))*60);
msg2 = sprintf('Douhr Time %2.2d:%2.2d',HR,MN);
disp(msg2);
HR = fix(AsrTime);
MN = fix((AsrTime-fix(AsrTime))*60);
msg3 = sprintf('Asr Time %2.2d:%2.2d',HR,MN);
disp(msg3);
HR = fix(SunSet);
MN = fix((SunSet-fix(SunSet))*60);
msg4 = sprintf('Magrib Time %2.2d:%2.2d',HR,MN);
disp(msg4);
HR = fix(IshaTime);
MN = fix((IshaTime-fix(IshaTime))*60);
msg5 = sprintf('Isha Time %2.2d:%2.2d',HR,MN);
disp(msg5);
Priere = strvcat(msg1,msg2,msg3,msg4,msg5);
msgbox(Priere);
