fit-wlinear.eps

开放gsl矩阵运算

%!PS-Adobe-3.0 EPSF-3.0
%%Creator: GNU libplot drawing library 1.6
%%Title: PostScript plot
%%CreationDate: Sun Aug  6 11:26:34 2000
%%DocumentData: Clean7Bit
%%LanguageLevel: 1
%%Pages: 1
%%PageOrder: Ascend
%%Orientation: Portrait
%%BoundingBox: 85 180 500 576
%%DocumentNeededResources: 
%%DocumentSuppliedResources: procset GNU_libplot 1.0 0
%%EndComments

%%BeginDefaults
%%PageResources: 
%%EndDefaults

%%BeginProlog
%%EndProlog

%%BeginSetup
/DrawDict 50 dict def
DrawDict begin
%%BeginResource procset GNU_libplot 1.0 0
/none null def
/numGraphicParameters 17 def
/stringLimit 65535 def
/arrowHeight 8 def
/arrowWidth 4 def

/Begin { save numGraphicParameters dict begin } def
/End { end restore } def

/SetB {
dup type /nulltype eq {
pop
false /brushRightArrow idef
false /brushLeftArrow idef
true /brushNone idef
} {
/brushDashOffset idef
/brushDashArray idef
0 ne /brushRightArrow idef
0 ne /brushLeftArrow idef
/brushWidth idef
false /brushNone idef
} ifelse
} def

/SetCFg {
/fgblue idef
/fggreen idef
/fgred idef
} def

/SetCBg {
/bgblue idef
/bggreen idef
/bgred idef
} def

/SetF {
/printSize idef
/printFont idef
} def

/SetP {
dup type /nulltype eq {
pop true /patternNone idef
} {
/patternGrayLevel idef
patternGrayLevel -1 eq {
/patternString idef
} if
false /patternNone idef
} ifelse
} def

/BSpl {
0 begin
storexyn
newpath
n 1 gt {
0 0 0 0 0 0 1 1 true subspline
n 2 gt {
0 0 0 0 1 1 2 2 false subspline
1 1 n 3 sub {
/i exch def
i 1 sub dup i dup i 1 add dup i 2 add dup false subspline
} for
n 3 sub dup n 2 sub dup n 1 sub dup 2 copy false subspline
} if
n 2 sub dup n 1 sub dup 2 copy 2 copy false subspline
patternNone not brushLeftArrow not brushRightArrow not and and { ifill } if
brushNone not { istroke } if
0 0 1 1 leftarrow
n 2 sub dup 1 sub dup rightarrow
} if
end
} dup 0 4 dict put def

/Circ {
newpath
0 360 arc
closepath
patternNone not { ifill } if
brushNone not { istroke } if
} def

/CBSpl {
0 begin
dup 2 gt {
storexyn
newpath
n 1 sub dup 0 0 1 1 2 2 true subspline
1 1 n 3 sub {
/i exch def
i 1 sub dup i dup i 1 add dup i 2 add dup false subspline
} for
n 3 sub dup n 2 sub dup n 1 sub dup 0 0 false subspline
n 2 sub dup n 1 sub dup 0 0 1 1 false subspline
patternNone not { ifill } if
brushNone not { istroke } if
} {
Poly
} ifelse
end
} dup 0 4 dict put def

/Elli {
0 begin
newpath
4 2 roll
translate
scale
0 0 1 0 360 arc
closepath
patternNone not { ifill } if
brushNone not { istroke } if
end
} dup 0 1 dict put def

/Line {
0 begin
2 storexyn
newpath
x 0 get y 0 get moveto
x 1 get y 1 get lineto
brushNone not { istroke } if
0 0 1 1 leftarrow
0 0 1 1 rightarrow
end
} dup 0 4 dict put def

/MLine {
0 begin
storexyn
newpath
n 1 gt {
x 0 get y 0 get moveto
1 1 n 1 sub {
/i exch def
x i get y i get lineto
} for
patternNone not brushLeftArrow not brushRightArrow not and and { ifill } if
brushNone not { istroke } if
0 0 1 1 leftarrow
n 2 sub dup n 1 sub dup rightarrow
} if
end
} dup 0 4 dict put def

/Poly {
3 1 roll
newpath
moveto
-1 add
{ lineto } repeat
closepath
patternNone not { ifill } if
brushNone not { istroke } if
} def

/Rect {
0 begin
/t exch def
/r exch def
/b exch def
/l exch def
newpath
l b moveto
l t lineto
r t lineto
r b lineto
closepath
patternNone not { ifill } if
brushNone not { istroke } if
end
} dup 0 4 dict put def

/Text {
ishow
} def

/idef {
dup where { pop pop pop } { exch def } ifelse
} def

