# Question:how to plot inside the for loop

## Question:how to plot inside the for loop

Maple 9.5

Below is my code:

> restart:
> q:=1.6*10^(-19):e0:=8.85*10^(-12):k:=1.38*10^(-23):t:=300:ni:=1.45*10^(16):eox:=3.9*e0:esi:=11.7*e0:
> tox:=2*10^(-9):tsi:=25*10^(-9):L:=1*10^(-6):cox:=(eox/tox):csi:=esi/tsi:a:=sqrt(q^2*ni/(2*k*t*esi)):V:=0:vgs:=2:u:=0.04:
> mm:=(cos((q*cox*(vgs-psi11)/(esi*k*t*a))*(exp(q*(V-psi11)/(2*k*t))))-((2*q*cox*(vgs-psi11)/(esi*k*t*a^2*(tsi)))*(exp(-q*(psi11-V)/(k*t))))):
> psi1:=solve(mm=0,psi11);
psi1 := 0.5434311697 - 0.003168514913 I
> mm:=(cos((q*cox*(vgs-psi22)/(esi*k*t*a))*(exp(q*(vds-psi22)/(2*k*t))))-((2*q*cox*(vgs-psi22)/(esi*k*t*a^2*(tsi)))*(exp(-q*(psi22-vds)/(k*t))))):
> qq2:=seq(solve(mm=0,psi22),vds=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2]):
> psi2:=(vector([qq2]));vds:=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2]:
psi2 := [0.5434311697 - 0.003168514913 I, 0.6398902333 - 0.003560933110 I,

0.7360983262 - 0.003949922699 I, 0.8586260376, 0.9544714953, 1.066786341,

1.145104060, 1.239782678, 1.333965608, 1.427671475, 1.521739902,

1.612717888 + 0.004670871971 I, 1.699606898 + 0.007295737134 I,

1.782802467 - 0.009645080463 I, 1.913366188, 1.969622161, 1.996631837,

1.999909512, 1.999998089, 1.999999960, 1.999999999]
> for i from 1 to 21 do
> Qgg[i] := Re( tsi*cox*L*(1/6*psi2[i]^3-1/2*psi2[i]^2-1/2*psi2[i]^2*vgs-1/2*psi2[i]^2*k*t/q+1/2*psi2[i]*vgs^2+psi2[i]*k*t*vgs/q+psi2[i]*vgs-8*psi2[i]*esi*k*t/(q*tsi*cox)-64*esi^2*k^2*t^2*ln(-8*esi*k*t-q*tsi*cox*vgs+q*tsi*cox*psi2[i])/(q^2*tsi^2*cox^2)-1/6*psi1^3+1/2*psi1^2+1/2*psi1^2*vgs+1/2*psi1^2*k*t/q-1/2*psi1*vgs^2-psi1*k*t*vgs/q-psi1*vgs+8*psi1*esi*k*t/(q*tsi*cox)+64*esi^2*k^2*t^2*ln(-8*esi*k*t-q*tsi*cox*vgs+q*tsi*cox*psi1)/(q^2*tsi^2*cox^2))/(k*t*(vg-psi2[i])/q+1/8*(vgs-psi2[i])^2+4*k^2*t^2*esi*ln(cox(vgs-psi2[i])+8*k*t*esi/(q*tsi))/(q^2*tsi*cox)-k*t*(vg-psi1)/q-1/8*(vgs-psi1)^2-4*k^2*t^2*esi*ln(cox(vgs-psi1)+8*k*t*esi/(q*tsi))/(q^2*tsi*cox)));
> end do:
> Qg:=simplify((Qgg));
/[ -15
Qg := TABLE\[1 = Float(infinity), 2 = -2.726896811 10 ,

-15 -15 -15
3 = -2.681258392 10 , 5 = -2.579829239 10 , 4 = -2.623431373 10 ,

-15 -15 -15
7 = -2.495397349 10 , 6 = -2.529657206 10 , 10 = -2.378078479 10 ,

-15 -15 -15
11 = -2.341729539 10 , 8 = -2.454901580 10 , 9 = -2.415751244 10 ,

-15 -15 -15
15 = -2.212746094 10 , 14 = -2.250797897 10 , 13 = -2.277803789 10 ,

-15 -15 -15
12 = -2.308132152 10 , 21 = -2.192572935 10 , 20 = -2.192572943 10 ,

-15 -15 -15
16 = -2.198907782 10 , 17 = -2.193213670 10 , 18 = -2.192589892 10 ,

-15]\
19 = -2.192573293 10 ]/

i have to plot Qg with respect to vds. the value of Qg are stored in a table which have mixed entries. also please suggest how can i set the x-axis of the plot in decimal notation?

Regards

﻿