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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
#include "stdafx.h"
% t9 c: \8 B8 a; J' N7 A' Q
#include <stdio.h>
7 m: o# c& G8 h#include <stdlib.h>
/ X9 }+ j( a5 n4 S# V#include <string.h>
  l& o. ^1 r8 u6 Z6 o
#include "mathlink.h"

: {- z  r8 O/ C+ A8 g/ `double f(MLINK lp,double x);

int main()
{
" Z! m3 d/ w' d9 }% W8 K        int argc=4;
; q' ?; T: w' P& \5 A        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
        env=MLInitialize(NULL);//初始化
        if(env==NULL)
, I4 S! D1 l, D7 m        {
4 K: D& j, B6 R* Y2 m                printf("initialization error!\n");
                return 1;
        }
; T8 T. t+ y* ]) V        lp=MLOpen(argc, argv);
8 @5 {, F2 e3 A  x4 B        if(lp==NULL)
        {
                printf("open error!\n");
                return 1;
        }
        for(int i=1;i<=10;i++)
: W/ i; U! B: U                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
5 A6 k. J& k6 [7 C2 u) H        MLDeinitialize(env);
; X2 c7 n) ^- H, I
        //printf("Hello World!\n");
        return 0;
& [% Z6 p# n0 f: y0 q( K9 a1 j}
' z* t' H$ H1 D! N1 G0 d: b

1 V5 O# }( l' f0 K9 r! Odouble f(MLINK lp, double x)
{
5 L; e% N* c! U        double res;
        int pkt;
' Q, b+ w4 ]( ]' ?        static int done=0;
        if(!done)
        {
2 x4 Q/ J4 @  O, B                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
: B! z& Y% @- v) P                MLPutFunction(lp,"ToExpression",1);
/ f+ ]& {/ Q' E9 _                MLPutString(lp,"cos(x)^5");
8 g" M3 C6 C  T/ u8 T, B+ J5 X2 ^             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
$ I  O0 b6 u1 o                MLNewPacket(lp);
                done=1;
1 ]0 G+ L) U9 T2 ~+ Q0 C* N        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
6 R8 j1 j* ^9 k" o        MLPutReal(lp,x);
5 ^; d2 i5 e5 q) G8 m' l$ }        MLEndPacket(lp);
5 v7 {; @* g( h* s( V9 W        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
3 f- ~+ A1 D; s+ M; N        MLNewPacket(lp);
        MLGetReal(lp, &res);
        return res;
% o5 J; h$ Q+ r! ?6 d  p}

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？