/* Programa riemann.calc Assunto: Somas de Riemann Descrição: cálculo de integrais aprox por somas de Riemann palavras chave: Riemann, integral, varredura por Tarcisio Praciano Pereira - C para matemáticos Sobral, Dezembro de 2005 - UeVA */ define f1(x) { return( re(power(x,2)*cos(2*x))); } define f2(x) { return( re(power(x,2)*sin(2*x))); } define f3(x) { return( re(x*cos(2*x))); } define f4(x) { return( re(x*sin(2*x))); } define entrada_float(msg, x) { printf("%s %f ",msg, x); scanf("%f", x); return(x); } define Riemann(inicio, fim, delta) { local soma1=0, soma2=0, soma3=0, soma4=0,x=inicio; local sp1 = "I1= ", sp2 = "I2= ",sp3 = "I3= ",sp4 = "I4= "; local I1,I2,I3,I4,I5; while (x <= fim) { soma1 = soma1 + f1(x); soma2 = soma2 + f2(x); soma3 = soma3 + f3(x); soma4 = soma4 + f4(x); x = x + delta; } I1 = (cos(12)-cos(2))*soma1*delta; ## power(x,2)*cos(2*x) I2 = (sin(12)+sin(2))*soma2*delta; ## power(x,2)*sin(2*x); I3 = (- 10*cos(12)- 4*cos(2)+ sin(12)+ sin(2) )*soma3*delta; ## x*cos(2*x) I4 = ( - 10*sin(12)+ 4*sin(2)- cos(12)+ cos(2) )*soma4*delta; ## x*sin(2*x) print I1 + I2 + I3 + I4 ; return(0); } Riemann(-4, 4, 0.000001) quit;