pro dlafitsys, d,ss,r

s = indgen(120)/2 - 30
a12 = ( s[ss[0] + 60.0])
a14 = ss[1]
a15 = ss[2]
a16 = ss[3]
b12 = (-s[ss[0] + 61.0])
b24 = -ss[4]
b25 = -ss[5]
b26 = -ss[6]
a46 = -ss[8]
b56 = -ss[9]
a45 = ( s[ss[7]+60.0])
b45 = (-s[ss[7]+61.0])

c = fltarr(10)
c[0] = d[0] - d[2] + a14
c[1] = d[0] - d[3] + a15
c[2] = d[0] - d[4] + a16
c[3] = d[1] - d[2] + b24
c[4] = d[1] - d[3] + b25
c[5] = d[1] - d[4] + b26
c[6] = d[1] - d[0] + b12 - a12
c[7] = d[2] - d[4] + a46
c[8] = d[3] - d[4] + b56
c[9] = d[3] - d[2] + b45 - a45
M = [[ 1, 0,-1, 0, 0], $
     [ 1, 0, 0,-1, 0], $
     [ 1, 0, 0, 0,-1], $
     [ 0, 1,-1, 0, 0], $
     [ 0, 1, 0,-1, 0], $
     [ 0, 1, 0, 0,-1], $
     [-1, 1, 0, 0, 0], $
     [ 0, 0, 1, 0,-1], $
     [ 0, 0, 0, 1,-1], $
     [ 0, 0,-1, 1, 0]]


svdc,m, w, u, v,/double
n = n_elements(w)
wp = fltarr(n,n)
for k = 0, n-1 do $
   if abs(w(k)) ge 1.0e-5 then wp(k,k) = 1.0/w(k)
t = v ## wp ## TRANSPOSE(U) ## C

r = t[0,*] - t[0,4]

end