1
subroutine res22(t,y,ydot,delta,ires,rpar,ipar)
3
double precision t,y,ydot,delta,rpar
6
dimension y(*), ydot(*), delta(*)
8
call f22(neq,t,y,delta)
10
delta(i) = ydot(i) - delta(i)
14
subroutine f22 (neq, t, y, ydot)
15
double precision t,y,ydot
18
dimension y(*), ydot(*)
20
ydot(2) = 100.0d0*(1.0d0 - y(1)*y(1))*y(2) - y(1)