<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">On 05/06/2015 07:41 AM, Andre Oliveira
wrote:<br>
</div>
<blockquote
cite="mid:907000342.1014472.1430912491034.JavaMail.yahoo@mail.yahoo.com"
type="cite">
<div style="color:#000; background-color:#fff; font-family:bookman
old style, new york, times, serif;font-size:16px">
<div id="yui_3_16_0_1_1430910871251_6719"><span
id="yui_3_16_0_1_1430910871251_6814">Prezados,</span></div>
<div id="yui_3_16_0_1_1430910871251_6931" dir="ltr"><span
id="yui_3_16_0_1_1430910871251_6814">meu interesse neste
momento é usar o R para uma apresentação com enfoque cálculo
numérico em uma escola. Fiz alguns testes aqui com</span>o
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! <br>
</div>
<div id="yui_3_16_0_1_1430910871251_6760">
<div id="yui_3_16_0_1_1430910871251_7014" dir="ltr"> <br>
</div>
<div style="" class="" id="yui_3_16_0_1_1430910871251_6931"
dir="ltr"> </div>
</div>
<div id="yui_3_16_0_1_1430910871251_6759">
<div id="yui_3_16_0_1_1430910871251_6758">
<br>
<br>
André Oliveira Souza. <br>
Graduação em Matemática, mestrado em estatística
aplicada.Instituto Federal de Educação, Ciência e Tecnologia
do Espirito Santo. IFES
<br>
<br>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
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.<br>
<br>
#include <stdio.h><br>
#include <math.h><br>
<br>
int main(void) {<br>
int i;<br>
float sum,j,a=1.0/3.0;<br>
<br>
printf ("Quantas vezes pretende somar? ");<br>
scanf ("%d",&i);<br>
<br>
sum = 0;<br>
<br>
for (j=0;j<=i;j++) {<br>
printf ("ITERAÇÃO %.0f SOMA %.10f MULTI %.10f ERRO
%g%%\n",j,sum,j*(1.0/3.0),(((j)*(a) - sum)/((j)*(a)))*100);<br>
sum = sum + a;<br>
}<br>
<br>
return(0);<br>
}<br>
<br>
ITERAÇÃO 999992 SOMA 333827.4375000000 MULTI 333330.6666666666 ERRO
-0.149026%<br>
ITERAÇÃO 999993 SOMA 333827.7812500000 MULTI 333331.0000000000 ERRO
-0.149035%<br>
ITERAÇÃO 999994 SOMA 333828.1250000000 MULTI 333331.3333333333 ERRO
-0.149035%<br>
ITERAÇÃO 999995 SOMA 333828.4687500000 MULTI 333331.6666666666 ERRO
-0.149035%<br>
ITERAÇÃO 999996 SOMA 333828.8125000000 MULTI 333332.0000000000 ERRO
-0.149044%<br>
ITERAÇÃO 999997 SOMA 333829.1562500000 MULTI 333332.3333333333 ERRO
-0.149044%<br>
ITERAÇÃO 999998 SOMA 333829.5000000000 MULTI 333332.6666666666 ERRO
-0.149044%<br>
ITERAÇÃO 999999 SOMA 333829.8437500000 MULTI 333333.0000000000 ERRO
-0.149053%<br>
ITERAÇÃO 1000000 SOMA 333830.1875000000 MULTI 333333.3333333333 ERRO
-0.149053%<br>
<br>
<br>
<pre class="moz-signature" cols="72">--
---
Zhu</pre>
</body>
</html>