witam, mam takie zadanko, i chciałabym je jakoś ugryźć i jak zaczynam pisac w maple program to juz na samym poczatku wyskakują mi błędy, co robie nie tak?
Predkosc obiektu w spadku swobodnym blisko powierzchni Ziemi opisuje rownanie
dv/dt=-g
gdzie v-predkosc, g- przyspieszenie ziemskie
Napisz program ktory uzywa metode Eulera do obliczenia tegoz rownania rozniczkowego (tzn oblicza funkcje v(t)). Dla uproszczenia zakladamy, ze v(0)=0 i obliczamy rownanie na przedziale czasowym od t=0 do t=10. Wykonaj obliczenia dla kilku roznych time step ( tzn ilosc iteracji) i porownaj swoje rozwiazanie z rozwiazaniem analitycznym ( takim obliczonym na papierze
> dv/dt = -g;
dv
-- = -g
dt
> f := proc (v, t) options operator, arrow; (-1)*9.81 end proc;
(v, t) -> (-1) 9.81
> a := 0.; b := 10; n := 100; h := (b-a)/n;
0.
10
100
0.1000000000
> x := array(0 .. n, [seq(a+i*h, i = 0 .. n)]);
> y := array(0 .. n); y[0] := 0;
0
> for i from 0 to n-1 do;
Error, invalid loop
Typesetting:-mambiguous(for i from 0 to n - 1 do,
Typesetting:-merror("invalid loop"))
> y[i+1] := y[i]+h*f(x[i], y[i]);
Error, unevaluated array index in assignment
>