Problem with If statement in a Do or a While loop

*To*: mathgroup at smc.vnet.net*Subject*: [mg39437] Problem with If statement in a Do or a While loop*From*: A.Reinstaller at MERIT.unimaas.nl*Date*: Sat, 15 Feb 2003 01:42:08 -0500 (EST)*Sender*: owner-wri-mathgroup at wolfram.com

There is a problem I continue to have, perhaps one knows what's wrong. I try to check the values generated by an interpolating function. I do so by using a Do loop in which I ask whether the function has reached some value within a small intervall and if so to write the time tag of the first hit out. Now, when running the program it does not step into the IF evaluations. I have added an else branch to be sure that even if the condition was not true something would be written out. I have also checked whether the values of the interpolating function are correctly extracted as well as whether they were really lying in the intervall asked for in the IF statement. In top of that I varied the stepsize to a minimal and increased the margins of the intervall. I used also a While instead of a Do loop. But in any case neither the THEN nor the ELSE branch "fired", telling me that the IF statement is not checked. WHY???? > Is there some restriction for the use of IF clauses in loops in > Mathematica? Otherwise I could not explain why nothing happens. I checked > in the manual but nothing was mentioned. > > Thanks for your forthcoming help. Andreas > > Here is the code: > > solutionIncome= > NDSolve[Join[veblen07,IvaluesIncome],{x1,x2,y1,y2,py},{t,0,Tmax}]; > > S={(q1*y1[t]+q2*y2[t])/.solutionIncome/.t->Tmax}; > c10=S*0.01; c90=S*0.9; > t10=0;t90=0; \[CapitalDelta]t=0; > tol=0.1; stept=Tmax*0.0001; > > Do[ { > y=(q1*y1[t]+q2*y2[t])/.solutionIncome/.t->\[Theta]; > If[ ( (y>(c10-(c10*tol) )|| y<(c10+(c10*tol)) ) && t10 ==0 > ) ,( t10=\[Theta]; Print["t10 found"]; Continue[] ) , (Print["t10 not > found"]; Continue[]) ]; > If[ ( ( y>(c90-(c90*tol)) || y<(c90+(c90*tol)) ) && t90 > ==0 ) ,(t90=\[Theta]; Print["t90 found"]; Continue[]), (Print["t90 not > found"]; Continue[]) ]; > If[ (t10>0) && (t90>0),(\[CapitalDelta]t=t90-t10; > Continue[]) ]; > }, > {\[Theta],0,Tmax,stept}]; > > >