1#!/usr/bin/env bash 2# Testing gmt trend1d 3 4ps=trend1.ps 5 6# Created a fake y(x) = (x/50)^2 + (x/60) + 4 + cos(2*pi*x/25) + noise and fit LS model 7#gmt math -T10/110/1 T 50 DIV 2 POW 2 MUL T 60 DIV ADD 4 ADD 0 0.25 NRAND ADD T 25 DIV 2 MUL PI MUL COS 2 MUL 2 ADD ADD = data1.txt 8gmt trend1d data1.txt -Np2,F1+o0+l25 -Fxm > model.txt 9gmt psxy -R0/120/0/20 -JX6i/2.5i -Baf data1.txt -Sc0.1c -Gblack -P -K -Xc > $ps 10gmt psxy -R -J -O -K -W1p,blue model.txt >> $ps 11# Created a fake y(x) = 4*(x/50)^3 + 3 + noise and fit LS model 12#gmt math -T-50/50/1 T 50 DIV 3 POW 4 MUL 3 ADD 0 0.25 NRAND ADD = data2.txt 13gmt trend1d data2.txt -NP0,P3 -Fxm > model.txt 14gmt psxy -R-50/50/0/6 -J -Baf data2.txt -Sc0.1c -Gblack -O -K -Y3.25i >> $ps 15gmt psxy -R -J -O -K -W1p,blue model.txt >> $ps 16# Created a fake y(x) = x + 3 + noise and fit LS model using old syntax 17#gmt math -T-50/50/1 T 50 DIV 3 ADD 0 0.25 NRAND ADD = data3.txt 18gmt trend1d data3.txt -N2 -Fxm > model.txt 19gmt psxy -R-50/50/0/6 -J -Baf data3.txt -Sc0.1c -Gblack -O -K -Y3.25i >> $ps 20gmt psxy -R -J -O -W1p,blue model.txt >> $ps 21