plot(VGS,ID1,VGS,ID2);
xlabel('V_(G_S)(V)');
ylabel('I_D(mA)');
legend('สมการเส้นโหลด','สมการ Shockley')

% กำหนดตัวแปร VGS และ ID เป็นตัวแปรของระบบสมการ
syms VGS ID;

% กำหนดค่าคงที่ IDSS, VP, VG และ RS
IDSS = 8*10^(-3);
VP = -4;
RS = 500;

% กำหนดสมการ (1) และ (2) ไว้ในตัวแปร eqn1 และ eqn2 สังเกตเครื่องหมาย == ด้วย
eqn1 = VGS == -ID*RS;
eqn2 = ID == IDSS*(1-VGS/VP)^2;

% แก้สมการ (1) และ (2) เพื่อหาค่า VGSQ และ IDQ ด้วยคำสั่ง vpasolve เก็บผลลัพธ์ไว้ในตัวแปร sol
sol = vpasolve([eqn1 eqn2]);
VGSQ = sol.VGS
IDQ = sol.ID

% สังเกตว่าค่าของ VGSQ และ IDQ จะมีอย่างละ 2 ค่า ซึ่งมีเพียงคู่เดียวที่ถูกต้องและใช้งานได้สำหรับ FET นั้น ๆ
% ทำการตรวจสอบค่าที่ได้และเลือกคู่ค่า VGSQ และ IDQ ที่สอดคล้องกับจุดตัดกราฟในข้อ 3 เพื่อจะใช้ในข้อ 4 ต่อไป
% ค่าของ VGSQ มีสองค่า เราเลือกได้โดยใช้ VGSQ(1) หรือ VSGQ(2) และเป็นเช่นเดียวกันสำหรับ IDQ
% ให้นักศึกษาเพิ่ม code ตรงนี้เพื่อแสดงค่าจุด Q-point เพียงคู่เดียวที่ถูกต้อง
% โดย print แสดงใน Command Window ให้เหมาะสมเช่น
% IDQ มีค่าเท่ากับ xxx.xxxx
% VGSQ มีค่าเท่ากับ yyy.yyyy



หาจุดตัด
VGS = -4:0.1:0;
if VGSQ(1,1) < -4
     VGSQ2=VGSQ(2,1)
     IDQ2=IDQ(2,1)
end
if VGSQ(2,1) < -4
     VGSQ2=VGSQ(1,1)
     IDQ2=IDQ(1,1)
end
plot(VGS,ID1,VGS,ID2,VGSQ2,IDQ2,'go');
xlabel('V_(G_S)(V)');
ylabel('I_D(mA)');
legend('สมการเส้นโหลด','สมการ Shockley')