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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
& M  V1 U6 ~+ H% E#include "stdafx.h"
- E- p3 u  Z1 K! v
$ ]# b0 q9 W8 V#include <stdio.h>
; O1 E& q5 S6 d, e#include <stdlib.h>
#include <string.h>
# P% p3 [& v! V1 X
#include "mathlink.h"

" q' ]) T& @/ k0 m8 y3 X4 Ldouble f(MLINK lp,double x);

8 H& c0 x2 j8 m5 Yint main()
$ g: r% \$ f2 R* E{
        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
1 a% e4 M' @  w+ C' e% O        MLEnvironment env;//mathlink的环境变量
8 H% m! P  @: Z/ @/ u3 M  `        MLINK lp;
* p, B2 z4 c7 x1 [$ ?3 M' l6 W        env=MLInitialize(NULL);//初始化
* Z7 K+ ]( l, P% N" \4 J        if(env==NULL)
        {
                printf("initialization error!\n");
                return 1;
        }
( `( x7 k/ {2 d$ g  b. L        lp=MLOpen(argc, argv);
        if(lp==NULL)
( b  }2 J% U3 _% Z" q        {
$ I; a, ]% ]9 l9 c3 Y                printf("open error!\n");
                return 1;
        }
: }) Q. G; w! H, U5 f1 [        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
0 i$ h0 q1 y" S
        MLClose(lp);
% E7 K$ r9 o; f2 Z! @        MLDeinitialize(env);
! `4 f. g3 w( t. a
- ?/ Y! z$ f) B" I( E# x        //printf("Hello World!\n");
+ g, k" \" x1 \2 k        return 0;
4 z# S; }/ L- J& ?+ Y}
  I, E0 q% {8 L& M5 Y; h. A

double f(MLINK lp, double x)
- z1 V( }. o1 `4 e: \  `{
5 ?5 c4 p; j: o        double res;
        int pkt;
        static int done=0;
) m; s/ V' A+ ]4 g        if(!done)
/ a. r  }3 q  `+ d) a6 x! s        {
* k& l0 F3 e2 S' y                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
                MLPutFunction(lp,"ToExpression",1);
8 `& X$ a2 @( ^9 E/ |/ r  I: y                MLPutString(lp,"cos(x)^5");
' z8 P. m9 u5 P5 A) b& p" ~             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
7 ]" d& g" D( _: S* E6 V  F                MLNewPacket(lp);
                done=1;
; Y1 F0 e  |: h) c6 m( U        }
- q; V4 y- U7 @        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
8 o. n6 h: q5 ]- i7 l4 u        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
        MLEndPacket(lp);
+ Z# e, F+ a6 R* }        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
  b4 M" ]) M! }8 S        MLGetReal(lp, &res);
/ [4 H8 H5 Y0 O" M5 j+ \6 W" {9 O% M9 a        return res;
2 R) k# E, N8 j0 c! p}
/ S5 |. J0 E* S/ \0 F- w
, ]# T6 c- n: \, I! F' E" b- m

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？