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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
#include "stdafx.h"
/ v3 K5 v8 y# |
#include <stdio.h>
* k8 C2 b4 U0 Z* a5 j#include <stdlib.h>
" q7 ^$ p8 J0 u8 L% }#include <string.h>
+ q/ m2 g# P( c
#include "mathlink.h"
+ h  F+ B% X5 m9 s& O/ {
$ Y. i8 B, B% T) [- ]double f(MLINK lp,double x);
# ^- X. ~9 E4 ^0 c
int main()
5 u: T: a# g& k( e/ f4 K: _7 @$ [6 X{
        int argc=4;
; a$ j6 V; ~; L& c& T. L, t        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
$ u& Y, v. |! V* J$ s2 ]        MLEnvironment env;//mathlink的环境变量
* ]2 U0 ?4 a3 h( b9 F5 d, `        MLINK lp;
% L- C9 {1 g" O' K9 p; j& }        env=MLInitialize(NULL);//初始化
        if(env==NULL)
        {
                printf("initialization error!\n");
- R$ {9 n: c7 a" {& F                return 1;
; \, ~* f; m% @        }
0 D: J) \% J8 x        lp=MLOpen(argc, argv);
        if(lp==NULL)
5 e4 o. R; M- W7 E        {
4 B, {7 G" x7 S! x4 F" y) a                printf("open error!\n");
1 M# v  Z# I7 i5 T  F. S$ l1 G                return 1;
* x3 f2 a4 ]$ w0 f        }
) W& }& E$ P% i( P2 Y        for(int i=1;i<=10;i++)
) E% ^1 }& Z, n7 s                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
) A. f: f2 H. n' Z; K  f# U, ]0 q" A
        MLClose(lp);
        MLDeinitialize(env);

        //printf("Hello World!\n");
        return 0;
}
# P5 C5 `5 Q- Z$ y+ T, o1 ?

double f(MLINK lp, double x)
0 C# a- ]3 D: N{
        double res;
        int pkt;
        static int done=0;
& I3 ?( K/ R4 h" \2 f; [/ K* x        if(!done)
" z4 D4 V  l0 u* ~" x" ]4 y1 c- S        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
8 H5 \, }& ]6 L' o2 t$ c3 R6 p: s                MLPutFunction(lp,"ToExpression",1);
& _) E% }. a; s& Q                MLPutString(lp,"cos(x)^5");
; \) i4 a) O; n& x1 M0 ~+ ]% T. g             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
                done=1;
3 U; W: q' F, X  m( R        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
! p; m+ W2 D% R/ q6 x9 V        MLPutFunction(lp,"F",1);
, O6 z2 K) j9 q$ p! T# p        MLPutReal(lp,x);
8 e1 B0 \. R. F$ m( v% K: ?        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
        MLGetReal(lp, &res);
' Q0 |: F8 C7 N" ?        return res;
+ H! U. y4 I( ?, f6 O: Z6 V}
4 {$ \1 r6 `4 S2 N$ T& y7 M

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？