subroutine computing the value of
solout furnishes the solution y at the nr-th grid-point x (thereby the initial value is the first grid-point).
solout
y
nr
x