Autor |
Wiadomość |
sandra
Doctor
Dołączył: 20 Sie 2009
Posty: 846
Przeczytał: 0 tematów
Pomógł: 8 razy Ostrzeżeń: 0/5
Płeć: Inżynierka
|
|
lab 2 inf |
|
Hej, czy ktoś mógłby mi podać treści zdań do zrobienia z drugiej laborki?
z góry dzięki 
Post został pochwalony 0 razy
|
|
Sob 16:55, 20 Kwi 2013 |
|
 |
|
 |
piotrek90
Magister Inżynier
Dołączył: 14 Wrz 2010
Posty: 116
Przeczytał: 0 tematów
Pomógł: 1 raz Ostrzeżeń: 0/5
Płeć: Inżynier
|
|
|
|
Przecież są w instrukcji podane na początku w ramce
W naszym przypadku mamy obliczyć metodą eulera/rungego dla różniczki o wzorze :
y'(x) = (y(x))^2-1)/(x^2 - 1)
Post został pochwalony 0 razy
|
|
Sob 18:21, 20 Kwi 2013 |
|
 |
Martyna_G
Magister Inżynier
Dołączył: 18 Gru 2009
Posty: 113
Przeczytał: 0 tematów
Pomógł: 1 raz Ostrzeżeń: 0/5
Płeć: Inżynierka
|
|
|
|
^
warunki początkowe:
y(2)=2
x nalezy do [2,4]
rozwiązanie dokładne:
y(x)=x
krok:
h=0,1
h=0,05
h=0,025
h=0,0125
Post został pochwalony 0 razy
|
|
Sob 21:42, 20 Kwi 2013 |
|
 |
Gość
|
|
|
Mógłby ktoś poratować skryptem? Zagapiłam się i zapomniałam o tym...
|
|
Nie 13:31, 21 Kwi 2013 |
|
 |
Gość
|
|
|
EuLeR:
 |  |
y = zeros(((b-a)/h)+1,1);
x = zeros(length(y),1);
y(1) = y0;
x(1) = a;
for i = 2: length(y)
x(i) = a+(i-1)*h;
y(i) = y(i-1) + h*ef(x(i-1),y(i-1));
end
|
RK2:
 |  |
y = zeros(((b-a)/h)+1,1);
x = zeros(length(y),1);
y(1) = y0;
x(1) = a;
for i = 2: length(y)
x(i) = a+(i-1)*h;
y(i) = y(i-1) + 0.5*h*ef(x(i-1),y(i-1)) + 0.5*h*ef(x(i-1)+h,y(i-1)+h*ef(x(i-1),y(i-1)));
end
|
oraz funkcja do równania:
 |  |
function [ fx ] = ef( x,yx )
fx = (yx*yx-1)/(x*x-1);
end
|
|
|
Nie 14:33, 21 Kwi 2013 |
|
 |
Gość
|
|
|
DZIĘKUJĘ
|
|
Nie 17:21, 21 Kwi 2013 |
|
 |
Gość
|
|
|
Janku, siedzę i patrzę w Twój kod i nie rozumiem. Próbowałam po swojemu, nie wykreśla mi NIC. Ratuj!
|
|
Pon 13:29, 22 Kwi 2013 |
|
 |
|
|
Nie możesz pisać nowych tematów Nie możesz odpowiadać w tematach Nie możesz zmieniać swoich postów Nie możesz usuwać swoich postów Nie możesz głosować w ankietach
|
|
|