gmt-fork/test/gmtmath/lsfit.ps

%!PS-Adobe-3.0
%%BoundingBox: 0 0 612 792
%%HiResBoundingBox: 0 0 612.0000 792.0000             
%%Title: GMT v6.0.0_831abca_2019.06.27 [64-bit] Document from psxy
%%Creator: GMT6
%%For: unknown
%%DocumentNeededResources: font Helvetica
%%CreationDate: Thu Jun 27 12:16:20 2019
%%LanguageLevel: 2
%%DocumentData: Clean7Bit
%%Orientation: Portrait
%%Pages: 1
%%EndComments

%%BeginProlog
250 dict begin
/! {bind def} bind def
/# {load def}!
/A /setgray #
/B /setdash #
/C /setrgbcolor #
/D /rlineto #
/E {dup stringwidth pop}!
/F /fill #
/G /rmoveto #
/H /sethsbcolor #
/I /setpattern #
/K /setcmykcolor #
/L /lineto #
/M /moveto #
/N /newpath #
/P /closepath #
/R /rotate #
/S /stroke #
/T /translate #
/U /grestore #
/V /gsave #
/W /setlinewidth #
/Y {findfont exch scalefont setfont}!
/Z /show #
/FP {true charpath flattenpath}!
/MU {matrix setmatrix}!
/MS {/SMat matrix currentmatrix def}!
/MR {SMat setmatrix}!
/edef {exch def}!
/FS {/fc edef /fs {V fc F U} def}!
/FQ {/fs {} def}!
/O0 {/os {N} def}!
/O1 {/os {P S} def}!
/FO {fs os}!
/Sa {M MS dup 0 exch G 0.726542528 mul -72 R dup 0 D 4 {72 R dup 0 D -144 R dup 0 D} repeat pop MR FO}!
/Sb {M dup 0 D exch 0 exch D neg 0 D FO}!
/SB {MS T /BoxR edef /BoxW edef /BoxH edef BoxR 0 M
  BoxW 0 BoxW BoxH BoxR arct BoxW BoxH 0 BoxH BoxR arct 0 BoxH 0 0 BoxR arct 0 0 BoxW 0 BoxR arct MR FO}!
/Sc {N 3 -1 roll 0 360 arc FO}!
/Sd {M 4 {dup} repeat 0 G neg dup dup D exch D D FO}!
/Se {N MS T R scale 0 0 1 0 360 arc MR FO}!
/Sg {M MS 22.5 R dup 0 exch G -22.5 R 0.765366865 mul dup 0 D 6 {-45 R dup 0 D} repeat pop MR FO}!
/Sh {M MS dup 0 G -120 R dup 0 D 4 {-60 R dup 0 D} repeat pop MR FO}!
/Si {M MS dup neg 0 exch G 60 R 1.732050808 mul dup 0 D 120 R 0 D MR FO}!
/Sj {M MS R dup -2 div 2 index -2 div G dup 0 D exch 0 exch D neg 0 D MR FO}!
/Sn {M MS dup 0 exch G -36 R 1.175570505 mul dup 0 D 3 {-72 R dup 0 D} repeat pop MR FO}!
/Sp {N 3 -1 roll 0 360 arc fs N}!
/SP {M {D} repeat FO}!
/Sr {M dup -2 div 2 index -2 div G dup 0 D exch 0 exch D neg 0 D FO}!
/SR {MS T /BoxR edef /BoxW edef /BoxH edef BoxR BoxW -2 div BoxH -2 div T BoxR 0 M
  BoxW 0 BoxW BoxH BoxR arct BoxW BoxH 0 BoxH BoxR arct 0 BoxH 0 0 BoxR arct 0 0 BoxW 0 BoxR arct MR FO}!
