?????????????? ?????? ? ????????? ?????????? ??????? ?????

Ucg(1)=Uc0; iLg(1)=iL0; t(1)=0.
do n=2, Nt ! Задание цикла по n
! Расчетные формулы
iL=(E-Uc0+L/dt*iL0)/(R1+dt/C+L/dt+R2)
Uc=dt/C*iL+Uc0
! присвоение новых начальных значений для напряжения и тока
iL0=iL; Uc0=Uc;
! формирование массивов тока, напряжения и времени для вывода на график
iLg(n)=100 * iL; Ucg(n)=Uc; t(n)=t(n-1)+dt
end do ! конец цикла по n
! Организация вывода результатов на график
call metafl('XWIN') ; !'XWIN' - вывод в окно, 'WMF' — вывод в графический файл
call disini ; ! инициализация библиотеки DISLIN
call titlin ( '100*iL', 2); ! вывод заголовка
call qplot ( t, iLg, Nt ); ! отображение графика
call disini
call titlin ( 'Uc', 1); ! вывод 2-го заголовка
call qplot ( t, Ucg, Nt ); ! отображение 2-го графика
call disfin ! завершение работы с DISLIN
end program rc ! завершение программы
В результате, после расчета появятся окна с графиками. Сначала появится
окно с графиком тока в индуктивности, затем, после нажатия на
правую кнопку мыши, появится окно с графиком напряжения на емкости.
Вид последнего окна будет следующим:
Вид окна DISLIN с графиком напряжения на емкости для варианта
расчета с шагом 0,1мс.
32