referring to fig1 a method and system for evaluating the effectiveness of a person to perform a specific task according to the present invention is illustrated . a user can employ the system to predict modulation of performance based on regulation of a sleep pattern . the factors utilized to regulate sleep are broadly indicated in the sleep regulation section 40 , shown in fig1 . these factors include sleep intensity 42 , sleep debt 44 , sleep accumulation 46 , and sleep quality or fragmentation 48 . performance modulation section 50 contains factors that are used for such performance modulation , including circadian rhythm 64 , sleep inertia 52 , effectiveness 56 and performance use 54 , as shown in fig1 . a circadian process 60 influences both the performance modulation 50 and sleep regulation 40 . sleep regulation 40 is dependent on the hours of sleep and wakefulness , current sleep debt , the circadian process , which is represented by the circadian oscillators 60 , and fragmentation 48 , namely awakenings during periods of sleep . performance modulation 50 depends on the current balance of the sleep reservoir , the circadian process , and sleep inertia . through the implementation of mathematical modeling the system and method of the present invention can predict changes in cognitive performance . one embodiment of the present invention makes use of a multi - oscillator circadian process , a circadian sleep propensity process , a sleep fragmentation process , and a circadian phase adjusting feature 66 for time zone changes . these processes can be implemented on a general - purpose digital computer . performance while awake and the drive to sleep are both controlled , in part , by a circadian process 60 . numerous studies of performance , reaction time , alertness ratings , measures of the tendency to fall asleep , and body temperature indicate that the underlying circadian process is not a simple sine wave . performance and alertness reach a major peak in the early evening , about 2000 hours , and fall to a minimum about 0400 hours . there is a secondary minimum in the early afternoon , about 1400 hours , and a morning peak at about 1000 hours . correlated with this pattern is a rising tendency to fall asleep that reaches a peak at about the same time performance and alertness reach a minimum . the existence of both a major and a minor peak in performance and two corresponding minima at other times , suggest that at least two oscillators are involved in the circadian process . these multiple oscillators are accounted for in the method for evaluating the effectiveness of the present invention . one embodiment of the present invention incorporates a circadian process that is composed of the sum of two cosine waves , one with a period of 24 hours , and one with a period of 12 hours . the two sinusoidal oscillators are out of phase thereby producing a predicted variation in arousal that closely parallels known patterns of body temperature . referring to fig2 a graph of the complex circadian oscillator is illustrated using units of percent change in effectiveness in one variation of the invention . the circadian process within the method for evaluating the effectiveness generates an arousal function that mirrors circadian changes in oral temperature . this arousal oscillator drives both variations in predicted cognitive effectiveness and sleep propensity . these two translations of the oscillator have identical frequency and phase components and differ only in amplitude and sign ; a rise in arousal produces an increase in performance and a decrease in propensity to sleep . the circadian process 60 is depicted in fig1 . the control of sleep and the influence of sleep on cognitive capacity is a homeostatic process . at the center of the homeostatic process is the sleep reservoir 70 , shown in fig1 . a fully rested person has optimal performance capacity , indicated as the reservoir capacity r c . while awake , the reservoir is depleted according to a performance use function , indicated by the arrow 54 pointing away from the reservoir 70 . while asleep , the reservoir is filled to replenish the capacity to perform and be alert as indicated by the arrow 46 representing sleep accumulation . the rate of accumulation for each minute of sleep is called sleep intensity , and is driven by two factors : 1 ) the circadian variation in sleep propensity 62 , and 2 ) the current sleep deficit 44 , which is the reservoir capacity r c minus the current level of the reservoir r t at time t . this deficit is constantly changing as a person sleeps and replenishes the reservoir 70 , or is awake and depletes the reservoir 70 . the effectiveness 56 is determined from the circadian oscillators 60 and performance rhythm 64 , the sleep reservoir 70 , and the sleep inertia 52 . the circadian oscillators 60 are represented by the following equation : t = time of day , p = 24 hr phase in hours , p ′= 12 hr relative phase in hrs , and β = rel . amplitude of 12 hr cycle . m = mesor ( a point around which a sine wave oscillates ), and a s = sleep propensity amplitude . where a p = a 1 + a 2 ( r c − r t )/ r c , performance rhythm amplitude . the connection of the sleep reservoir to the sleep debt calculation 44 forms a feedback loop with the sleep intensity 42 , the sleep accumulation 46 , and back to the sleep reservoir 70 . the sleep debt is represented by the following equation : r c = reservoir capacity , r t = current reservoir balance , ( r c − r t )= current reservoir deficit , and f = amplitude of feedback . the reservoir balance of the sleep reservoir 70 is represented by the following function : r t = r t − 1 + s − p , and represents the total sleep units at time interval t . the sleep intensity 42 is represented by the following equation : the sleep inertia 52 is represented by the following equation : where max = i max , and i = inertia time constant for two hr after awakening . sleep fragmentation 48 caused by a poor quality sleep environment results in a pause in sleep accumulation after returning to sleep following each awakening from a sleep interval . the performance use 54 can be represented by the following linear equation : the prediction of the model is cognitive effectiveness 56 interpreted as percent of baseline cognitive throughput and calculated by the equation : table 1 shows the default values for the variables in the foregoing equations . a method according to the present invention is also illustrated in the flow chart shown in fig1 . in step 200 the sleep schedule ( past and / or future ) is input to the simulation . at each time epoch , t , of the schedule , the sleep state ( asleep or awake ) is determined and the simulation is updated accordingly . in step 202 the circadian oscillators are modeled . when asleep , the circadian oscillators are used to calculate the sleep propensity 204 , which partly determines sleep intensity in step 206 . after any pause caused by fragmentation 208 , sleep accumulation is calculated in step 210 , which is used to update the amount of effective sleep in the sleep reservoir 70 in step 218 . the sleep debt is determined in step 222 based on the sleep reservoir calculated in step 218 and is used in step 206 to calculate the sleep intensity . when awake , performance use is calculated in step 216 and also triggers a sleep fragmentation determination in step 208 . at the start of each awake period , sleep inertia is calculated in step 220 and contributes to the cognitive effectiveness calculation at 230 . the sleep inertia from step 220 , the amount of effective sleep in the reservoir 70 , calculated in step 218 , and the circadian oscillators , modeled in step 202 are used to calculate the performance rhythm at step 214 , are all used to calculate the effectiveness to perform the task in step 230 . steps 218 , 222 , 206 , 208 and 210 form a feedback loop . fig2 is a graph representing circadian oscillation in performance corresponding to the temperature and arousal rhythm as a function of the time of day . the arousal oscillator 60 drives both variations in predicted performance ( cognitive effectiveness ) and sleep propensity . these two translations of the oscillator have identical frequency and phase components and differ only in amplitude and sign . a rise in arousal produces an increase in performance and a decrease in propensity to sleep . fig2 illustrates an inventive simulation of the performance rhythm for the first day of a fully rested person and after three days of sleep deprivation . the amplitude of the performance rhythm increases with the accumulated sleep debt on the third day . as can be seen in fig3 the circadian rhythm , as represented by oral temperature , and the subjective alertness are closely related , but not identical since alertness and performance are also sensitive to level of the sleep reservoir . the oscillation in the reservoir level is called the sleep - wake cycle and represents the current sleep debt ( r c − r t ). sleep accumulation does not start immediately upon retiring to sleep . following a period of wakefulness there is a minimal delay of approximately 5 minutes required to achieve a restful sleep state . this factor accounts for the penalty during recuperation that is caused by sleep in an environment that leads to frequent interruptions . these components of the sleep accumulation function are indicated as 42 and 48 , respectively , in fig1 to the left of the sleep reservoir 70 feeding into the sleep accumulation function 46 . the level of the reservoir 70 at time t + 1 is the level at time t , r t , plus sleep accumulation ( s ) while asleep and minus performance use ( p ) while awake . the units of the reservoir 70 are minutes of effective sleep . the method of the present invention can easily accommodate a complex pattern of sleep and waking . while asleep , the simulation adds to the reservoir 70 ; while awake the simulation depletes the reservoir 70 . a schedule can oscillate between these states as often as once a minute and the simulation will keep account of the net effects on performance capacity as the balance in the reservoir 70 , similar to the balance in a check book . the total reservoir capacity is 2880 units ( nominal minutes ) of effective sleep . this value is based on the following considerations . the average person is assumed to require 8 hours of sleep and is awake for 16 hours in the typical day . to remain in balance , then , sleep units must accumulate at twice the rate it is used during performance . hence , the rate of performance use , κ , is 0 . 5 relative to the average rate of sleep accumulation . studies of total sleep deprivation indicate that cognitive capacity depletes at a rate of about 25 % per day . hence , the reservoir has the capacity to sustain performance for four days . this translates into 2880 units of sleep : 4 days × 24 hours × 60 min / hr × 0 . 5 . the outcome of the reservoir process , according to the method of the present invention , during continuous sleep converges to an exponential accumulation function , if one ignores the circadian effects on sleep intensity . one embodiment of the method of the present invention is based on minute - by - minute additions to the reservoir 70 during sleep , with the size of these increments proportional to the reservoir deficit ( the feedback process ). integrated over time , this iterative process is described by an approximate exponential function , but is not an exponential function ; rather , it is a moment - by - moment simulation of the effects of sleep on the reservoir balance . therefore , the method can easily accommodate a momentary interruption in sleep ( fragmentation ) caused by a poor quality sleep environment . the incremental process is interrupted for the duration of the awakening and the reservoir 70 is depleted for that period of time by the performance function . upon return to sleep after an interruption , there is a delay before resumption of sleep accumulation , with the delay set at 5 min in one embodiment of the method . the result of this process is directly tied to real world events that drive the process , not to an a priori mathematical equation . the feedback process of the method of the present invention is used to determine the effects of long schedules of less than optimal sleep . such schedules deplete the reservoir 70 , and increase the intensity of sleep when sleep occurs . eventually , the greater average intensity of sleep permits the person to adjust to such a schedule and find a new equilibrium of sleep and stable performance , within limits . performance will not be as effective as it might be with a full eight hours of sleep , but performance does not necessarily degrade indefinitely . this is what is meant by a homeostatic sleep and performance process . it is much like a person adjusting to a restricted diet ; the person loses weight and conserves energy so that a new equilibrium stable weight is reached under the limited input of calories . according to the method of the present invention , cognitive effectiveness and alertness are primarily dependent on variations in the two processes just described : the endogenous circadian rhythm ( reflected in oral temperature ) and current sleep reservoir balance resulting from the sleep - wake cycle . a third factor is a temporary disturbance in performance that may occur immediately following awakening , called sleep inertia i . the variations in measured alertness , shown in fig3 as a dashed line , do not exactly match the variations in oral temperature . this is because , according to one aspect of the present invention , simulated cognitive effectiveness is computed as the sum of these three factors : the relative level of the sleep reservoir in percent units ( 100 × rt / rc ) plus the effects of the circadian oscillator , c , and minus sleep inertia , i . sleep inertia is represented as an exponential decay function of sleep intensity at the time of awakening and lasts for at most 2 hours . the predictions are normally in terms of changes in cognitive effectiveness , expressed as percent of baseline performance for a person when well rested . this measure corresponds to performance of a standard serial add - subtract task or the average of a range of standard cognitive tests , described below in greater detail . in addition , the parameters of the performance calculation can be adjusted to predict other components of performance , such as reaction time , lapses in attention , and target error . the method of the present invention can be used to make a number of predictions , such as for example , performance and alertness . the average person is assumed to require eight hours of sleep per day to be fully effective and to avoid accumulation of sleep debt . based on the joint interaction of the endogenous circadian oscillator , c , and the sleep - wake cycle , performance is predicted to have two peaks in percent effectiveness at approximately 1000 hours and 2000 hours , a minor dip in performance at about 1400 hours , and a major trough in effectiveness during the early morning hours when the person is normally asleep . fig4 displays the predicted variations in performance for a person sleeping eight hours per day from midnight to 0800 hours , heavy line . since the individual is asleep between midnight and 0800 , predicted performance is shown as a fine line indicating potential effectiveness , if awake . the nighttime pattern reveals a major trough in performance at about 0300 hours and a minor trough in performance at about 1400 hours . referring to fig3 the empirical pattern of alertness closely parallels the prediction of the method of the present invention with two peaks in alertness , a mid - afternoon dip in alertness , and a major trough in alertness at 0600 hours . note that what is important is the pattern of performance and alertness , not the exact time of peaks and troughs , which is a parameter of the model . a number of studies agree with the bimodal pattern of performance shown in fig3 and 4 . lavie ( 1991 ) reported the results of a study of sleep related traffic accidents in israel between 1984 and 1989 . in fig5 the connected stars reveal two peaks in sleep related accidents , a major peak at about 0300 hours and a minor peak at about 1500 hours in the afternoon . these correspond to the dips in performance predicted by the model in fig4 . fig6 shows the results from a variety of studies of performance from industrial settings . all performances are scaled so that “ good ” performance was high on the y - axis and “ bad ” performance was low on the y - axis . generally , there were two dips in performance , one at about 0300 hours and a second at about 1400 hours . all these results are consistent with the predictions of the model shown in fig4 . the basis for the performance pattern shown in fig4 is illustrated in fig7 . within the model , performance is the sum of two major factors , the oscillation of the sleep reservoir balance , top curve , and the oscillation of the circadian arousal process , bottom curve . the two curves are roughly in phase so that as sleep debt is accumulated while awake ( top curve ), the circadian arousal process is increasing and largely offsets the change in performance . performance is not perfectly constant , however , for two reasons . first , the two processes are not exactly in phase and this causes a strong early morning trough when the two cycles are both decreasing . second , the temperature rhythm displays an afternoon plateau and is responsible for the mid - afternoon dip in performance when sleep debt accumulates while temperature is roughly constant . as described earlier , the circadian process produces an oscillation in sleep propensity , shown as the solid line in fig8 . this rhythm is the negative of the arousal rhythm shown in fig2 and scaled in sleep units . sleep propensity sp combines with the current sleep debt sd resulting from the sleep - wake cycle , shown as the dashed line in fig8 . the summation of this sleep - wake cycle and the circadian process is called sleep intensity si and is diagrammed in fig9 . for a person taking a normal 8 hours sleep from midnight to 0800 hours , sleep is most intense in the early morning at about 0300 hours . there is a mid - afternoon increase in sleep propensity at about 1500 hours that coincides with the mid - afternoon dip in alertness , fig3 and increases in sleep related traffic accidents , fig5 . the method for evaluating effectiveness according to the present invention can be used to predict equilibrium states . a homeostatic representation of sleep regulation leads to an important implication that is seldom recognized . if a subject is scheduled to take less than an optimal amount of sleep each night , for example , four hours per day , the reservoir 70 initially loses more units during the awake period than are made up during the sleep period . this results in a sleep debt at the end of the sleep period that accumulates over days . however , since the rate of sleep accumulation increases with sleep debt , eventually , the rate of sleep accumulation increases such that four hours of sleep makes up for twenty hours awake . at this point , the reservoir 70 reaches an equilibrium state and no further debt is accumulated , although the initial deficit remains as long as the person remains on this schedule . the result of this process is shown in fig1 , in which a schedule of one day with 8 hours of sleep followed by six days of 4 hours of sleep per day is illustrated , with each sleep period starting at midnight . by the sixth day of the restricted sleep schedule , cognitive performance oscillates about a stable level well below the baseline level achieved with 8 hours of sleep . minimum effectiveness is about 60 % on the seventh day . fig1 also illustrates the operation of sleep inertia si during the first two hours following awakening according to the simulation . with only four hours sleep per night , the average intensity of sleep is high compared to a normal 8 hours sleep period . as a consequence , sleep inertia , which is driven by sleep intensity at time of awakening , is relatively high . this is evident as the initial “ notch ” in performance immediately after awakening in this seven - day record of performance . progressive sleep debt under extreme schedules can also be predicted using the method for evaluating the effectiveness according to the present invention . the sleep homeostat is not infinitely elastic ; there is a limit to the rate of sleep accumulation ( sleep intensity ). any schedule that provides less than 3 hours of sleep per day ( for the average person ) will not reach an equilibrium state and performance capacity will gradually deplete to zero , although the rate of depletion slows over the first week of restriction as sleep intensity rises to its maximum level . fig1 is a graph illustrating a simulation of an extended record of sleep and performance under a schedule of only 2 hours of sleep per day . under such a schedule , minimum performance declines to about 17 % at the end of the seventh day . the performance effects of sleep timing can be predicted . the method for evaluating the effectiveness of the present invention is sensitive to the time of day of the sleep period . fig1 illustrates the predicted performance of an individual given eight hours of sleep per day , starting at 1200 hours ( noon ) each day , 12 hours out of phase from that shown in fig4 . although performance reaches a peak of 100 % at the start of each awake period ( 2000 hours ), predicted performance then rapidly declines during the late night and early morning hours to a strong dip at about 0500 hours . minimum predicted performance under this schedule is predicted to be as low as 69 % compared to minimum performance under a normal sleep schedule of 86 %. note that for this chart , the phase correction algorithm was disabled so that the fall effects of out of phase sleep could be visualized . this alteration in pattern results from two factors . first , sleep intensity is initially less for sleep periods starting at noon . this results in accumulation of a small debt that is quickly offset by the homeostatic sleep mechanism . the second , more persistent pattern is the combined effect of the circadian oscillation of performance that reaches its minimum in the early morning hours and the sixteen hours of accumulated sleep debt at the end of the awake period that predicts a strong dip in effectiveness toward the end of the awake period . this pattern has strong implications for performance under shift schedules that require daytime sleep . it is well documented that most mistakes on the night shift occur during the early morning hours and this outcome is predicted by the present method for evaluating effectiveness . the system and method for evaluating effectiveness of the present invention can predict changes in cognitive capacity as measured by standard laboratory tests of cognitive performance . it is assumed that these tests measure changes in the fundamental capacity to perform a variety of tasks that rely , more or less , on the cognitive skills of discrimination , reaction time , mental processing , reasoning , and language comprehension and production . however , specific tasks , such as specific military tasks vary in their reliance on these skills , and deficits in cognitive capacity may not produce identical reductions in the capacity to perform all military tasks . it is reasonable to assume , however , that the changes in military task performance would correlate with changes in the underlying cognitive capacity . in other words , if one were to plot changes in military task performance as a function of measured changes in cognitive capacity , there would be a monotonic relationship between the two variables . therefore , if these two sets of data were available from a test population subjected to sleep deprivation , linear ( or non - linear ) regression techniques could be applied to derive a transform function ; this transform translates predicted cognitive changes into changes in military task performance . based on this reasoning , the method for evaluating the effectiveness , discussed previously as the cognitive effectiveness 56 , can be extended to predict variations in any task or component of a task ( given appropriate test data ) using the generalized task effectiveness ( te ) expression as follows : te = a ( r t / r c )+ b + c 1 [ cos ( 2π ( t − p )/ 24 )+ c 2 ( cos ( 4π ( t − p − p ′)/ 24 ))]+ i , a = linear component slope , b = linear component intercept , c 1 = circadian weighting factor , c 2 = 12 hr weighting factor , and p = acrophase of the task . with regard to “ jet lag ” the system and method for evaluating effectiveness can predict the adaptation of performance to changes in time zones that might accompany trans - meridian flights or that might occur if the subject shifts to a regular schedule of nighttime work . fig1 is a graph illustrating the adjustment of performance to two flights , an east - bound flight across six time zones and a west bound flight across 6 time zones . predicted performance while awake is more disrupted for a longer period of time by the east bound flight compared to the west bound flight , a commonly reported difference in “ jet lag ” for east and west plane travel . these effects on performance are a logical and inherent outcome of the interplay of the various processes in the model and do not require a special “ jet lag ” algorithm . the method for evaluating the effectiveness of the present invention has logic to detect such a change in work pattern and to readjust the phase of the circadian rhythm to the new work pattern indicative of the shift in time zone . the performance prediction was a natural outgrowth of the shifting circadian phase . when a person moves to another time zone or alters work patterns so that sleep and work occur at different times of day , the internal circadian oscillator that controls body temperature and alertness shifts to this new schedule . during the period of adjustment , a person experiences performance degradation , disrupted mood , and feelings of dysphoria , called circadian desynchronization or “ jet lag ”. the method for evaluating the effectiveness of the present invention can simulate this process and automatically adjust the phase of the circadian rhythm to coincide with the activity pattern of the person , at step 212 , fig1 . such simulation is important in order to accurately predict the effects of moving to a new time zone or changing to a new and regular work pattern , such as changing from the day shift to the night shift . the method for evaluating the effectiveness detects the average time of the period of wakefulness , and gradually adjusts the time of the peak in the circadian rhythm in relation to the new average awake hour . the details of the phase adjusting process 212 are charted in fig1 . the peak of the circadian rhythm has a reliable relationship to the timing of the period of wakefulness . when a person moves to a new work schedule or a new time zone , the change in average awake time 301 and 302 ( relative to a reference time zone ) is detected and a new “ goal phase ” is computed , 304 . for example , when moving from the central us time zone to germany , the awake time of the subject advances six hours . instead of waking at , say , 0600 central time , the subject awakens at 0000 central time , which is 0600 german time . this causes a shift of 6 hours in the “ goal phase ” of the person . however , the human physiological system does not adapt immediately to such a shift . in general , the speed of adjustment depends on whether the change requires a phase advance or a phase delay . a phase advance ( eastward time change ) takes about 1 . 5 days per hour of shift . one embodiment of the present invention adjusts to the new “ goal phase ” gradually over the course of about nine days for a six hour advance , steps 306 , 308 , 310 , 318 , 319 , 320 , 322 . during that time , the performance of the subject will show degradation due to the desynchronization of the internal circadian rhythm from the new rhythm of work and sleep . likewise , westerly travel causes a phase delay in the circadian rhythm and research shows that phase delays take less time for adjustment , about one day per hour of shift , or six days for a six hour time change , steps 306 , 308 , 312 , 313 , 314 , 316 . a basic premise of the method is that periods of wakefulness provide a strong entrainment stimulus that serves as the basis of phase adjustments to the circadian rhythm . as a result , the time of the period of wakefulness is the primary factor in determining the phase of the circadian rhythm . at the end of each period of wakefulness the model computes a new “ running average awake hour ,” which is the average of that period of wakefulness averaged with the average awake hour of the prior two periods of wakefulness , 301 and 302 . although the wakefulness period is used , this is not to imply that there are not other entrainment stimuli for adjusting the phase of the circadian rhythm such as bright light , activity , etc . also , brief awakenings during the night lasting only a few minutes ( for example , less than 60 minutes ) are disregarded . based on the running average awake hour , the system computes a new “ goal phase ” as the running average awake hour plus , for example , 3 hours , 304 . for an awake time from 08 to 24 ( average of 16 ), that gives an acrophase of 19 , which is a current default . the system adjusts the current acrophase each minute by comparing the goal to the prior current acrophase . if the difference is greater than 1 hr (+/− 1 ), it starts adjusting the acrophase by { fraction ( 1 / 1440 )} of an hour for each minute of delay ( 314 ) or { fraction ( 1 / 2160 )} of an hour for each minute of advance ( 320 ). the system adds or subtracts , depending on whether the difference is + or −. if the difference is less than an hour , it adjusts by that difference times { fraction ( 1 / 1440 )} or { fraction ( 1 / 2160 )} each minute , 316 or 322 . it checks to see if the difference is greater than 12 hours , in which case it is shorter to go in the opposite direction . in that case , the current acrophase is adjusted by + or − 24 , and the direction of change is reversed . this method ensures that the current acrophase is adjusted along the shortest path to the goal phase . the result is a model that adjusts about one hour a day per hour of phase delay and one hour per 1 . 5 days for each hour of phase advance . it adjusts smoothly and continuously so that it will also deal with rotating shifts , in which case the goal phase shifts back and forth and the adjusting acrophase keeps trying to find the goal at the prescribed rate , but never gets there unless the schedule stabilizes . other factors , such as bright light or dietary factors may be incorporated as factors that change the rate of adjustment of the physiological rhythm to the predominant schedule represented by the goal phase . as used herein a number of terms have specific meanings . in particular , the term average awake hour 301 refers to the average awake hour of any awake period based on the average clock time of the awake period , provided the awake period is one hour or longer . therefore , a period from 0800 hrs to midnight would have an average clock time of 1600 hrs ( the middle of the period ). the running average awake hour 302 is the average of the average awake hours for the last three periods of wakefulness . the three periods are weighted so that the period just ended has greater weight than either of the two previous periods and the previous period greater weight than the period prior to it , for example ( 1x + 0 . 67y + 0 . 33z )/ 2 . the term goal phase is simply the running average awake hour plus a constant displacement , which has a default value of + 3 hours , 304 . this gives a ‘ standard ’ acrophase of 1900 hrs for a 16 hour period of wakefulness from 0800 hrs to midnight . the term adjusting acrophase ( steps 306 to 322 ) refers to a calculation that adjusts at a rate of one hour per day for phase delays or ⅔ hour per day for phase advance until it is within 1 hour of the goal , and then gradually approaches it from that point to zero deviation . the adjusting acrophase is needed to account for the fact that the biological phase cannot shift in one big jump to the goal phase . the following is a breakdown of the determination of the running average awake hour according to the present invention . ah n = average awake hour for awake period n ( 301 ). at the end of each awake period calculate new av , where n = 0 is just ended awake period , n = 1 is previous awake period , and n = 2 is awake period prior to that . equation 1a , ah n , average awake hour for awake period n : ah n = if ( tae −( ca / 60 )& lt ; 0 ), then ( tae + 24 +( tae + 24 −( ca / 60 )))/ 2 ), av =( 0 . 33 * ah 2 + 0 . 67 * ah 1 + w * ah 0 )/( w + 1 ), ( tae −( ca / 60 )) is the time of the start of the awake period in hours , and ( tae +( tae −( ca / 60 )))/ 2 calculates the average awake hour 301 as the simple average of the beginning time and ending time . for unusually long awake periods greater than 24 hours , the average hour of the last portion of the awake period less than 24 hrs is averaged with a time 12 hrs after the start of the awake period for each proceeding 24 hour period awake . av 302 is calculated as the weighted average of the last three awake periods , such that the just ended period has weight equal to half the duration in hours of the just ended awake period , the two prior periods have a combined weight of one , and the prior period has a weight twice that of the period prior to it . given : ra = a parameter that sets the relative acrophase based on the average awake hour . for example an ra = 3 gives an acrophase of 19 by adding it to the average awake hour of 16 for an awake period from 0800 to 2400 . this is a parameter . then : at the end of each awake period calculate a new goal phase : the acrophase of the arousal / temperature rhythm is set relative to the average awake hour . that is the essence of using the awake time as the entrainment stimulus for the phase shifting model . as average awake hour shifts with a new schedule or time zone , then the goal phase adjusts also . the actual phase of the temperature rhythm does not adjust immediately to this new goal phase , but gradually moves toward it according to the adjustment algorithm discussed next . pc a = amount of phase change in minutes for each hour of advance , default value is 2160 minutes per hour of advance ( eastward flight direction ), pc d = amount of phase change in minutes for each hour of delay , default value is 1440 minutes per hour of delay ( westward flight direction ). then : at each minute of schedule , adjust the current acrophase according to the following algorithm : if ga − cp & gt ; 12 , then new cp = old cp + 24 − 1 / pc a else , new cp = old cp +( ga − cp )* 1 / pc d , 316 if cp − ga & gt ; 12 , then new cp = old cp − 24 + 1 / pc d else , new cp = old cp +( ga − cp )* 1 / pc a , 322 the initial if statements 308 and 310 check to see which is greater , goal phase or current acrophase . that determines which way to move the current acrophase , i . e . phase delay or advance . the next level of if statements determines if the difference is greater than one hour , the basic unit of change , 312 and 318 . if it is greater than one , then the acrophase is moved by +/− 1 / pc hours ; if it is less than one , the acrophase is moved by the difference times 1 / pc . this gives an average rate of change of about 1 / pc and then approaches the goal value gradually . the next level of if statements determines if the difference is greater than 12 hours , 313 and 319 . if it is greater than 12 hours , then 24 is subtracted or added as appropriate and moved in the opposite direction . otherwise , cp is adjusted by the delay or advance rate as appropriate , 314 or 320 . this determines the shortest distance and adjusts the direction of change accordingly . note that the method of change when the difference is less than 1 hour is the same expression , regardless of which is greater , ga or cp . this is because the ( ga − cp ) term will have the appropriate sign to move in the correct direction . fig1 illustrates a block diagram of a system for predicting task effectiveness according to present invention . data is entered using input devices 400 . these input devices 400 may include a keyboard and mouse , a data storage medium from which the data maybe extracted , and sensors that provide time and sleep / wakefulness data . the system 402 may include a microprocessor that uses the data entered through the input devices 400 to provide results either to the display 404 or the print out 406 . these results may be similar to the graphs shown in fig4 , 11 , and 12 . the system can utilize actual data for a particular person to provide task effectiveness data that is particular to that person , or can be used to predict task effectiveness for the general population or a sub - group of the population . for example , fig1 and 11 are graphs illustrating sleep and predictions of cognitive effectiveness with four and two hours of sleep per night , respectively , while fig1 shows a graph of the predicted performance of an individual given eight hours of sleep per day , starting at 1200 hours ( noon ) each day . these figures illustrate the predictive capacity of the system of the present invention . the raw sleep / wake data is entered , most likely , through a keyboard so that the system can produce the general predictive results , such as shown in fig1 - 12 . for a particular person , the actual sleep data may be entered for that person . the data may be entered to the system 402 via a keyboard and / or through sensors that provide the sleep / wake data for that person . similarly , the data concerning the cognitive skills required for a particular task may be entered using the keyboard or may be downloaded from a storage medium . having described several embodiments of the system and method for predicting task effectiveness in accordance with the present invention , it is believed that other modifications , variations and changes will be suggested to those skilled in the art in view of the description set forth above . it is therefore to be understood that all such variations , modifications and changes are believed to fall within the scope of the invention as defined in the appended claims .