anchor.diff:                         (subseq rate-str
                                  1 (1- (length rate-str)))
                        rate-str))))
  -             (compute-and-write-parts (rate)
  +             (write-functions (i rate)
                  (let ((rate-strs (split-rate-string-equally (rate-str-without-outer-parens rate))))                  
                    (princ #\.)
                    (loop for rate-str in rate-strs
                          for p = (incf npart)
                          collect p into parts
  -                       do (format file "~%~4Tpart~S := ~A;" p rate-str)
  -                       finally (return parts)))))
  +                       do (format file "~%~4Tpart(~S) = ~A;" p rate-str)
  +                       finally (format file "~%~4T$c(~A) = ~:[0~;~{part(~A)~^+~}~];" (1+ i) parts parts))))
  +             (scribble (obj pos)
  +               (let ((end (file-position file)))
  +                 (file-position file pos)
  +                 (prin1 obj file)
  +                 (file-position file end))))
   
       (format t "~&; Writing model block")
       (format file "MODEL ~A~%~
  
positive.code:           for i being the hash-values of kvars-table
            unless (eq kvar :last-index)
            do (setf (svref kvars (1- i)) kvar))
      (format t "~&; &Writing parameters")
      (format file "~2TWITHIN m DO")
      (with-print-context t
        (loop for k across kvars 
  
negative.code:     (format file "~A) OF REAL~%~
                    ~2TVARIABLE~%~
                      ~4Tc as ARRAY(~S) of Concentration~%~
                      ~4Tpart as Array("
              (map 'string (constantly #\space) (format nil "~A" (length allvars)))
              (length ode-vars))
  
      (setf part-pos (file-position file))
      (format file "~A) of Concentration~%~
                    ~2TEQUATION"
              (map 'string (constantly #\space) (format nil "~A" (length allvars))))
      ;; write the rate fns:
      (loop for v in ode-vars
            for i =  0 then (1+ i)
  