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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
1 H( w; H5 K. }0 L- W4 b#include "stdafx.h"
0 v( }1 \7 a% A
#include <stdio.h>
#include <stdlib.h>
) e$ d/ o, d" v8 {  i% o: t#include <string.h>
. x0 |8 U: H5 j
8 ^. z& \9 s2 j# P#include "mathlink.h"

' ~0 r" b7 ?8 O( ^3 Z% w3 bdouble f(MLINK lp,double x);
9 _, j+ w6 O0 q% }! D" y/ Z
) `$ z6 X% c6 f) {( w; G: B! W) ~int main()
+ I: |7 H8 X. l( u{
        int argc=4;
( u0 \0 c% J; x3 F5 b5 q9 t( c        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
' a  o4 _$ d6 m! {3 e: A# J3 `8 G        MLEnvironment env;//mathlink的环境变量
2 o1 j! H8 v- s7 ?! T        MLINK lp;
& c- i3 _6 }/ R& E6 B3 S        env=MLInitialize(NULL);//初始化
        if(env==NULL)
        {
; {3 \' A; c' p) w0 n                printf("initialization error!\n");
6 ]  s5 c: I* d# S0 w. t) K7 E                return 1;
        }
1 E9 l0 i) g# F6 w5 D( E        lp=MLOpen(argc, argv);
& q  t) h, M  h        if(lp==NULL)
        {
) f7 H/ z- k; X( P' ~' O                printf("open error!\n");
+ Y$ e& Z; _6 _5 _                return 1;
        }
        for(int i=1;i<=10;i++)
4 S6 i9 s2 U7 c9 J6 E                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
( ^3 O: }: U8 x8 W1 U
( Q8 M7 H1 j# d0 ?  }        MLClose(lp);
2 u, b$ k- j8 ?8 `4 e$ V        MLDeinitialize(env);

) Z4 v3 `4 ^. Y4 h+ C! q3 W+ W0 M        //printf("Hello World!\n");
        return 0;
}
# M" h) i8 ^) U9 T4 b

) X! a. D+ {& b8 T& V, ydouble f(MLINK lp, double x)
{
# i, \% W$ L7 M, N        double res;
6 u3 t/ f# g6 a0 T        int pkt;
7 x' H. j  ^5 {$ R8 u' Q3 _  W        static int done=0;
        if(!done)
        {
$ b) }" \8 v6 M% N" o7 ], J  S7 c                //char* expression="cos(x)^5";
0 S# p  y8 G, j" t0 [, `                MLPutFunction(lp,"EvaluatrPacket",1);
7 M; t, ]' r9 Y! w* @" T                MLPutFunction(lp,"ToExpression",1);
& F5 j; M8 H: U% c2 r5 a7 k* `3 e4 k                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
                done=1;
8 y4 o/ i. H' M# o+ K/ n: ]4 x# p        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
4 ^  e+ W2 q5 x7 X6 ^        MLPutReal(lp,x);
        MLEndPacket(lp);
( t5 F6 n6 z8 }$ l, c& l        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
. a$ P9 W/ }6 a$ [        MLNewPacket(lp);
        MLGetReal(lp, &res);
        return res;
}

6 g1 p" B9 v+ B& w' x

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？