clear; clc

for i = 0:5:100
    x = 273.15 + i;
    y = satprs (x );
    plot ( x, y, '*r','Markersize',10);
    hold on;
end