/ifill {
0 begin
gsave
patternGrayLevel -1 ne {
fgred bgred fgred sub patternGrayLevel mul add
fggreen bggreen fggreen sub patternGrayLevel mul add
fgblue bgblue fgblue sub patternGrayLevel mul add setrgbcolor
eofill
} {
eoclip
originalCTM setmatrix
pathbbox /t exch def /r exch def /b exch def /l exch def
/w r l sub ceiling cvi def
/h t b sub ceiling cvi def
/imageByteWidth w 8 div ceiling cvi def
/imageHeight h def
bgred bggreen bgblue setrgbcolor
eofill
fgred fggreen fgblue setrgbcolor
w 0 gt h 0 gt and {
l b translate w h scale
w h true [w 0 0 h neg 0 h] { patternproc } imagemask
} if
} ifelse
grestore
end
} dup 0 8 dict put def

/istroke {
gsave
brushDashOffset -1 eq {
[] 0 setdash
1 setgray
} {
brushDashArray brushDashOffset setdash
fgred fggreen fgblue setrgbcolor
} ifelse
brushWidth setlinewidth
originalCTM setmatrix
stroke
grestore
} def

/ishow {
0 begin
gsave
fgred fggreen fgblue setrgbcolor
/fontDict printFont findfont printSize scalefont dup setfont def
/descender fontDict begin 0 [FontBBox] 1 get FontMatrix end
transform exch pop def
/vertoffset 1 printSize sub descender sub def {
0 vertoffset moveto show
/vertoffset vertoffset printSize sub def
} forall
grestore
end
} dup 0 3 dict put def

/patternproc {
0 begin
/patternByteLength patternString length def
/patternHeight patternByteLength 8 mul sqrt cvi def
/patternWidth patternHeight def
/patternByteWidth patternWidth 8 idiv def
/imageByteMaxLength imageByteWidth imageHeight mul
stringLimit patternByteWidth sub min def
/imageMaxHeight imageByteMaxLength imageByteWidth idiv patternHeight idiv
patternHeight mul patternHeight max def
/imageHeight imageHeight imageMaxHeight sub store
/imageString imageByteWidth imageMaxHeight mul patternByteWidth add string def
0 1 imageMaxHeight 1 sub {
/y exch def
/patternRow y patternByteWidth mul patternByteLength mod def
/patternRowString patternString patternRow patternByteWidth getinterval def
/imageRow y imageByteWidth mul def
0 patternByteWidth imageByteWidth 1 sub {
/x exch def
imageString imageRow x add patternRowString putinterval
} for
} for
imageString
end
} dup 0 12 dict put def

/min {
dup 3 2 roll dup 4 3 roll lt { exch } if pop
} def

/max {
dup 3 2 roll dup 4 3 roll gt { exch } if pop
} def

/midpoint {
0 begin
/y1 exch def
/x1 exch def
/y0 exch def
/x0 exch def
x0 x1 add 2 div
y0 y1 add 2 div
end
} dup 0 4 dict put def

/thirdpoint {
0 begin
/y1 exch def
/x1 exch def
/y0 exch def
/x0 exch def
x0 2 mul x1 add 3 div
y0 2 mul y1 add 3 div
end
} dup 0 4 dict put def

/subspline {
0 begin
/movetoNeeded exch def
y exch get /y3 exch def
x exch get /x3 exch def
y exch get /y2 exch def
x exch get /x2 exch def
y exch get /y1 exch def
x exch get /x1 exch def
y exch get /y0 exch def
x exch get /x0 exch def
x1 y1 x2 y2 thirdpoint
/p1y exch def
/p1x exch def
x2 y2 x1 y1 thirdpoint
/p2y exch def
/p2x exch def
x1 y1 x0 y0 thirdpoint
p1x p1y midpoint
/p0y exch def
/p0x exch def
x2 y2 x3 y3 thirdpoint
p2x p2y midpoint
/p3y exch def
/p3x exch def
movetoNeeded { p0x p0y moveto } if
p1x p1y p2x p2y p3x p3y curveto
end
} dup 0 17 dict put def

/storexyn {
/n exch def
/y n array def
/x n array def
n 1 sub -1 0 {
/i exch def
y i 3 2 roll put
x i 3 2 roll put
} for
} def

/arrowhead {
0 begin
transform originalCTM itransform
/taily exch def
/tailx exch def
transform originalCTM itransform
/tipy exch def
/tipx exch def
/dy tipy taily sub def
/dx tipx tailx sub def
/angle dx 0 ne dy 0 ne or { dy dx atan } { 90 } ifelse def
gsave
originalCTM setmatrix
tipx tipy translate
angle rotate
newpath
arrowHeight neg arrowWidth 2 div moveto
0 0 lineto
arrowHeight neg arrowWidth 2 div neg lineto
patternNone not {
originalCTM setmatrix
/padtip arrowHeight 2 exp 0.25 arrowWidth 2 exp mul add sqrt brushWidth mul
arrowWidth div def
/padtail brushWidth 2 div def
tipx tipy translate
angle rotate
padtip 0 translate
arrowHeight padtip add padtail add arrowHeight div dup scale
arrowheadpath
ifill
} if
brushNone not {
originalCTM setmatrix
tipx tipy translate
angle rotate
arrowheadpath
istroke
} if
grestore
end
} dup 0 9 dict put def

