323334353637383940414243444546
sinit[0]=150.0; sinit[1]=10.0; sinit[2]=10.0; sinit[3]=0.0; Stochsim stoch = new Stochsim(sinit); for (k=0;k<NVAR;k++) smax[k]=-1; for (i=0;i<M;i++) { t=stoch.step(); for (k=0;k<NVAR;k++) if (stoch.s[k] > smax[k]) smax[k]=stoch.s[k]; if (i % 10 == 0) { j=i/10; tsav[j]=t;