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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
5 a" Z! ]* e# q#include "stdafx.h"
( p, P! ~& E- k7 o
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
! G$ o* k5 {! K
  `7 M4 k. S3 F8 m#include "mathlink.h"

% N! B" p7 x% i) R: tdouble f(MLINK lp,double x);
- W2 F, i0 w1 z; J5 L  C
- M0 u* J; `* O: t8 X* Dint main()
* p+ K4 h% r! S$ ~. x{
        int argc=4;
9 a$ v6 R7 s+ i5 M0 d7 C; T        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
+ p% \6 p5 w6 |& q( P3 T+ a        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
        env=MLInitialize(NULL);//初始化
/ q7 M1 q5 ~# j7 [" ^        if(env==NULL)
' t9 g: o) e( H& U        {
# A- r4 F9 p$ K/ J                printf("initialization error!\n");
8 }  D, g9 X' \- Y( j                return 1;
        }
2 f/ d. q$ l* X" ?- Z3 e9 t        lp=MLOpen(argc, argv);
        if(lp==NULL)
        {
- J4 f4 h* {7 T" G; V) s2 O  K  X                printf("open error!\n");
$ \- N/ j* S, A0 _                return 1;
        }
        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
0 i3 d5 I1 O  h  ^) x- j. r
% ~/ q3 d2 y& m* g        MLClose(lp);
: A$ Y* U$ G& J/ Y8 U        MLDeinitialize(env);

6 P* Q% D' j; D8 s' N; Y        //printf("Hello World!\n");
2 S& i! N  }+ a! z& b( R( k3 h* K7 u        return 0;
0 F8 s( t) }$ y$ u}


double f(MLINK lp, double x)
$ R0 E5 m0 n8 x{
        double res;
0 R5 W% i7 p% Y4 X! [        int pkt;
        static int done=0;
$ z% p8 U( s7 P  T+ ?% R        if(!done)
        {
                //char* expression="cos(x)^5";
+ h$ Z) O2 Y5 P2 l" ?: `; {                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
, T5 F4 e  n& Y" p% j3 s                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
" J" A) J% y: ]                MLNewPacket(lp);
                done=1;
        }
! Q' {1 R3 s9 l9 y- p" X% |8 Z# ^        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
0 |& a- ?! M; U% m+ a& I5 ~& S        MLPutReal(lp,x);
        MLEndPacket(lp);
2 n/ b& B' K9 S3 z7 z) g        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
4 y0 ?  O, F* H6 q: a' ~        MLGetReal(lp, &res);
        return res;
}

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？