Euler method: Difference between revisions

Line 1,735:
k := 0.07:
TR := 20:
Euler(diff(T(t), t) = -k*(T(t) - TR), T(0) = 100, t = 100, numsteps = 50); # step size = 2
Euler(diff(T(t), t) = -k*(T(t) - TR), T(0) = 100, t = 100, numsteps = 20); # step size = 5
Euler(diff(T(t), t) = -k*(T(t) - TR), T(0) = 100, t = 100, numsteps = 10); # step size = 10</lang>
{{out}}
<pre>
Line 1,760:
 
# step size = 2
printf("Step Size := %a\n", 2);
EulerMethod(f, 0, 100, 100, 2);
 
# step size = 5
printf("\nStep Size := %a\n", 5);
EulerMethod(f, 0, 100, 100, 5);
 
# step size = 10
printf("\nStep Size := %a\n", 10);
EulerMethod(f, 0, 100, 100, 10);</lang>
{{out}}
Anonymous user