rlc.inp

来自「Differential equation in TKSL」· INP 代码 · 共 8 行

var u1, u2, ur, uc, uL, i; 
const dt=1e-3, eps=1e-20, tmax=1e-1, c=1e-6, L=1, r=100, omega=0.1; 
system 
u1 = 1*sin(omega*t);
ur=R*i;
uc'= 1/c*i &0; {z韘k醤o ze vztahu i(t) = C.uC'(t) pro proud kondenz醫orem}
uL = u1-R*i-uc; {z韘k醤o z 2. Kirchhofova z醟ona}
i'= uL/L &0; {z韘k醤o ze vztahu uL(t) = L.i'(t) pro nap靦