szukanie zaawansowane
 [ Posty: 2 ] 
Autor Wiadomość
Mężczyzna
PostNapisane: 9 mar 2019, o 21:38 
Użytkownik

Posty: 4
Lokalizacja: Blachownia
Witam. Mam rozwiązać rówananie różniczkowe załączone na zdjęciu.
https://zapodaj.net/14d3679210cb4.png.html
Tylko nie wiem jak zaimplementować w kod Te wspolczynniki K2,K3,K4.
Kod wygląda tak:

Kod:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

%% Parameters
d = 0.01;
h = 0.1;
t0 = 0;
tEnd = 2/d;
N = (tEnd - t0)/h;

%% Inicialization
T = [t0:h:tEnd];
y = zeros(N+1, 1);
y(0) = d;


%% Runge-Kutta Method

for i = 1:N
    f =  y(i)^2 - y(i)^3;
   
    K1 = f;
    K2 =
    K3 =
    K4 =
   
   
    y(i+1) = y(i) + h/6*(K1 + 2K2 + 2K3 + K4);

end

%% Drawing diagram

plot(T, y, 'r-');
title("Euler Method (Explicit)");
xlabel("t");
ylabel("y");
grid on;



Dzięki za każdą pomoc !
Uniwersytet Wrocławski Instytut Matematyczny - rekrutacja 2019
Góra
Mężczyzna
PostNapisane: 9 mar 2019, o 22:36 
Użytkownik

Posty: 4787
Kod w Matlabie

Kod:
1
2
3
4
5
6
7
8
9
10
11
h=(b-a)/n;
y=y0;
for i = i:n
t= a+(i-1)*h;
K1 = feval(f, t,y);
K2 = feval(f, t+h/2, y+h*K1/2);
K3 = feval(f, t+h/2, y+K2/2);
K4 = feval (f, t+h, y+h*K3);
y = y+h*(K1+2*K2+2*K3+K4)/6;
t = t+h;
Góra
Utwórz nowy temat Odpowiedz w temacie  [ Posty: 2 ] 


 Zobacz podobne tematy
 Tytuł tematu   Autor   Odpowiedzi 
 Równanie różniczkowe - metoda Cauchy'ego  majkz  3
 Metoda rozdzielenia zmiennych - zadanie 2  hubertwojtowicz  1
 Równanie II stopnia - metoda przewidywania  panzam  1
 Równanie II rzędu - metoda przewidywań  zisary  2
 rownanie rozniczkowe 2 rzedu - metoda przewidywan  alina007  11
 
Atom [Regulamin Forum] [Instrukcja LaTeX-a] [Poradnik] [F.A.Q.] [Reklama] [Kontakt]
Copyright (C) Karpatka.pl