Здравствуйте, уважаемые форумчане.
Возникла необходимость выполнить расчет цепи при помощи диф.уравнений. Необходимо использовать фунцию rfixed.
Написал небольшой расчет с диффиренциальными уравнениями. Вроде бы все написал верно, но тем не менее расчет не происходит. Начальные условия указал, шаг, диапазон, исходное уравнение - все есть.

Может ли кто-нибудь найти ошибку в выложенном мною файле?

Нажмите для просмотра прикрепленного файла

Сейчас дошел до следующей проблемы. Уравнения упростил. Необходимо добиться, чтобы маткад посчитал диф.уравнения. Пробовал применять кроме функции rkfixed функцию Bulstoer, но считает очень долго и требует уменьшить шаг расчета. Уменьшаю до максимума, но все равно не получаю результата.
Может кто-нибудь взглянуть - в чем ошибка и как их можно посчитать?
Нажмите для просмотра прикрепленного файла

Рискну предположить, что проблема в исходных данных - их сочетание вызывает расчетные трудности.

1. Судя по е(t)=37sin(11t) период колебаний, навязываемых источником - 0,57 секунды.
2. С другой стороны, величины L и C Такие, что резонансная частота должна быть 0,3 МГц при периоде 3*10е-6 секунды.
3. Однако сопротивление R1 = 6,25 кОм - процесс для собственных колебаний явно апериодический.

То есть для расчета хотя бы трех колебаний на частоте е(t) надо брать интервал времени 1,5...2 секунды.
Но из-за переходного процесса в LC цепи вроде бы надо брать интервал 10е-7 секунды.

Формулы посмотрел на бегу - похоже на правду. И второй вариант - не с контурными токами, а с общим током и напряжением на C.

Прикрепляю свои версии файлов с результатими расчета для предложенных данных и еще дополнительно со своими L, C, R1.

Нажмите для просмотра прикрепленного файла
Нажмите для просмотра прикрепленного файла
Нажмите для просмотра прикрепленного файла

Спасибо, посмотрел. Проблему решил иначе. Прочитал про функцию BDF. Она смогла посчитать значения емкости и индуктивности с номиналами нано и микро.