Euler method: Difference between revisions

=={{header|Run BASIC}}==

<lang rinbasic>x = euler(-0.07,-20, 100, 0, 100, 2)

x = euler-0.07,-20, 100, 0, 100, 5)

x = euler(-0.07,-20, 100, 0, 100, 10)

end

FUNCTION euler(da,db, y, a, b, s)

print "===== da:";da;" db:";db;" y:";y;" a:";a;" b:";b;" s:";s;" ==================="

t = a

WHILE t <= b

PRINT t;chr$(9);y

y = y + s * (da * (y + db))

t = t + s

WEND

END FUNCTION</lang>

<pre>===== da:-0.07 db:-20 y:100 a:0 b:100 s:2 ===================

0 100

2 88.8

4 79.168

6 70.88448

8 63.7606528

10 57.6341614

12 52.3653788

14 47.8342258

......

===== da:-0.07 db:-20 y:100 a:0 b:100 s:10 ===================

0 100

10 44.0

20 27.2

30 22.16

40 20.648

50 20.1944

60 20.05832

70 20.017496

80 20.0052488</pre>