0

R を使用して複雑な方程式を解き始めました。方程式を生成した後、 を使用して解こうとしましたRyacas。残念ながら、結果をRyacas返す代わりに、次を返します。

CommandLine(1) : 最大評価スタック深度に達しました。必要に応じて、MaxEvalDepth を使用してスタック サイズを増やしてください。
CommandLine(1) : 最大評価スタック深度に達しました。必要に応じて、MaxEvalDepth を使用してスタック サイズを増やしてください。

スタックサイズを増やす方法を教えてくださいRyacas。私はいろいろな方法でそれをやろうとしましたが、私にRyacas与えられたアドバイスをどのように利用するか本当にわかりません.

===== 編集 =======

これが、解きたい方程式を生成するためのコードです。

#define net and gross values
net=10000
gross=12563.49

#construct an array for cash flows
flows=matrix(nrow=1, ncol=60)

#populate the array with cash flows
flows[c(1:60)]=c(-297.21)

#generate the equation
#flows
eq1=NULL
for (i in 1:60) {
  eq1=paste(eq1," + ", toString(flows[i]),"/((1 + x)^(",i, "/60)", ") ", collapse="")
}
#complete
equation=paste(toString(net), eq1, " == ", toString(gross), collapse="")

私はそれを使用してそれを解決しようとしますSolve(equation, "x").

4

1 に答える 1

1

これは、APR の方程式のように見えます。代わりに、次のような単純な反復を試してください。

#inputs
instalments=60
net=12800
monthly=387.63
interest=0.1890

#function
CalculateAPR <- function(InitialPayout, InterestRate, N, MonthlyRepayment) {
  i <- InterestRate 
  repeat{
    DF <- sapply(1:N, function(N) { MonthlyRepayment/((1+i)^(N/12)) } )
    if(InitialPayout>=sum(DF)) break()
    i <- i + 0.00001
  }
return(i)
}

#output
ans=CalculateAPR(net, interest, instalments, monthly)
rm(list = c('instalments', 'interest', 'monthly', 'net'))
print(ans)

各反復に 0.001% を追加するだけの、このアルゴリズムよりも効率的なアルゴリズムを試してみることをお勧めします。

于 2013-09-24T12:32:14.140 に答える