m = input("lancer la saisie");
a = input(" Rentrer l''année d''origine : ");
b = input(" Rentrer la dernière année : ");
n = input(" Nb de relevés : "); 
w = input(" Maximum de la moyenne annuelle des températures :");
u = input(" Maximum de l''amplitude saisonnière des températures :");
r = input(" Période de variation de la moyenne annuelle des températures :");
s = input(" Période de variation de l''amplitude saisonnière des températures :");

t=linspace(a,b,n+1);
plot2d(t',(((u*sin((2*%pi*t)/r)).*sin(2*%pi*t))+((w*sin(((2*%pi)/s)*t))))');
xtitle('milankovitch ','temps en années','température');
xgrid();