function [neqn,nout,nsteps,t0,tf,abserr,relerr]=intpar
  neqn=2;
  nout=5;
  nsteps=1000;
  t0=0;
  tf=2.5;
  abserr=1.0e-06;
  relerr=1.0e-06;

