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

遥望科西嘉

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

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
* H1 P( V0 m6 [2 T: [
#include "mathlink.h"

1 e# Q/ Z1 H5 M# ]double f(MLINK lp,double x);

int main()
{
" M8 c0 p$ [! l% \4 {1 @: u        int argc=4;
, r: R; T* q( \, _: R! @: v' N- |        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
! g4 o/ L# J! N) J4 N: a        MLEnvironment env;//mathlink的环境变量
8 K( Q, a5 c' P/ h0 k) R9 \+ x        MLINK lp;
2 z9 R5 y, n5 q/ i. @6 G& k        env=MLInitialize(NULL);//初始化
        if(env==NULL)
* i+ {9 N/ ~  J, |        {
                printf("initialization error!\n");
                return 1;
8 {( t" S! F! V; s2 ]/ U8 a        }
/ w9 h% L& t+ |8 j0 z" H        lp=MLOpen(argc, argv);
        if(lp==NULL)
        {
                printf("open error!\n");
2 [1 ^* _9 S; L                return 1;
        }
        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）
7 A! B6 S! V# [" T7 z* h: u
& c. `4 R6 M, |; i) n        MLClose(lp);
9 A7 Y$ ?- j' W. _( y        MLDeinitialize(env);

        //printf("Hello World!\n");
        return 0;
6 x9 ^  S: ^, X( |; v' `}

4 c4 |6 G$ g: W: q( Z& m8 `
5 \) k; a5 |9 `7 [- N! ^' }double f(MLINK lp, double x)
* v3 w) _2 M+ t- z2 h{
8 w4 J& F) ]& r        double res;
$ c8 B$ J: V8 C# t. w        int pkt;
- a1 v. ?" g1 G7 t" _4 J        static int done=0;
        if(!done)
/ a  e& f# T. P+ ~1 X9 j; ]* B        {
1 M4 |. q* l' D' ~# C                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
  a, G0 l* T9 ]% p                MLPutFunction(lp,"ToExpression",1);
                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
* n5 S$ }% r+ o( f9 ~, s                MLNewPacket(lp);
                done=1;
% ?: n% ?6 v: p) }8 r6 U        }
$ t8 k, M3 c7 s  X        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
& I: C* k9 h! k8 c, T& y; R8 V2 J        MLPutFunction(lp,"F",1);
0 Y. P2 g. ^1 z* U' J        MLPutReal(lp,x);
        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
7 p) m8 @# ^  `) k, z4 X4 _        MLNewPacket(lp);
0 L" J) I. S/ ]        MLGetReal(lp, &res);
9 |8 ?# V, ^% T3 f, Z        return res;
/ N8 B) b9 X# I3 ]! d}

8 y# c: Q- Z7 l! U

遥望科西嘉

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

遥望科西嘉

这么大个论坛没人了么？