lab8.mw

> x := 45;
 

(Typesetting:-mprintslash)([x := 45], [45]) 

> N := 10;
 

(Typesetting:-mprintslash)([N := 10], [10]) 

> y := x*(Pi/180);
 

(Typesetting:-mprintslash)([y := 1/4*Pi], [1/4*Pi]) 

> sum(sin(i*y), i=1..N);
 

1+1/2*2^(1/2) 

> evalf(%);
 

1.707106781 

> sumsines(x,N);
 

sumsines(45, 10) 

> sumsines:=proc(x,N):
sum(sin(i*x), i=1..N):
evalf(%):
end proc;
  
 

(Typesetting:-mprintslash)([sumsines := proc (x, N) sum(sin(i*x), i = 1 .. N); evalf(%) end proc], [proc (x, N) sum(sin(i*x), i = 1 .. N); evalf(%) end proc]) 

> sumsines(1,1);
 

.8414709848 

> for N from 1 to 10 do
print( sumsines(1,N))
od;
 

.8414709848 

1.750768412 

1.891888420 

1.135085925 

.1761616503 

-.1032538479 

.5537327508 

1.543090997 

1.955209482 

1.411188371 

>