polecenie MN lab 5

Dane do zadania nr 5:

E = 5 V

R =50 Ohm

L = 0.1 H

C = 1e-5 F

x0 = [0;0]

Podpowiedzi:

- funkcja euler:

     function [x, xi] =euler(x0, fun, t, h)

            for i=1:length(t)
            x=x0+fun(x0)*h;
            x0=x;
            xi(:,i)=x;

            ...

Jeżeli przedział "t" dzielimy na mniejsze podprzedziały o szerokości "h" to możemy wykorzystać polecenie length - liczba podprzedziałów, na które dzielimy cały przedział t.

Sprawozdanie:

- pliki funkcja.m (plik wyznaczający wartość funkcji w każdym kolejnym kroku), euler.m (plik z procedurą eulera), driver.m (w pliku mają być dane do zadania, wywołanie funkcji euler oraz uruchomienie rysowania wykresów),

- printscreen z wszystkimi sześcioma wykresami,

- komentarz: czy uzyskane wyniki dla dwóch różnych kroków całkowania są poprawne czy też nie.