CMathODE *SpanNew()
{
  CMathODE *ode=ODENew(0,1e-4,ODEF_HAS_MASSES|ODEF_PAUSED);
  ode->derivative=&MyDerivative;
  ode->drag_v2=0.002;
  ode->drag_v3=0.00001;
  ode->acceleration_limit=5e3;
  QueIns(ode,Fs->last_ode);
  return ode;
}

U0 SpanDel(CMathODE *ode)
{
  if (ode) {
    QueRem(ode);
    QueDel(&ode->next_mass,TRUE);
    QueDel(&ode->next_spring,TRUE);
    ODEDel(ode);
  }
}

#define M_SIZE (sizeof(MyMass)	-offset(CMass.start_saved_area))
#define S_SIZE (sizeof(MySpring)-offset(CSpring.start_saved_area))

U8 *SpanSave(CMathODE *ode,I64 *_size=NULL)
{
  I64 cnt;
  U8 *result,*ptr;
  MyMass   *tempm;
  MySpring *temps;
  SpanHeader h;

  ODERenum(ode);
  h.version=SPAN_VERSION;
  if (ode->next_mass!=&ode->next_mass)
    h.num_masses=ode->last_mass->num+1;
  else
    h.num_masses=0;
  if (ode->next_spring!=&ode->next_spring)
    h.num_springs=ode->last_spring->num+1;
  else
    h.num_springs=0;

  cnt=sizeof(SpanHeader)+h.num_masses*M_SIZE+h.num_springs*S_SIZE;

  ptr=result=MAlloc(cnt);
  MemCpy(ptr,&h,sizeof(SpanHeader));
  ptr+=sizeof(SpanHeader);

  tempm=ode->next_mass;
  while (tempm!=&ode->next_mass) {
    MemCpy(ptr,&tempm->start_saved_area,M_SIZE);
    ptr+=M_SIZE;
    tempm=tempm->next;
  }

  temps=ode->next_spring;
  while (temps!=&ode->next_spring) {
    MemCpy(ptr,&temps->start_saved_area,S_SIZE);
    ptr+=S_SIZE;
    temps=temps->next;
  }
  if (_size) *_size=cnt;
  return result;
}

Bool SpanWrite(CMathODE *ode)
{
  U8 *name,*buf;
  I64 size;
  Bool result=FALSE,old_silent=Silent(ON);
  MkDir("HOME/Span");
  Silent(old_silent);
  if (name=PopUpFileName("HOME/Span/Game.DAT")) {
    if (buf=SpanSave(ode,&size)) {
      FileWrite(name,buf,size);
      Free(buf);
      result=TRUE;
    }
    Free(name);
  }
  return result;
}

U0 SpanLoad(CMathODE *ode,U8 *src)
{
  I64 i;
  MyMass   *tempm;
  MySpring *temps;
  SpanHeader h;

  if (!src) return;

  MemCpy(&h,src,sizeof(SpanHeader));
  src+=sizeof(SpanHeader);

  for (i=0;i<h.num_masses;i++) {
    tempm=CAlloc(sizeof(MyMass));
    MemCpy(&tempm->start_saved_area,src,M_SIZE);
    src+=M_SIZE;
    QueIns(tempm,ode->last_mass);
  }

  for (i=0;i<h.num_springs;i++) {
    temps=CAlloc(sizeof(MySpring));
    MemCpy(&temps->start_saved_area,src,S_SIZE);
    src+=S_SIZE;
    QueIns(temps,ode->last_spring);
    temps->end1=MassFindNum(ode,temps->end1_num);
    temps->end2=MassFindNum(ode,temps->end2_num);
  }
}

U8 *SpanRead()
{
  U8 *src=NULL,*name;
  Bool old_silent=Silent(ON);
  MkDir("HOME/Span");
  Silent(old_silent);
  if (name=PopUpPickFile("HOME/Span")) {
    src=FileRead(name);
    Free(name);
  }
  return src;
}