// As excel file of the same name on red usb
// Uses first go at probe with JRM amplifiers and Hamamatsu detectors
// At each temperature, signal from each channel and ratio
Data= [... 
-90	4.379	0.6338	0.144736241
-75	4.295	0.7312	0.17024447
-50	4.158	0.8922	0.214574315
-25	4.022	1.045	0.259820985
0	3.884	1.184	0.304840371
25	3.73	1.305	0.349865952
50	3.604	1.411	0.391509434
75	3.482	1.5	0.430786904
100	3.358	1.57	0.467540203
125	3.236	1.622	0.501236094
150	3.116	1.659	0.53241335
175	2.999	1.683	0.561187062
200	2.877	1.69	0.587417449
225	2.763	1.69	0.611653999
250	2.65	1.68	0.633962264
275	2.536	1.66	0.654574132
300	2.419	1.629	0.673418768
325	2.29	1.582	0.690829694
350	2.129	1.507	0.707844058
375	1.901	1.376	0.723829563
400	1.571	1.163	0.740292807
425	1.154	0.8719	0.755545927
450	0.7363	0.5683	0.771832134
475	0.417	0.3286	0.788009592
500	0.2203	0.1774	0.805265547
525	0.1186	0.0979	0.825463744
550	0.0616	0.05257	0.853409091
575	0.03329	0.02946	0.884950436
600	0.0187	0.01731	0.925668449
625	0.01084	0.01062	0.979704797
650	0.00651	0.006841	1.050844854
660	0.005324	0.005787	1.086964688
];
x=Data(:,1);
y1=Data(:,2);
y2=Data(:,3);
y3=Data(:,4);

clf()
f=gcf();
f.color_map = graycolormap(32)
y=[y1 y2];
linestyle=[1 1];
fs=3;
plot2d("nl",x,y,linestyle);
hl=legend("660 nm","630 nm",5,%f);
hl.font_size=fs;
xlabel("Temperature / °C","fontsize",fs)
ylabel("Signal / V","fontsize",fs)
a=gca();
a.font_size=fs;
a.children(1).links(1).line_style=2;//drill down to set dash line

axis_color = "gray60"
linestyle=color(axis_color);
y=y3;
h2=newaxes();
plot2d("nn",x,y,linestyle)
ylabel("Ratio","fontsize",fs,"color",axis_color)
h2.filled="off";
h2.axes_visible(1)="off";
h2.y_location="right";
h2.font_color=color(axis_color)
h2.foreground=color(axis_color)
h2.font_size=fs;
