从枫到C++的翻译

嘿，所以我有一个枫分支方法，我必须将其转换为C++。我试图根据枫论坛上的代码生成帮助转换它，但它一直抛出错误。我将不胜感激这方面的一些帮助。谢谢，从枫到C++的翻译

这里是枫

代码

使用二分法解决以下数学问题：一个。方程的最小正根

f(x):=evalf(1/x-evalf(Pi)*cos(evalf(Pi)*x));

与delta = 10^-5和EPS = 10^-6

plot(f(x),x=.05..10.0);

从图表上方我们可以得出结论，给定的方程已经最小的正实数根位于0.0之间

Bisect:=proc(funct_equation,ai,bi,Mi,epsfi,deltaxi) local k,M,a,b,u,v,w,c,e,epsf,deltax,feq, notsolved: M:=Mi: feq:=funct_equation: a:=ai: b:=bi: epsf:=epsfi: deltax:=deltaxi: notsolved:=true: u:=evalf(subs(x=a,feq)): v:=evalf(subs(x=b,feq)): printf("a=%+9.6f %+12.6e\nb=%+9.6f %+12.6e\n\n",a,u,b,v); e:=b-a; if (sign(u)<>sign(v)) then printf(" n  x   f\n"); for k from 1 by 1 while (k<M and notsolved) do: 
    e:=0.5*e; 
    c:=a+e; 
    w:=evalf(subs(x=c,feq)): 
    printf("%2d %+9.6f %+12.6e\n",k,c,w); 
    if (abs(e)<deltax or abs(w)<epsf) then 
     notsolved:=false: 
    else 
     if (sign(w) <> sign(u)) then 
     b:=c: v:=w: 
     else 
     a:=c: u:=w: 
     fi: 
    fi: od: printf("Root = %+9.6f function = %+12.6e\n",0.5*(a+b),evalf(subs(x=0.5*(a+b),feq))); fi: end: with(plots):

：和2.0

为了与需要我们调用平分法用根隔离间隔（0.01,2.0）的准确性得到它们的值

警告，更名COORDS已经被重新定义

Bisect(f(x),0.01,2.0,30,1.0e-6,1.0e-5):

来源

2011-04-23 blackzodiac

你可能想发布'转变成C++代码和错误。人们可能会发现帮助更容易。 – forsvarir 2011-04-23 19:43:41

evalf（subs（x = c，feq））''feq'是一个函数对象，这看起来特别不平凡。 – Potatoswatter 2011-04-23 20:03:36

没问题，所以根据代码生成论坛上枫叶 http://www.maplesoft.com/support/help/Maple/view.aspx?path=CodeGeneration 我得到的唯一错误是过早输入和呼叫，我使它转换为例如C（功能） – blackzodiac 2011-04-23 20:12:59

你不会需要一个subs电话，如果你把你的feq的过程。

restart: 
Bisect:=proc(func::procedure,ai,bi,Mi,epsfi,deltaxi) 
local k::integer, 
    M::integer, 
    a,b,u,v, 
    w::float, 
    c,e, 
    epsf::float, 
    deltax, 
    notsolved; 
    M:=Mi: 
    a:=ai: b:=bi: epsf:=epsfi: 
    deltax:=deltaxi: notsolved:=true: 
    u:=func(a); 
    v:=func(b); 
    printf("a=%+9.6f %+12.6e\nb=%+9.6f %+12.6e\n\n",a,u,b,v); 
    e:=b-a; 
    if (sign(u)<>sign(v)) then 
    printf(" n  x   f\n"); 
    for k from 1 by 1 while (k<M and notsolved) do 
     e:=0.5*e; 
     c:=a+e; 
     w:=func(c); 
     printf("%2d %+9.6f %+12.6e\n",k,c,w); 
     if (abs(e)<deltax or abs(w)<epsf) then 
     notsolved:=false: 
     else 
     if (sign(w) <> sign(u)) then 
     b:=c: v:=w: 
     else 
     a:=c: u:=w: 
     fi: 
    fi: 
    od: 
    printf("Root = %+9.6f function = %+12.6e\n",0.5*(a+b),func(0.5*(a+b),feq)); 
fi: 
0.5*(a+b); 
end: 

with(plots): 

f:=subs(Pi=evalf[16](Pi),proc(x::float) 1/x-Pi*cos(Pi*x); end proc); 

Bisect(f,0.01,2.0,30,1.0e-6,1.0e-5); 

f(%); 

CodeGeneration[C](f); 

CodeGeneration[C](Bisect);

此外，如果开始时为f的表达，可以随时使用unapply命令把它变成操作员（一种程序，但也可以是代码生成）。

例如，我也可以通过以下方式创建程序f。（注意，这些中的一个产生一个缺省10位近似裨在生成的C代码以及其他的16位数字近似。）

f_expression := 1/x-Pi*cos(Pi*x); 

f:=unapply(f_expression, [x::float]); 

CodeGeneration[C](f); 

f:=subs(Pi=evalf[16](Pi),unapply(f_expression, [x::float])); 

CodeGeneration[C](f);

来源

2011-04-26 20:48:30 acer

从枫到C++的翻译

回答

相关问题