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

遥望科西嘉

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

#include <stdio.h>
" l8 `5 o- i) P% `7 N#include <stdlib.h>
#include <string.h>

$ a) c" N1 ~, Z#include "mathlink.h"
; z* i* b- O. d( B) R- |
0 |* V1 K. P& rdouble f(MLINK lp,double x);
+ {# V7 U' j/ }9 p. f8 ~
int main()
{
! c8 j" V4 j8 n& A        int argc=4;
0 N. T: {1 a- R# }, E7 @5 x        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
        env=MLInitialize(NULL);//初始化
% P/ `) z6 @/ n) z: l        if(env==NULL)
  @' j" r3 x% H5 {" t4 h        {
                printf("initialization error!\n");
                return 1;
" o0 E+ n( Y4 m9 [3 \        }
        lp=MLOpen(argc, argv);
        if(lp==NULL)
# w1 Y$ k, P* Y* T& X        {
                printf("open error!\n");
3 s0 t. l, V, ~( R                return 1;
7 Z4 r" U8 ^0 t6 b( i2 F        }
' _  W' y* G: X* n        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
& {6 M, y( N# T/ M4 o- L9 A- P- E
. _  e5 R! Z" [( _' e4 g3 f# z        MLClose(lp);
        MLDeinitialize(env);

        //printf("Hello World!\n");
        return 0;
}

6 O: b0 ^/ g. C. j- C% L: C% l  }
double f(MLINK lp, double x)
{
$ C6 N. y0 A' z/ B        double res;
3 P# N  A* g3 s9 E        int pkt;
2 E# m3 R, ~6 K' N2 e# S. |        static int done=0;
        if(!done)
7 t: p1 X6 M9 ^$ e9 Y        {
                //char* expression="cos(x)^5";
$ B; L' f9 `& q6 {                MLPutFunction(lp,"EvaluatrPacket",1);
; ~# w  w5 a* R7 Q& s; {5 Y+ e                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
  [# k8 {9 d$ _             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
                MLNewPacket(lp);
! P: A3 h6 W4 C+ c/ e! n5 u                done=1;
        }
5 \/ E; B2 }& N8 T! G        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
5 M9 N8 A7 y3 p2 [! W( C* H, p        MLPutFunction(lp,"F",1);
. f, p1 I5 m+ v6 v6 j& f/ v+ n& \" g        MLPutReal(lp,x);
        MLEndPacket(lp);
3 j, Z& y. }' R( h) X' p% ~! O        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
        MLGetReal(lp, &res);
6 L0 T4 h7 c! {- a' r. \5 M        return res;
}

2 p$ w: ]  V( H% `. B% {/ E

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？