# SECTION 1.4# Example 1f := (x,y) -> evalf(x*sqrt(y)); #Define the function f

NiM+JSJmR2orNiQlInhHJSJ5RzYiNiQlKW9wZXJhdG9yRyUmYXJyb3dHRiktJSZldmFsZkc2IyomOSQiIiItJSVzcXJ0RzYjOSVGMkYpRilGKTYkIiIhRjg=

h := 0.1;

NiM+JSJoRyQiIiIhIiI=

x0 := 1.; y0 := 4.;

NiM+JSN4MEckIiIiIiIhNiM+JSN5MEckIiIlIiIh

x1 := x0 + h;

NiM+JSN4MUckIiM2ISIi

y1 := y0 + h*f(x0,y0);

NiM+JSN5MUckIisrKysrVSEiKg==

x2 := x1 + h;

NiM+JSN4MkckIiM3ISIi

y2 := y1 + h*f(x1,y1);

NiM+JSN5MkckIis8SFZEVyEiKg==

# Use a repetition statement ( a 'for' loop): x := x0: y := y0: #initialize x and yprint('n','x','y');print(0,x,y);

NiUlIm5HJSJ4RyUieUc=NiUiIiEkIiIiRiMkIiIlRiM=

# Review the syntax of the "for" statement?for for n from 1 to 5 doy := y + h*f(x,y):x := x + h:print(n,x,y);end do:

NiUiIiIkIiM2ISIiJCIrKysrK1UhIio=NiUiIiMkIiM3ISIiJCIrPEhWRFchIio=NiUiIiQkIiM4ISIiJCIrc00oeW4lISIqNiUiIiUkIiM5ISIiJCIrcUQvZlwhIio=NiUiIiYkIiM6ISIiJCIrIUcyM0YmISIq

# Example 2f := (x,y) -> evalf(y);

NiM+JSJmR2orNiQlInhHJSJ5RzYiNiQlKW9wZXJhdG9yRyUmYXJyb3dHRiktJSZldmFsZkc2IzklRilGKUYpNiQiIiFGMg==

x0 := 0; y0 := 1.;

NiM+JSN4MEciIiE=NiM+JSN5MEckIiIiIiIh

Numbersteps := [1,2,4,8,16];

NiM+JSxOdW1iZXJzdGVwc0c3JyIiIiIiIyIiJSIiKSIjOw==

print('N','h','Approximation');

NiUlIk5HJSJoRyUuQXBwcm94aW1hdGlvbkc=

for k from 1 to 5 doN := Numbersteps[k]:h := evalf(1/N):x := x0: y := y0:for n from 1 to N doy := y + h*f(x,y):x := x + h:end do:print(N,h,y); end do:

NiUiIiIkRiMiIiEkIiIjRiU=NiUiIiMkIisrKysrXSEjNSQiKysrK11BISIqNiUiIiUkIisrKysrRCEjNSQiK11pU1RDISIqNiUiIikkIisrKytdNyEjNSQiKzlYeWxEISIqNiUiIzskIisrKytdaSEjNiQiKypcR3pqIyEiKg==