setwd("c:/R/bme_pde/chap13") source("ge.R") tout=seq(from=0,to=20,by=2); nout=11;gt=rep(0,nout);ncase=2; for(it in 1:nout){ gt[it]=ge(tout[it]); cat(sprintf("\n t = %5.1f g(t) = %8.5f",tout[it],gt[it])); }