Olá, preciso de ajuda com uma atividade de linguagem c. O programa precisa ler o valor de certo bem, o número de prestações no qual foi pago, e o valor das prestações. A partir destes valores precisa calcular a taxa de juros pela raiz da função f(j)=somatório-(valor do bem/valor da prestação). Somatório é 1-(1/(1+j))^n/1-(1/(1+j)) quando j>0, sendo j a taxa de juros e n numero de prestações. Para achar a raiz da função é necessário usa o método de muller. Porém, o juros da astronômico, e já mexi em tudo e não sai. Se alguém poder ajudar a achar o erro.
O código é:
(ignorem as outras funções potencia, valor absoluto, erro relativo e raiz quadrada, é que não pode usar o math.h. O importante esta em negrito)
#include <stdio.h>
#include <stdlib.h>
#define TOL 0.000000001
double potencia(double x, int n)
{
int k;
double pot=1;
for(k = 1; k <= n; k++)
{
pot = pot*x;
}
return pot;
}
double valor_absoluto(double x)
{
if(x<0)
{
x=-1*x;
}
return x;
}
double erro_relativo(double x, double y)
{
double erro=1, valor;
if (y!=0)
{
valor=(x-y)/y;
erro=valor_absoluto (valor);
}
else
{
erro=1;
}
return erro;
}
double raiz_quadrada(double x)
{
double esp=1, x0, x1;
x0=x;
while(esp>=TOL)
{
x1=x0-(((x0*x0)-x)/(2*x0));
esp=valor_absoluto(((x1-x0)/x1));
x0=x1;
}
return x1;
}
double f(double j,double vbem, int nprest, double prest)
{
double somatorio, funcao;
if (j==0)
{
somatorio=nprest;
funcao=somatorio-(vbem/prest);
}
else
{
somatorio=(1-(potencia(1/(1+j),nprest)))/(1-(1/(1+j)));
funcao=somatorio-(vbem/prest);
}
return funcao;
}
double juros (double vbem, int nprest, double prest)
{
double x0=0, x1, fx0, fx1, x2, fx2, fx3, h1, h2, q2, a, b, c, x3, esp, fx=1, resp, qlinha, qlinhalinha;
x1=0.3*(prest/(vbem-prest));
x2=prest/(vbem-prest);
fx0=f(x0,vbem,nprest,prest);
fx1=f(x1,vbem,nprest,prest);
fx2=f(x2,vbem,nprest,prest);
h1=x1-x0;
h2=x2-x1;
q2=h2/h1;
a=(q2*fx2)-(q2*(1+q2)*fx1)+((potencia(q2,2))*fx0);
b=((2*q2+1)*fx2)-((potencia((1+q2),2))*fx1)-((potencia(q2,2))*fx0);
c=(1+q2)*fx2;
qlinha=(-2*c)/(b+(raiz_quadrada(b*b-(4*a*c))));
qlinhalinha=(-2*c)/(b-(raiz_quadrada(b*b-(4*a*c))));
if (valor_absoluto(qlinha)>valor_absoluto(qlinhalinha))
{
x3=x2+(h1*qlinhalinha);
}
else
{
x3=x2+(h1*qlinha);
}
fx3=f(x3,vbem,nprest,prest);
esp=erro_relativo(x3, x2);
x0=x1;
x1=x2;
x2=x3;
resp=x3;
if (fx0==0)
{
fx=fx0;
resp=x0;
}
if (fx1==0)
{
fx=fx1;
resp=x1;
}
if (fx2==0)
{
fx=fx2;
resp=x2;
}
while (esp>=TOL && fx!=0)
{
fx0=f(x0,vbem,nprest,prest);
fx1=f(x1,vbem,nprest,prest);
fx2=f(x2,vbem,nprest,prest);
h1=x1-x0;
h2=x2-x1;
q2=h2/h1;
a=(q2*fx2)-(q2*(1+q2)*fx1)+((potencia(q2,2))*fx0);
b=((2*q2+1)*fx2)-((potencia((1+q2),2))*fx1)-((potencia(q2,2))*fx0);
c=(1+q2)*fx2;
qlinha=(-2*c)/(b+(raiz_quadrada(b*b-(4*a*c))));
qlinhalinha=(-2*c)/(b-(raiz_quadrada(b*b-(4*a*c))));
if (valor_absoluto(qlinha)>valor_absoluto(qlinhalinha))
{
x3=x2+(h1*qlinhalinha);
}
else
{
x3=x2+(h1*qlinha);
}
fx3=f(x3,vbem,nprest,prest);
esp=erro_relativo(x3, x2);
x0=x1;
x1=x2;
x2=x3;
resp=x3;
if (fx0==0)
{
fx=fx0;
resp=x0;
}
if (fx1==0)
{
fx=fx1;
resp=x1;
}
if (fx2==0)
{
fx=fx2;
resp=x2;
}
}
return resp;
}
int main()
{
int nprest=1;
double vbem=1,prest=1;
while (vbem>0)
{
printf("Digite valor do bem,numero de prestacoes e valor das prestacoes=");
scanf("%lf %lf %lf", &vbem, &nprest, &prest);
printf("juros= %lf\n", juros(vbem,nprest,prest));
}
system("PAUSE");
return 0;
}