/Ss {M 1.414213562 mul dup dup dup -2 div dup G 0 D 0 exch D neg 0 D FO}!
/St {M MS dup 0 exch G -60 R 1.732050808 mul dup 0 D -120 R 0 D MR FO}!
/SV {0 exch M 0 D D D D D 0 D FO}!
/Sv {0 0 M D D 0 D D D D D 0 D D FO}!
/Sw {2 copy M 5 2 roll arc FO}!
/Sx {M 1.414213562 mul 5 {dup} repeat -2 div dup G D neg 0 G neg D S}!
/Sy {M dup 0 exch G dup -2 mul dup 0 exch D S}!
/S+ {M dup 0 G dup -2 mul dup 0 D exch dup G 0 exch D S}!
/S- {M dup 0 G dup -2 mul dup 0 D S}!
/sw {stringwidth pop}!
/sh {V MU 0 0 M FP pathbbox N 4 1 roll pop pop pop U}!
/sd {V MU 0 0 M FP pathbbox N pop pop exch pop U}!
/sH {V MU 0 0 M FP pathbbox N exch pop exch sub exch pop U}!
/sb {E exch sh}!
/bl {}!
/bc {E -2 div 0 G}!
/br {E neg 0 G}!
/ml {dup 0 exch sh -2 div G}!
/mc {dup E -2 div exch sh -2 div G}!
/mr {dup E neg exch sh -2 div G}!
/tl {dup 0 exch sh neg G}!
/tc {dup E -2 div exch sh neg G}!
/tr {dup E neg exch sh neg G}!
/mx {2 copy lt {exch} if pop}!
/PSL_xorig 0 def /PSL_yorig 0 def
/TM {2 copy T PSL_yorig add /PSL_yorig edef PSL_xorig add /PSL_xorig edef}!
/PSL_reencode {findfont dup length dict begin
  {1 index /FID ne {def}{pop pop} ifelse} forall
  exch /Encoding edef currentdict end definefont pop
}!
/PSL_eps_begin {
  /PSL_eps_state save def
  /PSL_dict_count countdictstack def
  /PSL_op_count count 1 sub def
  userdict begin
  /showpage {} def
  0 setgray 0 setlinecap 1 setlinewidth
  0 setlinejoin 10 setmiterlimit [] 0 setdash newpath
  /languagelevel where
  {pop languagelevel 1 ne {false setstrokeadjust false setoverprint} if} if
}!
/PSL_eps_end {
  count PSL_op_count sub {pop} repeat
  countdictstack PSL_dict_count sub {end} repeat
  PSL_eps_state restore
}!
/PSL_transp {
  /.setopacityalpha where {pop .setblendmode .setopacityalpha}{
  /pdfmark where {pop [ /BM exch /CA exch dup /ca exch /SetTransparency pdfmark}
  {pop pop} ifelse} ifelse
}!
/Standard+_Encoding [
/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef
/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef
/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef	/.notdef
/.notdef	/threequarters	/threesuperior	/trademark	/twosuperior	/yacute		/ydieresis	/zcaron
/space		/exclam		/quotedbl	/numbersign	/dollar		/percent	/ampersand	/quoteright
/parenleft	/parenright	/asterisk	/plus		/comma		/hyphen		/period		/slash
/zero		/one		/two		/three		/four		/five		/six		/seven
/eight		/nine		/colon		/semicolon	/less		/equal		/greater	/question
/at		/A		/B		/C		/D		/E		/F		/G
/H		/I		/J		/K		/L		/M		/N		/O
/P		/Q		/R		/S		/T		/U		/V		/W
/X		/Y		/Z		/bracketleft	/backslash	/bracketright	/asciicircum	/underscore
/quoteleft	/a		/b		/c		/d		/e		/f		/g
/h		/i		/j		/k		/l		/m		/n		/o
/p		/q		/r		/s		/t		/u		/v		/w
/x		/y		/z		/braceleft	/bar		/braceright	/asciitilde	/florin
/Atilde		/Ccedilla	/Eth		/Lslash		/Ntilde		/Otilde		/Scaron		/Thorn
/Yacute		/Ydieresis	/Zcaron		/atilde		/brokenbar	/ccedilla	/copyright	/degree
/divide		/eth		/logicalnot	/lslash		/minus		/mu		/multiply	/ntilde
/onehalf	/onequarter	/onesuperior	/otilde		/plusminus	/registered	/scaron		/thorn
/.notdef	/exclamdown	/cent		/sterling	/fraction	/yen		/florin		/section
/currency	/quotesingle	/quotedblleft	/guillemotleft	/guilsinglleft	/guilsinglright	/fi		/fl
/Aacute		/endash		/dagger		/daggerdbl	/periodcentered	/Acircumflex	/paragraph	/bullet
/quotesinglbase	/quotedblbase	/quotedblright	/guillemotright	/ellipsis	/perthousand	/Adieresis	/questiondown
/Agrave		/grave		/acute		/circumflex	/tilde		/macron		/breve		/dotaccent
/dieresis	/Eacute		/ring		/cedilla	/Ecircumflex	/hungarumlaut	/ogonek		/caron
/emdash		/Edieresis	/Egrave		/Iacute		/Icircumflex	/Idieresis	/Igrave		/Oacute
/Ocircumflex	/Odieresis	/Ograve		/Uacute		/Ucircumflex	/Udieresis	/Ugrave		/aacute
/acircumflex	/AE		/adieresis	/ordfeminine	/agrave		/eacute		/ecircumflex	/edieresis
/egrave		/Oslash		/OE		/ordmasculine	/iacute		/icircumflex	/idieresis	/igrave
/oacute		/ae		/ocircumflex	/odieresis	/ograve		/dotlessi	/uacute		/ucircumflex
/udieresis	/oslash		/oe		/germandbls	/ugrave		/Aring		/aring		/ydieresis
] def
/PSL_font_encode 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 39 array astore def
/F0 {/Helvetica Y}!
/F1 {/Helvetica-Bold Y}!
/F2 {/Helvetica-Oblique Y}!
/F3 {/Helvetica-BoldOblique Y}!
/F4 {/Times-Roman Y}!
/F5 {/Times-Bold Y}!
/F6 {/Times-Italic Y}!
/F7 {/Times-BoldItalic Y}!
/F8 {/Courier Y}!
/F9 {/Courier-Bold Y}!
/F10 {/Courier-Oblique Y}!
/F11 {/Courier-BoldOblique Y}!
/F12 {/Symbol Y}!
/F13 {/AvantGarde-Book Y}!
/F14 {/AvantGarde-BookOblique Y}!
/F15 {/AvantGarde-Demi Y}!
/F16 {/AvantGarde-DemiOblique Y}!
/F17 {/Bookman-Demi Y}!
/F18 {/Bookman-DemiItalic Y}!
/F19 {/Bookman-Light Y}!
/F20 {/Bookman-LightItalic Y}!
/F21 {/Helvetica-Narrow Y}!
/F22 {/Helvetica-Narrow-Bold Y}!
/F23 {/Helvetica-Narrow-Oblique Y}!
/F24 {/Helvetica-Narrow-BoldOblique Y}!
/F25 {/NewCenturySchlbk-Roman Y}!
/F26 {/NewCenturySchlbk-Italic Y}!
/F27 {/NewCenturySchlbk-Bold Y}!
/F28 {/NewCenturySchlbk-BoldItalic Y}!
/F29 {/Palatino-Roman Y}!
/F30 {/Palatino-Italic Y}!
/F31 {/Palatino-Bold Y}!
/F32 {/Palatino-BoldItalic Y}!
/F33 {/ZapfChancery-MediumItalic Y}!
/F34 {/ZapfDingbats Y}!
/F35 {/Ryumin-Light-EUC-H Y}!
/F36 {/Ryumin-Light-EUC-V Y}!
/F37 {/GothicBBB-Medium-EUC-H Y}!
/F38 {/GothicBBB-Medium-EUC-V Y}!
/PSL_pathtextdict 26 dict def
/PSL_pathtext
  {PSL_pathtextdict begin
    /ydepth exch def
    /textheight exch def
    /just exch def
    /offset exch def
    /str exch def
    /pathdist 0 def
    /setdist offset def
    /charcount 0 def
    /justy just 4 idiv textheight mul 2 div neg ydepth sub def
    V flattenpath
	{movetoproc} {linetoproc}
	{curvetoproc} {closepathproc}
	pathforall
    U N
    end
  } def
PSL_pathtextdict begin
/movetoproc
  { /newy exch def /newx exch def
    /firstx newx def /firsty newy def
    /ovr 0 def
    newx newy transform
    /cpy exch def /cpx exch def
  } def
/linetoproc
  { /oldx newx def /oldy newy def
    /newy exch def /newx exch def
    /dx newx oldx sub def
    /dy newy oldy sub def
    /dist dx dup mul dy dup mul add sqrt def
    dist 0 ne
    { /dsx dx dist div ovr mul def
      /dsy dy dist div ovr mul def
      oldx dsx add oldy dsy add transform
      /cpy exch def /cpx exch def
      /pathdist pathdist dist add def
      {setdist pathdist le
	  {charcount str length lt
	      {setchar} {exit} ifelse}
	  { /ovr setdist pathdist sub def
	    exit}
	  ifelse
      } loop
    } if
  } def
/curvetoproc
  { (ERROR: No curveto's after flattenpath!)
    print
  } def
/closepathproc
  {firstx firsty linetoproc
    firstx firsty movetoproc
  } def
/setchar
  { /char str charcount 1 getinterval def
    /charcount charcount 1 add def
    /charwidth char stringwidth pop def
    V cpx cpy itransform T
      dy dx atan R
      0 justy M
      char show
      0 justy neg G
      currentpoint transform
      /cpy exch def /cpx exch def
    U /setdist setdist charwidth add def
  } def
end
/PSL_set_label_heights
{
  /PSL_n_labels_minus_1 PSL_n_labels 1 sub def
  /PSL_heights PSL_n_labels array def
  0 1 PSL_n_labels_minus_1
  { /psl_k exch def
    /psl_label PSL_label_str psl_k get def
    PSL_label_font psl_k get cvx exec
    psl_label sH /PSL_height edef
    PSL_heights psl_k PSL_height put
  } for
} def
/PSL_curved_path_labels
{ /psl_bits exch def
  /PSL_placetext psl_bits 2 and 2 eq def
  /PSL_clippath psl_bits 4 and 4 eq def
  /PSL_strokeline false def
  /PSL_fillbox psl_bits 128 and 128 eq def
  /PSL_drawbox psl_bits 256 and 256 eq def
  /PSL_n_paths1 PSL_n_paths 1 sub def
  /PSL_usebox PSL_fillbox PSL_drawbox or def
  PSL_clippath {clipsave N clippath} if
  /psl_k 0 def
  /psl_p 0 def
  0 1 PSL_n_paths1
  { /psl_kk exch def
    /PSL_n PSL_path_n  psl_kk get def
    /PSL_m PSL_label_n psl_kk get def
    /PSL_x PSL_path_x psl_k PSL_n getinterval def
    /PSL_y PSL_path_y psl_k PSL_n getinterval def
    /PSL_node_tmp PSL_label_node psl_p PSL_m getinterval def
    /PSL_angle_tmp PSL_label_angle psl_p PSL_m getinterval def
    /PSL_str_tmp PSL_label_str psl_p PSL_m getinterval def
    /PSL_fnt_tmp PSL_label_font psl_p PSL_m getinterval def
    PSL_curved_path_label
    /psl_k psl_k PSL_n add def
    /psl_p psl_p PSL_m add def
  } for
  PSL_clippath {PSL_eoclip} if N
} def
/PSL_curved_path_label
{
  /PSL_n1 PSL_n 1 sub def
  /PSL_m1 PSL_m 1 sub def
  PSL_CT_calcstringwidth
  PSL_CT_calclinedist
  PSL_CT_excludelabels
  PSL_CT_addcutpoints
  /PSL_nn1 PSL_nn 1 sub def
  /n 0 def
  /k 0 def
  /j 0 def
  /PSL_seg 0 def
  /PSL_xp PSL_nn array def
  /PSL_yp PSL_nn array def
  PSL_xp 0 PSL_xx 0 get put
  PSL_yp 0 PSL_yy 0 get put
  1 1 PSL_nn1
  { /i exch def
    /node_type PSL_kind i get def
    /j j 1 add def
    PSL_xp j PSL_xx i get put
    PSL_yp j PSL_yy i get put
    node_type 1 eq
    {n 0 eq
      {PSL_CT_drawline}
      {	PSL_CT_reversepath
	PSL_CT_textline} ifelse
      /j 0 def
      PSL_xp j PSL_xx i get put
      PSL_yp j PSL_yy i get put
    } if
  } for
  n 0 eq {PSL_CT_drawline} if
} def
/PSL_CT_textline
{ PSL_fnt k get cvx exec
  /PSL_height PSL_heights k get def
  PSL_placetext	{PSL_CT_placelabel} if
  PSL_clippath {PSL_CT_clippath} if
  /n 0 def /k k 1 add def
} def
/PSL_CT_calcstringwidth
{ /PSL_width_tmp PSL_m array def
  0 1 PSL_m1
  { /i exch def
    PSL_fnt_tmp i get cvx exec
    PSL_width_tmp i PSL_str_tmp i get stringwidth pop put
  } for
} def
/PSL_CT_calclinedist
{ /PSL_newx PSL_x 0 get def
  /PSL_newy PSL_y 0 get def
  /dist 0.0 def
  /PSL_dist PSL_n array def
  PSL_dist 0 0.0 put
  1 1 PSL_n1
  { /i exch def
    /PSL_oldx PSL_newx def
    /PSL_oldy PSL_newy def
    /PSL_newx PSL_x i get def
    /PSL_newy PSL_y i get def
    /dx PSL_newx PSL_oldx sub def
    /dy PSL_newy PSL_oldy sub def
    /dist dist dx dx mul dy dy mul add sqrt add def
    PSL_dist i dist put
  } for
} def
/PSL_CT_excludelabels
{ /k 0 def
  /PSL_width PSL_m array def
  /PSL_angle PSL_m array def
  /PSL_node PSL_m array def
  /PSL_str PSL_m array def
  /PSL_fnt PSL_m array def
  /lastdist PSL_dist PSL_n1 get def
  0 1 PSL_m1
  { /i exch def
    /dist PSL_dist PSL_node_tmp i get get def
    /halfwidth PSL_width_tmp i get 2 div PSL_gap_x add def
    /L_dist dist halfwidth sub def
    /R_dist dist halfwidth add def
    L_dist 0 gt R_dist lastdist lt and
    {
      PSL_width k PSL_width_tmp i get put
      PSL_node k PSL_node_tmp i get put
      PSL_angle k PSL_angle_tmp i get put
      PSL_str k PSL_str_tmp i get put
      PSL_fnt k PSL_fnt_tmp i get put
      /k k 1 add def
    } if
  } for
  /PSL_m k def
  /PSL_m1 PSL_m 1 sub def
} def
/PSL_CT_addcutpoints
{ /k 0 def
  /PSL_nc PSL_m 2 mul 1 add def
  /PSL_cuts PSL_nc array def
  /PSL_nc1 PSL_nc 1 sub def
  0 1 PSL_m1
  { /i exch def
    /dist PSL_dist PSL_node i get get def
    /halfwidth PSL_width i get 2 div PSL_gap_x add def
    PSL_cuts k dist halfwidth sub put
    /k k 1 add def
    PSL_cuts k dist halfwidth add put
    /k k 1 add def
  } for
  PSL_cuts k 100000.0 put
  /PSL_nn PSL_n PSL_m 2 mul add def
  /PSL_xx PSL_nn array def
  /PSL_yy PSL_nn array def
  /PSL_kind PSL_nn array def
  /j 0 def
  /k 0 def
  /dist 0.0 def
  0 1 PSL_n1
  { /i exch def
    /last_dist dist def
    /dist PSL_dist i get def
    k 1 PSL_nc1
    { /kk exch def
      /this_cut PSL_cuts kk get def
      dist this_cut gt
      { /ds dist last_dist sub def
	/f ds 0.0 eq {0.0} {dist this_cut sub ds div} ifelse def
	/i1 i 0 eq {0} {i 1 sub} ifelse def
	PSL_xx j PSL_x i get dup PSL_x i1 get sub f mul sub put
	PSL_yy j PSL_y i get dup PSL_y i1 get sub f mul sub put
	PSL_kind j 1 put
	/j j 1 add def
	/k k 1 add def
      } if
    } for
    dist PSL_cuts k get le
    {PSL_xx j PSL_x i get put PSL_yy j PSL_y i get put
      PSL_kind j 0 put
      /j j 1 add def
    } if
  } for
} def
/PSL_CT_reversepath
{PSL_xp j get PSL_xp 0 get lt
  {0 1 j 2 idiv
    { /left exch def
      /right j left sub def
      /tmp PSL_xp left get def
      PSL_xp left PSL_xp right get put
      PSL_xp right tmp put
      /tmp PSL_yp left get def
      PSL_yp left PSL_yp right get put
      PSL_yp right tmp put
    } for
  } if
} def
/PSL_CT_placelabel
{
  /PSL_just PSL_label_justify k get def
  /PSL_height PSL_heights k get def
  /psl_label PSL_str k get def
  /psl_depth psl_label sd def
  PSL_usebox
  {PSL_CT_clippath
    PSL_fillbox
    {V PSL_setboxrgb fill U} if
    PSL_drawbox
    {V PSL_setboxpen S U} if N
  } if
  PSL_CT_placeline psl_label PSL_gap_x PSL_just PSL_height psl_depth PSL_pathtext
} def
/PSL_CT_clippath
{
  /H PSL_height 2 div PSL_gap_y add def
  /xoff j 1 add array def
  /yoff j 1 add array def
  /angle 0 def
  0 1 j {
    /ii exch def
    /x PSL_xp ii get def
    /y PSL_yp ii get def
    ii 0 eq {
      /x1 PSL_xp 1 get def
      /y1 PSL_yp 1 get def
      /dx x1 x sub def
      /dy y1 y sub def
    }
    { /i1 ii 1 sub def
      /x1 PSL_xp i1 get def
      /y1 PSL_yp i1 get def
      /dx x x1 sub def
      /dy y y1 sub def
    } ifelse
    dx 0.0 eq dy 0.0 eq and not
    { /angle dy dx atan 90 add def} if
    /sina angle sin def
    /cosa angle cos def
    xoff ii H cosa mul put
    yoff ii H sina mul put
  } for
  PSL_xp 0 get xoff 0 get add PSL_yp 0 get yoff 0 get add M
  1 1 j {
    /ii exch def
    PSL_xp ii get xoff ii get add PSL_yp ii get yoff ii get add L
  } for
  j -1 0 {
    /ii exch def
    PSL_xp ii get xoff ii get sub PSL_yp ii get yoff ii get sub L
  } for P
} def
/PSL_CT_drawline
{
  /str 20 string def
  PSL_strokeline
  {PSL_CT_placeline S} if
  /PSL_seg PSL_seg 1 add def
  /n 1 def
} def
/PSL_CT_placeline
{PSL_xp 0 get PSL_yp 0 get M
  1 1 j { /ii exch def PSL_xp ii get PSL_yp ii get L} for
} def
/PSL_draw_path_lines
{
  /PSL_n_paths1 PSL_n_paths 1 sub def
  V
  /psl_start 0 def
  0 1 PSL_n_paths1
  { /psl_k exch def
    /PSL_n PSL_path_n psl_k get def
    /PSL_n1 PSL_n 1 sub def
    PSL_path_pen psl_k get cvx exec
    N
    PSL_path_x psl_start get PSL_path_y psl_start get M
    1 1 PSL_n1
    { /psl_i exch def
      /psl_kk psl_i psl_start add def
      PSL_path_x psl_kk get PSL_path_y psl_kk get L
    } for
    /psl_xclose PSL_path_x psl_kk get PSL_path_x psl_start get sub def
    /psl_yclose PSL_path_y psl_kk get PSL_path_y psl_start get sub def
    psl_xclose 0 eq psl_yclose 0 eq and { P } if
    S
    /psl_start psl_start PSL_n add def
  } for
  U
} def
/PSL_straight_path_labels
{
  /psl_bits exch def
  /PSL_placetext psl_bits 2 and 2 eq def
  /PSL_rounded psl_bits 32 and 32 eq def
  /PSL_fillbox psl_bits 128 and 128 eq def
  /PSL_drawbox psl_bits 256 and 256 eq def
  /PSL_n_labels_minus_1 PSL_n_labels 1 sub def
  /PSL_usebox PSL_fillbox PSL_drawbox or def
  0 1 PSL_n_labels_minus_1
  { /psl_k exch def
    PSL_ST_prepare_text
    PSL_usebox
    {  PSL_rounded
        {PSL_ST_textbox_round}
        {PSL_ST_textbox_rect}
      ifelse
      PSL_fillbox {V PSL_setboxrgb fill U} if
      PSL_drawbox {V PSL_setboxpen S U} if
      N
    } if
    PSL_placetext {PSL_ST_place_label} if
  } for
} def
/PSL_straight_path_clip
{
  /psl_bits exch def
  /PSL_rounded psl_bits 32 and 32 eq def
  /PSL_n_labels_minus_1 PSL_n_labels 1 sub def
  N clipsave clippath
  0 1 PSL_n_labels_minus_1
  { /psl_k exch def
    PSL_ST_prepare_text
    PSL_rounded
      {PSL_ST_textbox_round}
      {PSL_ST_textbox_rect}
    ifelse
  } for
  PSL_eoclip N
} def
/PSL_ST_prepare_text
{
  /psl_xp PSL_txt_x psl_k get def
  /psl_yp PSL_txt_y psl_k get def
  /psl_label PSL_label_str psl_k get def
  PSL_label_font psl_k get cvx exec
  /PSL_height PSL_heights psl_k get def
  /psl_boxH PSL_height PSL_gap_y 2 mul add def
  /PSL_just PSL_label_justify psl_k get def
  /PSL_justx PSL_just 4 mod 1 sub 2 div neg def
  /PSL_justy PSL_just 4 idiv 2 div neg def
  /psl_SW psl_label stringwidth pop def
  /psl_boxW psl_SW PSL_gap_x 2 mul add def
  /psl_x0 psl_SW PSL_justx mul def
  /psl_y0 PSL_justy PSL_height mul def
  /psl_angle PSL_label_angle psl_k get def
} def
/PSL_ST_textbox_rect
{
  psl_xp psl_yp T psl_angle R psl_x0 psl_y0 T
  PSL_gap_x neg PSL_gap_y neg M
  0 psl_boxH D psl_boxW 0 D 0 psl_boxH neg D P
  psl_x0 neg psl_y0 neg T psl_angle neg R psl_xp neg psl_yp neg T
} def
/PSL_ST_textbox_round
{
  /psl_BoxR PSL_gap_x PSL_gap_y lt {PSL_gap_x} {PSL_gap_y} ifelse def
  /psl_xd PSL_gap_x psl_BoxR sub def
  /psl_yd PSL_gap_y psl_BoxR sub def
  /psl_xL PSL_gap_x neg def
  /psl_yB PSL_gap_y neg def
  /psl_yT psl_boxH psl_yB add def
  /psl_H2 PSL_height psl_yd 2 mul add def
  /psl_W2 psl_SW psl_xd 2 mul add def
  /psl_xR psl_xL psl_boxW add def
  /psl_x0 psl_SW PSL_justx mul def
  psl_xp psl_yp T psl_angle R psl_x0 psl_y0 T
  psl_xL psl_yd M
  psl_xL psl_yT psl_xR psl_yT psl_BoxR arct psl_W2 0 D
  psl_xR psl_yT psl_xR psl_yB psl_BoxR arct 0 psl_H2 neg D
  psl_xR psl_yB psl_xL psl_yB psl_BoxR arct psl_W2 neg 0 D
  psl_xL psl_yB psl_xL psl_yd psl_BoxR arct P
  psl_x0 neg psl_y0 neg T psl_angle neg R psl_xp neg psl_yp neg T
} def
/PSL_ST_place_label
{
    V psl_xp psl_yp T psl_angle R
    psl_SW PSL_justx mul psl_y0 M
    psl_label dup sd neg 0 exch G show
    U
} def
/PSL_nclip 0 def
/PSL_clip {clip /PSL_nclip PSL_nclip 1 add def} def
/PSL_eoclip {eoclip /PSL_nclip PSL_nclip 1 add def} def
/PSL_cliprestore {cliprestore /PSL_nclip PSL_nclip 1 sub def} def
%%EndProlog

%%BeginSetup
/PSLevel /languagelevel where {pop languagelevel} {1} ifelse def
PSLevel 1 gt { << /PageSize [612 792] /ImagingBBox null >> setpagedevice } if
%%EndSetup

%%Page: 1 1

%%BeginPageSetup
V 0.06 0.06 scale
%%EndPageSetup

/PSL_page_xsize 10200 def
/PSL_page_ysize 13200 def
/PSL_completion {} def
/PSL_movie_completion {} def
%PSL_End_Header
gsave
0 A
FQ
O0
-3900 PSL_xorig sub PSL_page_xsize 2 div add 1200 TM

% PostScript produced by:
%@GMT: gmt psxy -R-3/3/0/4 -JX6.5i/4i -P -Baf -BWSne @stepdata.txt -Sc0.15c -Gblue -K -Xc
%@PROJ: xy -3.00000000 3.00000000 0.00000000 4.00000000 -3.000 3.000 0.000 4.000 +xy
%GMTBoundingBox: -234 72 468 288
%%BeginObject PSL_Layer_1
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
25 W
2 setlinecap
N 0 4800 M 0 -4800 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 0 M -83 0 D S
N 0 1200 M -83 0 D S
N 0 2400 M -83 0 D S
N 0 3600 M -83 0 D S
N 0 4800 M -83 0 D S
/PSL_AH0 0
/MM {neg exch M} def
PSL_font_encode 0 get 0 eq {Standard+_Encoding /Helvetica /Helvetica PSL_reencode PSL_font_encode 0 1 put} if
200 F0
(0) sw mx
(1) sw mx
(2) sw mx
(3) sw mx
(4) sw mx
def
/PSL_A0_y PSL_A0_y 83 add def
0 PSL_A0_y MM
(0) mr Z
1200 PSL_A0_y MM
(1) mr Z
2400 PSL_A0_y MM
(2) mr Z
3600 PSL_A0_y MM
(3) mr Z
4800 PSL_A0_y MM
(4) mr Z
/PSL_A0_y PSL_A0_y PSL_AH0 add def
N 0 240 M -42 0 D S
N 0 480 M -42 0 D S
N 0 720 M -42 0 D S
N 0 960 M -42 0 D S
N 0 1440 M -42 0 D S
N 0 1680 M -42 0 D S
N 0 1920 M -42 0 D S
N 0 2160 M -42 0 D S
N 0 2640 M -42 0 D S
N 0 2880 M -42 0 D S
N 0 3120 M -42 0 D S
N 0 3360 M -42 0 D S
N 0 3840 M -42 0 D S
N 0 4080 M -42 0 D S
N 0 4320 M -42 0 D S
N 0 4560 M -42 0 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
7800 0 T
25 W
N 0 4800 M 0 -4800 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 0 M 83 0 D S
N 0 1200 M 83 0 D S
N 0 2400 M 83 0 D S
N 0 3600 M 83 0 D S
N 0 4800 M 83 0 D S
N 0 240 M 42 0 D S
N 0 480 M 42 0 D S
N 0 720 M 42 0 D S
N 0 960 M 42 0 D S
N 0 1440 M 42 0 D S
N 0 1680 M 42 0 D S
N 0 1920 M 42 0 D S
N 0 2160 M 42 0 D S
N 0 2640 M 42 0 D S
N 0 2880 M 42 0 D S
N 0 3120 M 42 0 D S
N 0 3360 M 42 0 D S
N 0 3840 M 42 0 D S
N 0 4080 M 42 0 D S
N 0 4320 M 42 0 D S
N 0 4560 M 42 0 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
-7800 0 T
25 W
N 0 0 M 7800 0 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 0 M 0 -83 D S
N 1300 0 M 0 -83 D S
N 2600 0 M 0 -83 D S
N 3900 0 M 0 -83 D S
N 5200 0 M 0 -83 D S
N 6500 0 M 0 -83 D S
N 7800 0 M 0 -83 D S
/PSL_AH0 0
/MM {neg M} def
(”3) sh mx
(”2) sh mx
(”1) sh mx
(0) sh mx
(1) sh mx
(2) sh mx
(3) sh mx
def
/PSL_A0_y PSL_A0_y 83 add PSL_AH0 add def
0 PSL_A0_y MM
(”3) bc Z
1300 PSL_A0_y MM
(”2) bc Z
2600 PSL_A0_y MM
(”1) bc Z
3900 PSL_A0_y MM
(0) bc Z
5200 PSL_A0_y MM
(1) bc Z
6500 PSL_A0_y MM
(2) bc Z
7800 PSL_A0_y MM
(3) bc Z
N 260 0 M 0 -42 D S
N 520 0 M 0 -42 D S
N 780 0 M 0 -42 D S
N 1040 0 M 0 -42 D S
N 1560 0 M 0 -42 D S
N 1820 0 M 0 -42 D S
N 2080 0 M 0 -42 D S
N 2340 0 M 0 -42 D S
N 2860 0 M 0 -42 D S
N 3120 0 M 0 -42 D S
N 3380 0 M 0 -42 D S
N 3640 0 M 0 -42 D S
N 4160 0 M 0 -42 D S
N 4420 0 M 0 -42 D S
N 4680 0 M 0 -42 D S
N 4940 0 M 0 -42 D S
N 5460 0 M 0 -42 D S
N 5720 0 M 0 -42 D S
N 5980 0 M 0 -42 D S
N 6240 0 M 0 -42 D S
N 6760 0 M 0 -42 D S
N 7020 0 M 0 -42 D S
N 7280 0 M 0 -42 D S
N 7540 0 M 0 -42 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
0 4800 T
25 W
N 0 0 M 7800 0 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 0 M 0 83 D S
N 1300 0 M 0 83 D S
N 2600 0 M 0 83 D S
N 3900 0 M 0 83 D S
N 5200 0 M 0 83 D S
N 6500 0 M 0 83 D S
N 7800 0 M 0 83 D S
N 260 0 M 0 42 D S
N 520 0 M 0 42 D S
N 780 0 M 0 42 D S
N 1040 0 M 0 42 D S
N 1560 0 M 0 42 D S
N 1820 0 M 0 42 D S
N 2080 0 M 0 42 D S
N 2340 0 M 0 42 D S
N 2860 0 M 0 42 D S
N 3120 0 M 0 42 D S
N 3380 0 M 0 42 D S
N 3640 0 M 0 42 D S
N 4160 0 M 0 42 D S
N 4420 0 M 0 42 D S
N 4680 0 M 0 42 D S
N 4940 0 M 0 42 D S
N 5460 0 M 0 42 D S
N 5720 0 M 0 42 D S
N 5980 0 M 0 42 D S
N 6240 0 M 0 42 D S
N 6760 0 M 0 42 D S
N 7020 0 M 0 42 D S
N 7280 0 M 0 42 D S
N 7540 0 M 0 42 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
0 -4800 T
0 setlinecap
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
4 W
V
{0 0 1 C} FS
35 0 821 Sc
35 130 854 Sc
35 260 1085 Sc
35 390 929 Sc
35 520 936 Sc
35 650 971 Sc
35 780 1162 Sc
35 910 814 Sc
35 1040 827 Sc
35 1170 806 Sc
35 1300 884 Sc
35 1430 975 Sc
35 1560 960 Sc
35 1690 1031 Sc
35 1820 812 Sc
35 1950 1123 Sc
35 2080 884 Sc
35 2210 1129 Sc
35 2340 1041 Sc
35 2470 1178 Sc
35 2600 1184 Sc
35 2730 1086 Sc
35 2860 1139 Sc
35 2990 1218 Sc
35 3120 1053 Sc
35 3250 1150 Sc
35 3380 1040 Sc
35 3510 930 Sc
35 3640 1292 Sc
35 3770 1227 Sc
35 3900 1179 Sc
35 4030 1021 Sc
35 4160 1046 Sc
35 4290 1126 Sc
35 4420 1276 Sc
35 4550 2622 Sc
35 4680 3830 Sc
35 4810 3649 Sc
35 4940 3654 Sc
35 5070 3639 Sc
35 5200 3809 Sc
35 5330 3702 Sc
35 5460 3594 Sc
35 5590 3830 Sc
35 5720 3907 Sc
35 5850 3781 Sc
35 5980 3587 Sc
35 6110 3841 Sc
35 6240 3974 Sc
35 6370 3914 Sc
35 6500 3971 Sc
35 6630 4029 Sc
35 6760 3813 Sc
35 6890 4001 Sc
35 7020 3933 Sc
35 7150 3754 Sc
35 7280 4128 Sc
35 7410 4001 Sc
35 7540 4100 Sc
35 7670 4083 Sc
35 7800 3941 Sc
U
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt psxy -R-3/3/0/4 -JX6.5i/4i -O -K stepfit_lsq.txt -W2p
%@PROJ: xy -3.00000000 3.00000000 0.00000000 4.00000000 -3.000 3.000 0.000 4.000 +xy
%%BeginObject PSL_Layer_2
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
33 W
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
0 852 M
4420 368 D
260 2512 D
3120 259 D
S
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt psxy -R-3/3/0/4 -JX6.5i/4i -O -K stepfit_svd.txt -Sc2p -Gred
%@PROJ: xy -3.00000000 3.00000000 0.00000000 4.00000000 -3.000 3.000 0.000 4.000 +xy
%%BeginObject PSL_Layer_3
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
4 W
V
{1 0 0 C} FS
17 0 852 Sc
17 130 863 Sc
17 260 874 Sc
17 390 885 Sc
17 520 895 Sc
17 650 906 Sc
17 780 917 Sc
17 910 928 Sc
17 1040 939 Sc
17 1170 950 Sc
17 1300 960 Sc
17 1430 971 Sc
17 1560 982 Sc
17 1690 993 Sc
17 1820 1004 Sc
17 1950 1014 Sc
17 2080 1025 Sc
17 2210 1036 Sc
17 2340 1047 Sc
17 2470 1058 Sc
17 2600 1068 Sc
17 2730 1079 Sc
17 2860 1090 Sc
17 2990 1101 Sc
17 3120 1112 Sc
17 3250 1122 Sc
17 3380 1133 Sc
17 3510 1144 Sc
17 3640 1155 Sc
17 3770 1166 Sc
17 3900 1177 Sc
17 4030 1187 Sc
17 4160 1198 Sc
17 4290 1209 Sc
17 4420 1220 Sc
17 4550 2476 Sc
17 4680 3732 Sc
17 4810 3743 Sc
17 4940 3753 Sc
17 5070 3764 Sc
17 5200 3775 Sc
17 5330 3786 Sc
17 5460 3797 Sc
17 5590 3807 Sc
17 5720 3818 Sc
17 5850 3829 Sc
17 5980 3840 Sc
17 6110 3851 Sc
17 6240 3861 Sc
17 6370 3872 Sc
17 6500 3883 Sc
17 6630 3894 Sc
17 6760 3905 Sc
17 6890 3915 Sc
17 7020 3926 Sc
17 7150 3937 Sc
17 7280 3948 Sc
17 7410 3959 Sc
17 7540 3970 Sc
17 7670 3980 Sc
17 7800 3991 Sc
U
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt pstext -R-3/3/0/4 -JX6.5i/4i -O -K -F+f12p+jRB+cRB -Dj0.2i
%@PROJ: xy -3.00000000 3.00000000 0.00000000 4.00000000 -3.000 3.000 0.000 4.000 +xy
%%BeginObject PSL_Layer_4
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
7560 240 M PSL_font_encode 0 get 0 eq {Standard+_Encoding /Helvetica /Helvetica PSL_reencode PSL_font_encode 0 1 put} if
200 F0
V MU 0 0 M (Fit y\(x\) = a + b*x + c*H\(x-0.5\) + ) FP 200 F12 (e) FP 200 F0 (\(x\)) FP pathbbox N pop exch pop add U neg 0 G
(Fit y\(x\) = a + b*x + c*H\(x-0.5\) + ) Z
200 F12 (e) Z
200 F0 (\(x\)) Z
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt pslegend -R-3/3/0/4 -JX6.5i/4i -O -K -DjTL+w1.75i+jTL+o0.1i/0.1i -F+p
%@PROJ: xy -3.00000000 3.00000000 0.00000000 4.00000000 -3.000 3.000 0.000 4.000 +xy
%%BeginObject PSL_Layer_5
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
120 3887 T
25 W
O1
793 2100 1050 397 Sr
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt psxy -R0/6.5/0/4 -Jx1i -O -K -N -S @GMTAPI@-000001 --PROJ_LENGTH_UNIT=inch --GMT_HISTORY=false
%@PROJ: xy 0.00000000 6.50000000 0.00000000 4.00000000 0.000 6.500 0.000 4.000 +xy
%%BeginObject PSL_Layer_6
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
4 W
V
O1
33 W
90 187 617 S-
{1 0 0 C} FS
O0
4 W
17 187 397 Sc
{0 0 1 C} FS
35 187 177 Sc
U
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt pstext -R0/6.5/0/4 -Jx1i -O -K -N -F+f+j @GMTAPI@-000002 --GMT_HISTORY=false
%@PROJ: xy 0.00000000 6.50000000 0.00000000 4.00000000 0.000 6.500 0.000 4.000 +xy
%%BeginObject PSL_Layer_7
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
427 547 M PSL_font_encode 0 get 0 eq {Standard+_Encoding /Helvetica /Helvetica PSL_reencode PSL_font_encode 0 1 put} if
200 F0
(LSQFIT solution) bl Z
427 327 M (SVDFIT solution) bl Z
427 107 M (Input data) bl Z
%%EndObject
-120 -3887 T
%%EndObject
0 A
FQ
O0
0 6000 TM

% PostScript produced by:
%@GMT: gmt psxy -R0/2/-3/4 -JX6.5i/4i -O -Baf @sinusoiddata.txt -Sc0.15c -Gblue -K -Y5i
%@PROJ: xy 0.00000000 2.00000000 -3.00000000 4.00000000 0.000 2.000 -3.000 4.000 +xy
%%BeginObject PSL_Layer_8
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
25 W
2 setlinecap
N 0 4800 M 0 -4800 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 686 M -83 0 D S
N 0 2057 M -83 0 D S
N 0 3429 M -83 0 D S
N 0 4800 M -83 0 D S
/PSL_AH0 0
/MM {neg exch M} def
PSL_font_encode 0 get 0 eq {Standard+_Encoding /Helvetica /Helvetica PSL_reencode PSL_font_encode 0 1 put} if
200 F0
(”2) sw mx
(0) sw mx
(2) sw mx
(4) sw mx
def
/PSL_A0_y PSL_A0_y 83 add def
686 PSL_A0_y MM
(”2) mr Z
2057 PSL_A0_y MM
(0) mr Z
3429 PSL_A0_y MM
(2) mr Z
4800 PSL_A0_y MM
(4) mr Z
/PSL_A0_y PSL_A0_y PSL_AH0 add def
N 0 0 M -42 0 D S
N 0 1371 M -42 0 D S
N 0 2743 M -42 0 D S
N 0 4114 M -42 0 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
7800 0 T
25 W
N 0 4800 M 0 -4800 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 686 M 83 0 D S
N 0 2057 M 83 0 D S
N 0 3429 M 83 0 D S
N 0 4800 M 83 0 D S
/PSL_AH0 0
/MM {exch M} def
(”2) sw mx
(0) sw mx
(2) sw mx
(4) sw mx
def
/PSL_A0_y PSL_A0_y 83 add PSL_AH0 add def
686 PSL_A0_y MM
(”2) mr Z
2057 PSL_A0_y MM
(0) mr Z
3429 PSL_A0_y MM
(2) mr Z
4800 PSL_A0_y MM
(4) mr Z
N 0 0 M 42 0 D S
N 0 1371 M 42 0 D S
N 0 2743 M 42 0 D S
N 0 4114 M 42 0 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
-7800 0 T
25 W
N 0 0 M 7800 0 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 0 M 0 -83 D S
N 1950 0 M 0 -83 D S
N 3900 0 M 0 -83 D S
N 5850 0 M 0 -83 D S
N 7800 0 M 0 -83 D S
/PSL_AH0 0
/MM {neg M} def
(0.0) sh mx
(0.5) sh mx
(1.0) sh mx
(1.5) sh mx
(2.0) sh mx
def
/PSL_A0_y PSL_A0_y 83 add PSL_AH0 add def
0 PSL_A0_y MM
(0.0) bc Z
1950 PSL_A0_y MM
(0.5) bc Z
3900 PSL_A0_y MM
(1.0) bc Z
5850 PSL_A0_y MM
(1.5) bc Z
7800 PSL_A0_y MM
(2.0) bc Z
N 390 0 M 0 -42 D S
N 780 0 M 0 -42 D S
N 1170 0 M 0 -42 D S
N 1560 0 M 0 -42 D S
N 2340 0 M 0 -42 D S
N 2730 0 M 0 -42 D S
N 3120 0 M 0 -42 D S
N 3510 0 M 0 -42 D S
N 4290 0 M 0 -42 D S
N 4680 0 M 0 -42 D S
N 5070 0 M 0 -42 D S
N 5460 0 M 0 -42 D S
N 6240 0 M 0 -42 D S
N 6630 0 M 0 -42 D S
N 7020 0 M 0 -42 D S
N 7410 0 M 0 -42 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
0 4800 T
25 W
N 0 0 M 7800 0 D S
/PSL_A0_y 83 def
/PSL_A1_y 0 def
8 W
N 0 0 M 0 83 D S
N 1950 0 M 0 83 D S
N 3900 0 M 0 83 D S
N 5850 0 M 0 83 D S
N 7800 0 M 0 83 D S
/PSL_AH0 0
/MM {M} def
(0.0) sh mx
(0.5) sh mx
(1.0) sh mx
(1.5) sh mx
(2.0) sh mx
def
/PSL_A0_y PSL_A0_y 83 add def
0 PSL_A0_y MM
(0.0) bc Z
1950 PSL_A0_y MM
(0.5) bc Z
3900 PSL_A0_y MM
(1.0) bc Z
5850 PSL_A0_y MM
(1.5) bc Z
7800 PSL_A0_y MM
(2.0) bc Z
/PSL_A0_y PSL_A0_y PSL_AH0 add def
N 390 0 M 0 42 D S
N 780 0 M 0 42 D S
N 1170 0 M 0 42 D S
N 1560 0 M 0 42 D S
N 2340 0 M 0 42 D S
N 2730 0 M 0 42 D S
N 3120 0 M 0 42 D S
N 3510 0 M 0 42 D S
N 4290 0 M 0 42 D S
N 4680 0 M 0 42 D S
N 5070 0 M 0 42 D S
N 5460 0 M 0 42 D S
N 6240 0 M 0 42 D S
N 6630 0 M 0 42 D S
N 7020 0 M 0 42 D S
N 7410 0 M 0 42 D S
/PSL_LH 0 def /PSL_L_y PSL_A0_y PSL_A1_y mx def
0 -4800 T
0 setlinecap
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
4 W
V
{0 0 1 C} FS
35 0 2060 Sc
35 39 2379 Sc
35 78 2781 Sc
35 117 2661 Sc
35 156 3242 Sc
35 195 3458 Sc
35 234 3518 Sc
35 273 3919 Sc
35 312 3980 Sc
35 351 4264 Sc
35 390 4159 Sc
35 429 4172 Sc
35 468 4033 Sc
35 507 3920 Sc
35 546 3776 Sc
35 585 3560 Sc
35 624 3347 Sc
35 663 3124 Sc
35 702 2916 Sc
35 741 2676 Sc
35 780 2443 Sc
35 819 2302 Sc
35 858 2022 Sc
35 897 1760 Sc
35 936 1477 Sc
35 975 1657 Sc
35 1014 1368 Sc
35 1053 1372 Sc
35 1092 1488 Sc
35 1131 1423 Sc
35 1170 1512 Sc
35 1209 1753 Sc
35 1248 1840 Sc
35 1287 1858 Sc
35 1326 2425 Sc
35 1365 2668 Sc
35 1404 2756 Sc
35 1443 3008 Sc
35 1482 3240 Sc
35 1521 3458 Sc
35 1560 3791 Sc
35 1599 4101 Sc
35 1638 4096 Sc
35 1677 4301 Sc
35 1716 4157 Sc
35 1755 4037 Sc
35 1794 3923 Sc
35 1833 3734 Sc
35 1872 3892 Sc
35 1911 3515 Sc
35 1950 3541 Sc
35 1989 2637 Sc
35 2028 2853 Sc
35 2067 2626 Sc
35 2106 2171 Sc
35 2145 2046 Sc
35 2184 2010 Sc
35 2223 1741 Sc
35 2262 1267 Sc
35 2301 1459 Sc
35 2340 1465 Sc
35 2379 1537 Sc
35 2418 1378 Sc
35 2457 1533 Sc
35 2496 1504 Sc
35 2535 1846 Sc
35 2574 2019 Sc
35 2613 2369 Sc
35 2652 2704 Sc
35 2691 2571 Sc
35 2730 2770 Sc
35 2769 3433 Sc
35 2808 3465 Sc
35 2847 3564 Sc
35 2886 3876 Sc
35 2925 4073 Sc
35 2964 4205 Sc
35 3003 4025 Sc
35 3042 4232 Sc
35 3081 3947 Sc
35 3120 4192 Sc
35 3159 3771 Sc
35 3198 3546 Sc
35 3237 3366 Sc
35 3276 3343 Sc
35 3315 2741 Sc
35 3354 2635 Sc
35 3393 2451 Sc
35 3432 2395 Sc
35 3471 1921 Sc
35 3510 1653 Sc
35 3549 1699 Sc
35 3588 1565 Sc
35 3627 1139 Sc
35 3666 1433 Sc
35 3705 1525 Sc
35 3744 1218 Sc
35 3783 1618 Sc
35 3822 1963 Sc
35 3861 1675 Sc
35 3900 2278 Sc
35 3939 2371 Sc
35 3978 2693 Sc
35 4017 3113 Sc
35 4056 3097 Sc
35 4095 3180 Sc
35 4134 3698 Sc
35 4173 3868 Sc
35 4212 3848 Sc
35 4251 3978 Sc
35 4290 3969 Sc
35 4329 4071 Sc
35 4368 3952 Sc
35 4407 3983 Sc
35 4446 3637 Sc
35 4485 3719 Sc
35 4524 3665 Sc
35 4563 3263 Sc
35 4602 2840 Sc
35 4641 2725 Sc
35 4680 2550 Sc
35 4719 2209 Sc
35 4758 1938 Sc
35 4797 1827 Sc
35 4836 1715 Sc
35 4875 1532 Sc
35 4914 1223 Sc
35 4953 1298 Sc
35 4992 1357 Sc
35 5031 1597 Sc
35 5070 1412 Sc
35 5109 1656 Sc
35 5148 1810 Sc
35 5187 1949 Sc
35 5226 2153 Sc
35 5265 2578 Sc
35 5304 2680 Sc
35 5343 3246 Sc
35 5382 3317 Sc
35 5421 3406 Sc
35 5460 3698 Sc
35 5499 3735 Sc
35 5538 3933 Sc
35 5577 4149 Sc
35 5616 4211 Sc
35 5655 3968 Sc
35 5694 3668 Sc
35 5733 3986 Sc
35 5772 3744 Sc
35 5811 3799 Sc
35 5850 3331 Sc
35 5889 3409 Sc
35 5928 2786 Sc
35 5967 2657 Sc
35 6006 2263 Sc
35 6045 2281 Sc
35 6084 1878 Sc
35 6123 1709 Sc
35 6162 1473 Sc
35 6201 1709 Sc
35 6240 1759 Sc
35 6279 1491 Sc
35 6318 1198 Sc
35 6357 1646 Sc
35 6396 1651 Sc
35 6435 1945 Sc
35 6474 2118 Sc
35 6513 2150 Sc
35 6552 2617 Sc
35 6591 2775 Sc
35 6630 3045 Sc
35 6669 3208 Sc
35 6708 3272 Sc
35 6747 3772 Sc
35 6786 4181 Sc
35 6825 4249 Sc
35 6864 4318 Sc
35 6903 4055 Sc
35 6942 4099 Sc
35 6981 3943 Sc
35 7020 4088 Sc
35 7059 3936 Sc
35 7098 3595 Sc
35 7137 3277 Sc
35 7176 3170 Sc
35 7215 2919 Sc
35 7254 2425 Sc
35 7293 2492 Sc
35 7332 2129 Sc
35 7371 2017 Sc
35 7410 1747 Sc
35 7449 1796 Sc
35 7488 1558 Sc
35 7527 1606 Sc
35 7566 1201 Sc
35 7605 1402 Sc
35 7644 1481 Sc
35 7683 1777 Sc
35 7722 1724 Sc
35 7761 1965 Sc
35 7800 2281 Sc
U
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt psxy -R0/2/-3/4 -JX6.5i/4i -O -K sinusoidfit_lsq.txt -W2p -i0,2
%@PROJ: xy 0.00000000 2.00000000 -3.00000000 4.00000000 0.000 2.000 -3.000 4.000 +xy
%%BeginObject PSL_Layer_9
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
33 W
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
0 2170 M
39 241 D
117 761 D
39 235 D
39 212 D
39 181 D
39 144 D
39 102 D
39 56 D
39 7 D
39 -40 D
39 -87 D
39 -131 D
39 -170 D
39 -203 D
39 -229 D
156 -1004 D
39 -229 D
39 -202 D
39 -170 D
39 -130 D
39 -86 D
39 -40 D
39 8 D
39 57 D
39 102 D
39 144 D
39 182 D
39 212 D
39 236 D
117 761 D
39 241 D
39 221 D
39 192 D
39 157 D
39 116 D
39 71 D
39 24 D
39 -24 D
39 -72 D
39 -117 D
39 -158 D
39 -192 D
39 -222 D
39 -241 D
117 -761 D
39 -235 D
39 -212 D
39 -181 D
39 -144 D
39 -102 D
39 -55 D
39 -8 D
39 41 D
39 87 D
39 131 D
39 170 D
39 203 D
39 229 D
156 1003 D
39 229 D
39 202 D
39 170 D
39 130 D
39 87 D
39 40 D
39 -9 D
39 -56 D
39 -102 D
39 -145 D
39 -182 D
39 -212 D
39 -236 D
117 -760 D
39 -242 D
39 -220 D
39 -193 D
39 -157 D
39 -116 D
39 -71 D
39 -24 D
39 25 D
39 72 D
39 116 D
39 158 D
39 193 D
39 221 D
39 241 D
117 761 D
39 235 D
39 212 D
39 181 D
39 144 D
39 102 D
39 56 D
39 7 D
39 -40 D
39 -87 D
39 -131 D
39 -170 D
39 -203 D
39 -229 D
156 -1004 D
39 -229 D
39 -202 D
39 -170 D
39 -130 D
39 -86 D
39 -40 D
39 8 D
39 57 D
39 102 D
39 144 D
39 182 D
39 212 D
39 236 D
117 761 D
39 241 D
39 221 D
39 192 D
39 157 D
39 116 D
39 71 D
39 24 D
39 -24 D
39 -72 D
39 -117 D
39 -158 D
39 -192 D
39 -222 D
39 -241 D
117 -761 D
39 -235 D
39 -212 D
39 -181 D
39 -144 D
39 -102 D
39 -55 D
39 -8 D
39 41 D
39 87 D
39 131 D
39 170 D
39 203 D
39 229 D
156 1003 D
39 229 D
39 202 D
39 170 D
39 130 D
39 87 D
39 40 D
39 -9 D
39 -56 D
39 -102 D
39 -145 D
39 -182 D
39 -212 D
39 -236 D
117 -760 D
39 -242 D
39 -220 D
39 -193 D
39 -157 D
39 -116 D
39 -71 D
39 -24 D
39 25 D
39 72 D
39 116 D
39 158 D
39 193 D
39 221 D
S
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt psxy -R0/2/-3/4 -JX6.5i/4i -O -K sinusoidfit_svd.txt -Sc2p -Gred -i0,2
%@PROJ: xy 0.00000000 2.00000000 -3.00000000 4.00000000 0.000 2.000 -3.000 4.000 +xy
%%BeginObject PSL_Layer_10
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
4 W
V
{1 0 0 C} FS
17 0 2170 Sc
17 39 2411 Sc
17 78 2665 Sc
17 117 2921 Sc
17 156 3172 Sc
17 195 3407 Sc
17 234 3619 Sc
17 273 3800 Sc
17 312 3944 Sc
17 351 4046 Sc
17 390 4102 Sc
17 429 4109 Sc
17 468 4069 Sc
17 507 3982 Sc
17 546 3851 Sc
17 585 3681 Sc
17 624 3478 Sc
17 663 3249 Sc
17 702 3002 Sc
17 741 2747 Sc
17 780 2491 Sc
17 819 2245 Sc
17 858 2016 Sc
17 897 1814 Sc
17 936 1644 Sc
17 975 1514 Sc
17 1014 1428 Sc
17 1053 1388 Sc
17 1092 1396 Sc
17 1131 1453 Sc
17 1170 1555 Sc
17 1209 1699 Sc
17 1248 1881 Sc
17 1287 2093 Sc
17 1326 2329 Sc
17 1365 2580 Sc
17 1404 2836 Sc
17 1443 3090 Sc
17 1482 3331 Sc
17 1521 3552 Sc
17 1560 3744 Sc
17 1599 3901 Sc
17 1638 4017 Sc
17 1677 4088 Sc
17 1716 4112 Sc
17 1755 4088 Sc
17 1794 4016 Sc
17 1833 3899 Sc
17 1872 3741 Sc
17 1911 3549 Sc
17 1950 3327 Sc
17 1989 3086 Sc
17 2028 2832 Sc
17 2067 2576 Sc
17 2106 2325 Sc
17 2145 2090 Sc
17 2184 1878 Sc
17 2223 1697 Sc
17 2262 1553 Sc
17 2301 1451 Sc
17 2340 1396 Sc
17 2379 1388 Sc
17 2418 1429 Sc
17 2457 1516 Sc
17 2496 1647 Sc
17 2535 1817 Sc
17 2574 2020 Sc
17 2613 2249 Sc
17 2652 2495 Sc
17 2691 2751 Sc
17 2730 3006 Sc
17 2769 3252 Sc
17 2808 3481 Sc
17 2847 3683 Sc
17 2886 3853 Sc
17 2925 3983 Sc
17 2964 4070 Sc
17 3003 4110 Sc
17 3042 4101 Sc
17 3081 4045 Sc
17 3120 3943 Sc
17 3159 3798 Sc
17 3198 3616 Sc
17 3237 3404 Sc
17 3276 3168 Sc
17 3315 2918 Sc
17 3354 2661 Sc
17 3393 2408 Sc
17 3432 2166 Sc
17 3471 1946 Sc
17 3510 1753 Sc
17 3549 1596 Sc
17 3588 1480 Sc
17 3627 1409 Sc
17 3666 1385 Sc
17 3705 1410 Sc
17 3744 1482 Sc
17 3783 1598 Sc
17 3822 1756 Sc
17 3861 1949 Sc
17 3900 2170 Sc
17 3939 2411 Sc
17 3978 2665 Sc
17 4017 2921 Sc
17 4056 3172 Sc
17 4095 3407 Sc
17 4134 3619 Sc
17 4173 3800 Sc
17 4212 3944 Sc
17 4251 4046 Sc
17 4290 4102 Sc
17 4329 4109 Sc
17 4368 4069 Sc
17 4407 3982 Sc
17 4446 3851 Sc
17 4485 3681 Sc
17 4524 3478 Sc
17 4563 3249 Sc
17 4602 3002 Sc
17 4641 2747 Sc
17 4680 2491 Sc
17 4719 2245 Sc
17 4758 2016 Sc
17 4797 1814 Sc
17 4836 1644 Sc
17 4875 1514 Sc
17 4914 1428 Sc
17 4953 1388 Sc
17 4992 1396 Sc
17 5031 1453 Sc
17 5070 1555 Sc
17 5109 1699 Sc
17 5148 1881 Sc
17 5187 2093 Sc
17 5226 2329 Sc
17 5265 2580 Sc
17 5304 2836 Sc
17 5343 3090 Sc
17 5382 3331 Sc
17 5421 3552 Sc
17 5460 3744 Sc
17 5499 3901 Sc
17 5538 4017 Sc
17 5577 4088 Sc
17 5616 4112 Sc
17 5655 4088 Sc
17 5694 4016 Sc
17 5733 3899 Sc
17 5772 3741 Sc
17 5811 3549 Sc
17 5850 3327 Sc
17 5889 3086 Sc
17 5928 2832 Sc
17 5967 2576 Sc
17 6006 2325 Sc
17 6045 2090 Sc
17 6084 1878 Sc
17 6123 1697 Sc
17 6162 1553 Sc
17 6201 1451 Sc
17 6240 1396 Sc
17 6279 1388 Sc
17 6318 1429 Sc
17 6357 1516 Sc
17 6396 1647 Sc
17 6435 1817 Sc
17 6474 2020 Sc
17 6513 2249 Sc
17 6552 2495 Sc
17 6591 2751 Sc
17 6630 3006 Sc
17 6669 3252 Sc
17 6708 3481 Sc
17 6747 3683 Sc
17 6786 3853 Sc
17 6825 3983 Sc
17 6864 4070 Sc
17 6903 4110 Sc
17 6942 4101 Sc
17 6981 4045 Sc
17 7020 3943 Sc
17 7059 3798 Sc
17 7098 3616 Sc
17 7137 3404 Sc
17 7176 3168 Sc
17 7215 2918 Sc
17 7254 2661 Sc
17 7293 2408 Sc
17 7332 2166 Sc
17 7371 1946 Sc
17 7410 1753 Sc
17 7449 1596 Sc
17 7488 1480 Sc
17 7527 1409 Sc
17 7566 1385 Sc
17 7605 1410 Sc
17 7644 1482 Sc
17 7683 1598 Sc
17 7722 1756 Sc
17 7761 1949 Sc
17 7800 2170 Sc
U
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt pstext -R0/2/-3/4 -JX6.5i/4i -O -K -F+f12p+jRB+cRB -Dj0.2i
%@PROJ: xy 0.00000000 2.00000000 -3.00000000 4.00000000 0.000 2.000 -3.000 4.000 +xy
%%BeginObject PSL_Layer_11
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
clipsave
0 0 M
7800 0 D
0 4800 D
-7800 0 D
P
PSL_clip N
7560 240 M PSL_font_encode 0 get 0 eq {Standard+_Encoding /Helvetica /Helvetica PSL_reencode PSL_font_encode 0 1 put} if
200 F0
V MU 0 0 M (Fit y\(x\) = a + b*cos\(6) FP 200 F12 (p) FP 200 F0 (x - c\) + ) FP 200 F12 (e) FP 200 F0 (\(x\) = a + b) FP 0 -50 G 140 F0 (1) FP 0 50 G 200 F0 (*cos\(6) FP 200 F12 (p) FP 200 F0 (x\) + b) FP 0 -50 G 140 F0 (2) FP 0 50 G 200 F0 (*sin\(6) FP 200 F12 (p) FP 200 F0 (x\) + ) FP 200 F12 (e) FP 200 F0 (\(x\)) FP pathbbox N pop exch pop add U neg 0 G
(Fit y\(x\) = a + b*cos\(6) Z
200 F12 (p) Z
200 F0 (x - c\) + ) Z
200 F12 (e) Z
200 F0 (\(x\) = a + b) Z
0 -50 G 140 F0 (1) Z
0 50 G 200 F0 (*cos\(6) Z
200 F12 (p) Z
200 F0 (x\) + b) Z
0 -50 G 140 F0 (2) Z
0 50 G 200 F0 (*sin\(6) Z
200 F12 (p) Z
200 F0 (x\) + ) Z
200 F12 (e) Z
200 F0 (\(x\)) Z
PSL_cliprestore
%%EndObject
0 A
FQ
O0
0 0 TM

% PostScript produced by:
%@GMT: gmt psxy -R0/2/-3/4 -JX6.5i/4i -O -T
%@PROJ: xy 0.00000000 2.00000000 -3.00000000 4.00000000 0.000 2.000 -3.000 4.000 +xy
%%BeginObject PSL_Layer_12
0 setlinecap
0 setlinejoin
3.32551 setmiterlimit
%%EndObject

grestore
PSL_movie_completion /PSL_movie_completion {} def
%PSL_Begin_Trailer
%%PageTrailer
U
showpage

%%Trailer

end
%%EOF