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

遥望科西嘉

目的是算出cos(x)5次方的程序，整个程序为，结果却不对，调试了很长时间也没有搞定，向大神求助了
#include "stdafx.h"
" b4 `# L; ~# r$ ?4 c% r( D
+ ]6 G5 ], p  [- e7 R#include <stdio.h>
2 E2 M: b+ W" A#include <stdlib.h>
#include <string.h>
) w: R  W  A) L" s3 p8 ]$ j; Q
#include "mathlink.h"

& i: g$ A- W# x) adouble f(MLINK lp,double x);
# \$ }# M- ]# V' s! k$ a4 @$ k
! O1 ?! k9 ^4 P' B% Oint main()
{
        int argc=4;
1 r! S% G) k6 X# L        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
2 u8 ~' s$ d( B1 h$ i        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
4 x/ }+ u' D# V2 w' r. k# F4 b        env=MLInitialize(NULL);//初始化
6 b0 {# D* i2 L' t        if(env==NULL)
; B& O6 n' {  M" q, A7 s+ X, Q        {
                printf("initialization error!\n");
, ^  a  z% U1 G% T* p9 _0 g' L2 o                return 1;
. Y  G; q5 Z# Y* U& O+ |" _        }
        lp=MLOpen(argc, argv);
& z% t' o. x/ B9 |        if(lp==NULL)
        {
                printf("open error!\n");
" O$ `; |7 i' u, {% b                return 1;
9 v* x# y8 O( c0 H8 u        }
        for(int i=1;i<=10;i++)
2 H+ w  W; m+ x" s' m" f                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

        MLClose(lp);
% U- _2 |/ U) Q$ H/ R  ~/ d0 F9 b        MLDeinitialize(env);
' B4 h, l4 G! v' D( g. [
        //printf("Hello World!\n");
! E' p* n9 {2 `- s5 V        return 0;
  V/ V; S8 |  i0 A}


. I8 x! [. b" ?2 wdouble f(MLINK lp, double x)
{
        double res;
        int pkt;
+ j) J+ Y8 G6 \. m        static int done=0;
        if(!done)
        {
                //char* expression="cos(x)^5";
2 |3 E9 L) h9 |- o1 y                MLPutFunction(lp,"EvaluatrPacket",1);
, v+ @/ h7 P8 ^) V, b: A; q                MLPutFunction(lp,"ToExpression",1);
5 z+ T! W- O9 g6 t6 M* b/ c; u3 N                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
) ~  ]5 e( _8 |5 J0 E9 u" U                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
7 f! n: \6 B% r) i! n7 F                MLNewPacket(lp);
                done=1;
! U3 E/ t# [5 ]9 k5 X        }
0 s" _+ w/ \3 Z3 U  Y/ G+ l" i        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
/ p- [4 Y+ c  M* V, K* ?" g7 v        MLPutFunction(lp,"F",1);
# P# b  b6 j7 v' r) Q. A9 ?! c        MLPutReal(lp,x);
        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
3 t; A2 Y" o# ~4 G1 ?8 T/ N$ n  @; E3 C1 u        MLNewPacket(lp);
( {% }3 Z, l& u4 Q+ W4 _+ H        MLGetReal(lp, &res);
1 {- c" h$ s" u0 z" l, k        return res;
}

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？