/*calculate the "t"*/
                        t = (2 * pi * ((day - 1) / 365.0));
                        
                        /*calculte the fDeclination*/
 Declination = (0.322003- 22.971 * cos(t)- 0.357898 * cos(2*t)- 0.14398 * cos(3*t)+ 3.94638 * sin(t)+ 0.019334 * sin(2*t)+ 0.05928 * sin(3*t));

                        fDeclination=Declination*pi/180;
                        
                        /*calculate the altitude*/
fAltitude = asin(sin(fDeclination) * sin(latitude) + cos(fDeclination) * cos(latitude) * cos((15 * (time - 12)) * (pi / 180)));
Altitude =fAltitude*180/pi;

                           /*calculate the azimuth*/
fAzimuth = (acos((cos(latitude) * sin(fDeclination) - cos(fDeclination) * sin(latitude) * cos((15 * (time - 12)) * (pi / 180))) / cos(fAltitude)));

                      /* Adjust, to keep the range valid */
					if (time > 12) {
			Azimuth =-180+fAzimuth*180/pi;
						               } 
                        else {
			Azimuth =180-fAzimuth*180/pi;
			     }
                        
                        
                  /* Convert back to steps, from radians. */
fazimuth = (Azimuth / 1.5 * 180.0 / pi); // 1.5 IS THE STEPS PER DEGREE MOVEMENT//
faltitude = (Altitude / 1.8 * 180.0 / pi); // 1.8 IS THE STEPS PER DEGREE MOVEMENT//


  // FOR THE MOVEMENT OF THE RESPECTIVE MOTORS AS PER THE CALCULATED ANGLES//

	for (i= 0; i =fazimuth; i++) {
		movealtitude(CLOCKWISE, ALTITUDEONDELAY);
	                                }

	for (i = 0; i =faltitude; i++) {
		moveazimuth(CLOCKWISE, AZIMUTHONDELAY);
	                               }