Euler method: Difference between revisions

Content added Content deleted
Line 1,735: Line 1,735:
k := 0.07:
k := 0.07:
TR := 20:
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 = 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 = 20); # step size = 5
Euler(diff(T(t), t) = -k*(T(t) - TR), T(0) = 100, t = 100, numsteps = 10); # step size 10</lang>
Euler(diff(T(t), t) = -k*(T(t) - TR), T(0) = 100, t = 100, numsteps = 10); # step size = 10</lang>
{{out}}
{{out}}
<pre>
<pre>
Line 1,760: Line 1,760:


# step size = 2
# step size = 2
printf("Step Size : %a\n", 2);
printf("Step Size = %a\n", 2);
EulerMethod(f, 0, 100, 100, 2);
EulerMethod(f, 0, 100, 100, 2);


# step size = 5
# step size = 5
printf("\nStep Size : %a\n", 5);
printf("\nStep Size = %a\n", 5);
EulerMethod(f, 0, 100, 100, 5);
EulerMethod(f, 0, 100, 100, 5);


# step size = 10
# step size = 10
printf("\nStep Size : %a\n", 10);
printf("\nStep Size = %a\n", 10);
EulerMethod(f, 0, 100, 100, 10);</lang>
EulerMethod(f, 0, 100, 100, 10);</lang>
{{out}}
{{out}}