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

遥望科西嘉

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

#include <stdio.h>
#include <stdlib.h>
' T+ @4 u' l5 z0 i$ @5 u#include <string.h>

4 i! R! B6 o5 o2 W( P' F#include "mathlink.h"
; q/ |: f( e# P, k( C5 E3 B7 e2 }
2 G: P/ z8 N8 q% vdouble f(MLINK lp,double x);

' p" ?% l% J/ Wint main()
0 O9 f* k, e. w2 \. R( ]% k+ [{
        int argc=4;
8 t2 ~4 }/ _& B/ S3 u% I        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
        MLEnvironment env;//mathlink的环境变量
/ R0 [# O2 r! X7 H        MLINK lp;
        env=MLInitialize(NULL);//初始化
# j1 i2 K. G5 S2 L* {        if(env==NULL)
        {
                printf("initialization error!\n");
( X4 t$ C, Z# ~& C- K  W/ O/ L                return 1;
. [/ A/ e' O: e) U        }
( a0 r* x: H1 U2 D, b  [& l        lp=MLOpen(argc, argv);
, C  Y' J" I& C  S  a1 j        if(lp==NULL)
        {
5 i4 [2 W; i( n9 m0 s                printf("open error!\n");
                return 1;
        }
        for(int i=1;i<=10;i++)
! m, r" K2 B6 L! P3 z  \                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
/ v/ M. ]! N9 R
& D; `; j; B2 }5 U7 _8 B        MLClose(lp);
        MLDeinitialize(env);
1 |& B1 W4 U) I3 v2 P6 p
$ _" r# m5 e6 x: n# p! N7 F        //printf("Hello World!\n");
        return 0;
}
( ^$ c4 W5 e/ N, x; _2 r

, O( q* p) F1 h2 rdouble f(MLINK lp, double x)
; a0 Z# s! o( `{
) {% s- J# u9 N" ~  W; A  \        double res;
! o! `) P1 Z+ Q' r7 o7 y# m) T        int pkt;
# b9 U) N# ?. L8 f4 v        static int done=0;
5 a4 \# ~5 U0 h4 B2 s        if(!done)
0 v" S$ I+ H8 @1 U- ^6 A% w, W        {
3 A: Q4 |7 r& M; D) D# n( o                //char* expression="cos(x)^5";
' m$ P! `$ [  M                MLPutFunction(lp,"EvaluatrPacket",1);
0 Y3 r1 b. \  J+ }1 q                MLPutFunction(lp,"ToExpression",1);
+ j% M' G% V+ M+ a' O2 S                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
; @# J5 A( F; p6 j  e, u                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
                done=1;
, U2 B; R. f3 u1 c: y7 Q$ v        }
9 [7 W- O) k, u  z3 `- s  C        MLPutFunction(lp,"EvaluatePacket",1);
9 g( O. x9 h/ \9 G& p4 n& f        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
; @$ N! |: |( f7 z. q        MLEndPacket(lp);
/ E% J: U: C; E2 ^6 n* D% U5 g& g  ?" k        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
        MLGetReal(lp, &res);
        return res;
}
& p5 j, e' o% A- i5 ?  u

遥望科西嘉

这么大个论坛没人了么？

遥望科西嘉

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