/arrowheadpath {
newpath
arrowHeight neg arrowWidth 2 div moveto
0 0 lineto
arrowHeight neg arrowWidth 2 div neg lineto
} def

/leftarrow {
0 begin
y exch get /taily exch def
x exch get /tailx exch def
y exch get /tipy exch def
x exch get /tipx exch def
brushLeftArrow { tipx tipy tailx taily arrowhead } if
end
} dup 0 4 dict put def

/rightarrow {
0 begin
y exch get /tipy exch def
x exch get /tipx exch def
y exch get /taily exch def
x exch get /tailx exch def
brushRightArrow { tipx tipy tailx taily arrowhead } if
end
} dup 0 4 dict put def
%%EndResource
%%EndSetup

%%Page: 1 1
%%PageResources: 
%%PageBoundingBox: 85 180 500 576
%%BeginPageSetup
%I Idraw 8

Begin
%I b u
%I cfg u
%I cbg u
%I f u
%I p u
%I t
[ 1 0 0 1 0 0 ] concat
/originalCTM matrix currentmatrix def
/trueoriginalCTM matrix currentmatrix def
%%EndPageSetup

Begin %I Poly
[0.6776471 0 0 0.6776471 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
0 setlinecap 0 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 4
2304 2304
2304 9216
9216 9216
9216 2304
4 Poly
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 4
1930 2107
1952 2118
1985 2151
1985 1920
4 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 2
1974 2140
1974 1920
2 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 2
1930 1920
2029 1920
2 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 28
2260 2074
2249 2041
2227 2019
2194 2008
2183 2008
2150 2019
2128 2041
2117 2074
2117 2085
2128 2118
2150 2140
2183 2151
2205 2151
2238 2140
2260 2118
2271 2085
2271 2019
2260 1975
2249 1953
2227 1931
2194 1920
2161 1920
2139 1931
2128 1953
2128 1964
2139 1975
2150 1964
2139 1953
28 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 8
2183 2008
2161 2019
2139 2041
2128 2074
2128 2085
2139 2118
2161 2140
2183 2151
8 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 9
2205 2151
2227 2140
2249 2118
2260 2085
2260 2019
2249 1975
2238 1953
2216 1931
2194 1920
9 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 28
2469 2118
2458 2107
2469 2096
2480 2107
2480 2118
2469 2140
2447 2151
2414 2151
2381 2140
2359 2118
2348 2096
2337 2052
2337 1986
2348 1953
2370 1931
2403 1920
2425 1920
2458 1931
2480 1953
2491 1986
2491 1997
2480 2030
2458 2052
2425 2063
2414 2063
2381 2052
2359 2030
2348 1997
28 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 9
2414 2151
2392 2140
2370 2118
2359 2096
2348 2052
2348 1986
2359 1953
2381 1931
2403 1920
9 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 8
2425 1920
2447 1931
2469 1953
2480 1986
2480 1997
2469 2030
2447 2052
2425 2063
8 MLine
End

Begin %I Poly
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 16
2623 2151
2590 2140
2568 2107
2557 2052
2557 2019
2568 1964
2590 1931
2623 1920
2645 1920
2678 1931
2700 1964
2711 2019
2711 2052
2700 2107
2678 2140
2645 2151
16 Poly
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 10
2623 2151
2601 2140
2590 2129
2579 2107
2568 2052
2568 2019
2579 1964
2590 1942
2601 1931
2623 1920
10 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 10
2645 1920
2667 1931
2678 1942
2689 1964
2700 2019
2700 2052
2689 2107
2678 2129
2667 2140
2645 2151
10 MLine
End

Begin %I MLine
[0.6776471 0 0 0.6776471 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
0 setlinecap 0 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 2
2304 9216
2304 9078
2 MLine
End

Begin %I MLine
[0.6776471 0 0 0.6776471 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
0 setlinecap 0 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 2
2304 2304
2304 2442
2 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 4
3313 2107
3335 2118
3368 2151
3368 1920
4 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 2
3357 2140
3357 1920
2 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB
%I cfg Black
0 0 0 SetCFg
%I cbg White
1 1 1 SetCBg
%I p
none SetP
%I t
[0.05 0 0 0.05 18 108 ] concat
%I 2
3313 1920
3411 1920
2 MLine
End

Begin %I MLine
[0.7807418 0 0 0.7807418 0 0 ] trueoriginalCTM originalCTM
concatmatrix pop
1 setlinecap 1 setlinejoin
%I b 65535
1 0 0 [ ] 0 SetB