clear;
clc;

sim('id_vgs',1);
ids = Id{1}.Values.Data(65:120);
vgs = Vgs{1}.Values.Data(65:120);
p = polyfit(vgs(25:56), ids(25:56), 1);

line = polyval(p, vgs);
vt_minus_vds_2 = -p(2)/p(1) - 0.05;

figure(1)
plot(vgs, ids);
hold on
plot(vgs, polyval(p, vgs));
hold off
grid
xlabel('Vgs')
ylabel('Ids')


    