mathematica c语言混合编程的问题 大神求救

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
$ F0 V$ x  P. x  x2 g+ [) x#include "stdafx.h"

#include <stdio.h>
! V6 n; F8 l! o! R% F/ k#include <stdlib.h>
#include <string.h>

#include "mathlink.h"
, R% Z$ ^2 T- k4 |
double f(MLINK lp,double x);

% j& T3 T4 R( G/ ?  c5 }5 [' Bint main()
' V6 L2 c- ?8 `: C5 Q) S+ x) Y" M{
9 N3 m; O$ Y2 X% A& K) n, n+ y        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
9 R' v6 K+ z7 A0 H        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
        env=MLInitialize(NULL);//初始化
        if(env==NULL)
        {
5 }. B5 e" E! S/ r$ l% _                printf("initialization error!\n");
3 _+ J/ W) x7 T- \                return 1;
        }
        lp=MLOpen(argc, argv);
        if(lp==NULL)
3 P( F7 w. h/ O! t  B        {
: Y4 k: y" k3 l& S' b' w) s" e                printf("open error!\n");
                return 1;
. D/ {9 D" Q: H* I  h9 S  y        }
: {: t; B" y3 t; D4 U  M        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
( K' W; C! V7 c  J: N% r        MLDeinitialize(env);
6 W; ?0 W$ \8 `! }2 r
        //printf("Hello World!\n");
        return 0;
5 F4 A! |6 u7 x9 s$ t) c" O}
* M1 o$ O9 _3 a5 S! E

+ c9 A: }; i! U" ?1 e: Bdouble f(MLINK lp, double x)
0 R/ u* @2 x% E7 h' J8 P{
3 n' B) c: s. q; D7 s9 F        double res;
        int pkt;
$ d3 X4 w6 n, r        static int done=0;
* ?+ e3 U+ n6 w0 h; [; I$ X        if(!done)
' g+ l1 a$ E* r" J; H3 U        {
) x0 o( B) I( Q, {( f                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
  ?0 d. J7 w7 M: ?$ B7 n6 v             MLEndPacket(lp);
: ~3 R8 O2 @+ L$ h4 K  i                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
3 H; h. ^5 J2 S                MLNewPacket(lp);
( V0 x; t4 A+ t& j% A0 G) U                done=1;
        }
+ E6 K5 v% F9 f        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
* p6 Y4 G! P$ t) Z. }        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
8 k4 k+ Z' ]& o! x2 @0 M. u        MLGetReal(lp, &res);
2 _  X  a6 \8 k7 U1 B        return res;
}
, J+ \* P, P, N6 Z

遥望科西嘉

我觉得是不是因为在mathematica上面还要写个接口，提供c数据交换的接口

遥望科西嘉

这么大个论坛没人了么？