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

遥望科西嘉

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

#include <stdio.h>
#include <stdlib.h>
! R3 A7 v% P( P& p+ M6 ^* E#include <string.h>

#include "mathlink.h"
; f% y& {9 w) M1 Q5 M- A$ j. w
double f(MLINK lp,double x);

int main()
8 A, R) Z1 N& t" H) u8 `; l{
$ l" M9 \, `" u0 I! l        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
' x$ T) [* _# P# F0 Y5 W. p+ m        MLEnvironment env;//mathlink的环境变量
: K* N9 X6 F$ t        MLINK lp;
* b4 o4 z, k" l: B- ~5 o' O        env=MLInitialize(NULL);//初始化
4 B& {1 E& t( H* l        if(env==NULL)
# X8 \4 o0 B  @1 j  x        {
. N4 ?( r2 o! E# N& C                printf("initialization error!\n");
, c8 j8 F; n; L% ]$ A; _0 ~                return 1;
  f, Z( N6 o8 z3 T+ h        }
        lp=MLOpen(argc, argv);
: _3 L) A3 V. c" c7 {2 i        if(lp==NULL)
: ^$ Y; J7 x% H  Z2 k% @. E        {
                printf("open error!\n");
                return 1;
        }
        for(int i=1;i<=10;i++)
; A5 H4 r- ^0 \" D; l0 _                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

. s9 L) ?7 D# H7 {& M+ u        MLClose(lp);
0 U5 Z6 v& j9 N2 n$ E$ y% {3 ]1 Q        MLDeinitialize(env);

        //printf("Hello World!\n");
- _; C, {6 \" @) j1 [5 [1 r' `        return 0;
}


% R; w* h9 B3 @; S* D3 M: l" tdouble f(MLINK lp, double x)
( W2 }+ }' d  U* F3 J* ~{
        double res;
        int pkt;
6 V0 `( W8 Q4 t! r9 Z        static int done=0;
( ^  {2 m% D/ D* z+ u, N: n        if(!done)
" V) L! b2 I9 \; o        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
' M7 O4 i/ t7 W3 Y' \) W, [                MLPutFunction(lp,"ToExpression",1);
! @! Y+ P( ^0 c. ~* K                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
                done=1;
        }
' i1 B$ z# l6 I$ M* P        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
! l! S6 f; U* X        MLPutFunction(lp,"F",1);
! I; J+ Q4 D5 [2 e" L3 l% V        MLPutReal(lp,x);
- w" V' u6 D! n, b1 E+ S: [        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
% R- ~4 }# w# ]        MLNewPacket(lp);
# h$ ~8 }/ ]0 Z( m/ x9 L        MLGetReal(lp, &res);
5 o4 }( a( x( B$ }' y2 q        return res;
}
" E+ |1 W7 ?+ f9 l8 a. f

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？