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@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