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

遥望科西嘉

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

#include <stdio.h>
#include <stdlib.h>
9 @3 m+ d1 A8 @2 y2 a#include <string.h>

#include "mathlink.h"
! o( s. P8 p- v1 E
double f(MLINK lp,double x);

int main()
{
2 h' p- W0 S, ?        int argc=4;
9 ?3 e7 L6 p- b8 A* y# y, E1 I        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
( \9 I. W) u: R2 `' N. k        MLEnvironment env;//mathlink的环境变量
- B# C+ F% [# z( |2 L" ~& h        MLINK lp;
        env=MLInitialize(NULL);//初始化
- }7 q) _8 o- q8 m( B; z$ V        if(env==NULL)
6 d6 G# Y8 B( }2 O2 s4 A' I; o2 T2 i        {
( ?7 ?8 X! c$ _+ [                printf("initialization error!\n");
                return 1;
6 X( r% s7 Q; t0 J' G5 g( ]; v) R        }
        lp=MLOpen(argc, argv);
& o8 }( ^+ M5 I& V( C, ~0 S        if(lp==NULL)
        {
                printf("open error!\n");
                return 1;
: p  o: I1 W6 a* l8 M; M        }
+ M# ?, N4 d4 j. V  d2 Z7 f        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
        MLDeinitialize(env);

3 J# K$ E' e  x" Y* e( a; o        //printf("Hello World!\n");
$ {+ b5 O: G" n/ v! t* D        return 0;
9 d9 }5 o0 `9 U: e& O/ |}

( J- i/ z/ k' {+ W& i
- j% \4 N/ [- U% }: r9 M1 m( tdouble f(MLINK lp, double x)
{
6 a! U! d) {3 ?1 J+ `' f        double res;
        int pkt;
        static int done=0;
        if(!done)
& |! C1 ~5 ^- w, D" j  A) E7 w        {
                //char* expression="cos(x)^5";
, O" K$ d8 U! ~: r! b                MLPutFunction(lp,"EvaluatrPacket",1);
5 z" Y7 ~, `! d                MLPutFunction(lp,"ToExpression",1);
: ^+ ^, V  B, x                MLPutString(lp,"cos(x)^5");
8 _- X6 a! p) u4 P4 b; ^* a+ A; m/ H             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
- Y6 }/ L) {& V                MLNewPacket(lp);
                done=1;
        }
        MLPutFunction(lp,"EvaluatePacket",1);
& }' l5 V( A* n, T        MLPutFunction(lp,"N",1);
. G- O4 E( [# x+ f$ \        MLPutFunction(lp,"F",1);
$ r# h2 Y( _& a2 k2 F: Q        MLPutReal(lp,x);
0 @0 [, e0 E7 q/ f5 R$ j  b; ?        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
/ _& Y* n) g) f        MLGetReal(lp, &res);
+ {6 |2 f. \6 [8 Z& X0 Z9 Q" n        return res;
. d" y" ?/ k% }: n4 i9 y% R}

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？