Курсовая работа: Численное моделирование и анализ переходных процессов в электрической цепи

t1:=0.001;

t2:=0.007;

n:=100;

hx:=(t2-t1)/n;

a2:=2170;a1:=17.98;a0:=0.0007;

b2:=-880; b1:=36.38;b0:=-0.027;

c2:=-3515;c1:=62.485;c0:=-0.0917;

t:=t1;

S:=0;

fn:=sqr(a2*t1*t1+a1*t1+a0);

fk:=sqr(c2*t2*t2+c1*t2+c0);

repeat

if t<=0.003 then

f:=sqr(a2*t*t+a1*t+a0)

else if t<=0.005 then

f:=sqr(b2*t*t+b1*t+b0)

else f:=sqr(c2*t*t+c1*t+c0);

S:=S+f;

t:=t+hx;

until t>=t2;

S:=S-(fn+fk);

Int:=hx*(((fn+fk)/2)+S);

Q:=R4*Int;

writeln(' Int=',Int:2:8,' Q=',Q:2:7);

end.

Результат: Int =0.00007562

Q =0.0001422

5.3 Вычисление количества теплоты в пакете MathCAD