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

遥望科西嘉

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

#include <stdio.h>
! O$ w0 q7 k: T6 P. g#include <stdlib.h>
  e1 e; ~7 G" L8 \7 @! ^' b#include <string.h>
( M: {- C8 M) f! \
" R7 g  |( Z: y+ O3 X* t! b#include "mathlink.h"

double f(MLINK lp,double x);
9 p5 X" v" o( Y4 ~
# D+ i! t) D" i2 Kint main()
{
( k& I1 P* G6 G" u$ w0 a$ G! Q        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
3 [4 g0 H; a3 `3 P: j) g5 n0 E9 R8 ~% `        MLEnvironment env;//mathlink的环境变量
( z2 n6 b7 j+ ]% ~( w2 Z8 W+ i        MLINK lp;
        env=MLInitialize(NULL);//初始化
" Z. j( ~9 p/ P* E$ ^        if(env==NULL)
/ Z- L9 X9 H. _) j: B& L9 ?        {
                printf("initialization error!\n");
                return 1;
        }
  j- T) A( S  j: c+ R8 C        lp=MLOpen(argc, argv);
# C0 O( N% ~- Y" o: H        if(lp==NULL)
        {
                printf("open error!\n");
                return 1;
        }
2 C5 _' S6 v# I  @1 |& l        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

! S% |. y( D5 D/ _& H* I( [        MLClose(lp);
% E+ D3 w9 Q. z* Z; d        MLDeinitialize(env);
/ G' v# v+ G- ]/ E+ l
( ?6 n( j7 Z6 q; O8 q        //printf("Hello World!\n");
" [3 @" ^5 U8 S  C) O" Q! }+ P) H        return 0;
0 y& P! W* ?/ j3 q/ {0 X7 H0 y# X}


6 h& u6 T& t1 O; A/ u$ hdouble f(MLINK lp, double x)
{
        double res;
/ X  v1 [& i$ {% e2 _( H9 y, M        int pkt;
- J! L+ [& g- I: f        static int done=0;
, T( a7 t1 A, I4 C7 F6 p! e        if(!done)
        {
                //char* expression="cos(x)^5";
# {3 u" Q9 L7 v; X. K5 b  `- K                MLPutFunction(lp,"EvaluatrPacket",1);
1 ?. ]. o: R$ D( r                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
; }0 o; Y, _$ ^- b+ e             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
8 v9 `! z. }* W' ~/ Y                MLNewPacket(lp);
                done=1;
/ }8 [3 q2 X6 }: [. r        }
        MLPutFunction(lp,"EvaluatePacket",1);
& F& j0 e. o1 y        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
        MLPutReal(lp,x);
$ w, |8 _- A! v: X5 \: d7 ]- Y        MLEndPacket(lp);
' O8 P) _2 u3 s" ]: s# ^        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
        MLNewPacket(lp);
1 d* Q* T4 u, l        MLGetReal(lp, &res);
5 J- b8 {+ x) ^* S5 |/ i, ]        return res;
7 j6 t+ M1 L2 t% J( q; Q; O}

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？