[R-br] Cálulo numérico

[Rodolfo Timoteo da Silva]

On 05/06/2015 07:41 AM, Andre Oliveira wrote:
> Prezados,
> meu interesse neste momento é usar o R para uma apresentação com 
> enfoque cálculo numérico em uma escola. Fiz alguns testes aqui como 
> pacote library(animation), com a função optim e seus argumentos e 
> contruí um algoritmo de Newton Raphso em caso univariado. Gostaria de 
> saber se alguém do grupo tem experiência com R em cálculo numérico e 
> se tem algum material para optimização e algoritmos que possa 
> compartilhar? No caso pode ser uni e multivariado!
>
>
>
> André Oliveira Souza.
> Graduação em Matemática, mestrado em estatística aplicada.Instituto 
> Federal de Educação, Ciência e Tecnologia do Espirito Santo.  IFES
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Uma pergunta Sr. Andre: Como o R trabalha com propagação de erros do 
tipo abaixo? Você tem isso em mente em seu projeto de aula? Atualmente 
estou fazendo uma I.C. em relação a isso, mas utilizando C/C++/Matlab 
comparativamente. Futuramente estarei me dedicando a estatística já que 
é mais ligada ao meu curso em Controle Estatístico da Qualidade então é 
interessante ter uma noção de como o R trabalha com os problemas 
relacionados a IEEE 754.

#include <stdio.h>
#include <math.h>

int main(void) {
     int i;
     float sum,j,a=1.0/3.0;

     printf ("Quantas vezes pretende somar? ");
     scanf ("%d",&i);

     sum = 0;

     for (j=0;j<=i;j++) {
         printf ("ITERAÇÃO %.0f SOMA %.10f MULTI %.10f ERRO 
%g%%\n",j,sum,j*(1.0/3.0),(((j)*(a) - sum)/((j)*(a)))*100);
         sum = sum + a;
     }

     return(0);
}

ITERAÇÃO 999992 SOMA 333827.4375000000 MULTI 333330.6666666666 ERRO 
-0.149026%
ITERAÇÃO 999993 SOMA 333827.7812500000 MULTI 333331.0000000000 ERRO 
-0.149035%
ITERAÇÃO 999994 SOMA 333828.1250000000 MULTI 333331.3333333333 ERRO 
-0.149035%
ITERAÇÃO 999995 SOMA 333828.4687500000 MULTI 333331.6666666666 ERRO 
-0.149035%
ITERAÇÃO 999996 SOMA 333828.8125000000 MULTI 333332.0000000000 ERRO 
-0.149044%
ITERAÇÃO 999997 SOMA 333829.1562500000 MULTI 333332.3333333333 ERRO 
-0.149044%
ITERAÇÃO 999998 SOMA 333829.5000000000 MULTI 333332.6666666666 ERRO 
-0.149044%
ITERAÇÃO 999999 SOMA 333829.8437500000 MULTI 333333.0000000000 ERRO 
-0.149053%
ITERAÇÃO 1000000 SOMA 333830.1875000000 MULTI 333333.3333333333 ERRO 
-0.149053%


-- 

---
Zhu

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150506/44105aaa/attachment.html>