t=[30:5:100];
p=polyfit(t,R,2);
x1=0:.5:100;
y1=polyval(p,x1);
plot(t,R,'o');
axis([0 100 0 1500]);
xlabel('I(1℃)');
ylabel('Rp(1Ω)');
hold on;
plot(x1,y1);
grid on;
測試成功