x=0:0.1:20*pi;
y=sin(x);
plot(x,y);
% save test.mat x y
print -dpng -r300 testplot.png

ls

[a, b] = urlread('http://qaweb2.neptune.uvic.ca/QaqcResultService?operation=4&sensorId=4211&dateFrom=2001-01-01T00:00:00.000Z&dateTo=2099-01-01T00:00:00.000Z')