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

遥望科西嘉

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

#include <stdio.h>
% m+ d/ x( V% @, Q, O#include <stdlib.h>
#include <string.h>

% E6 v6 N; f4 N#include "mathlink.h"

double f(MLINK lp,double x);
$ l- n  {. F8 H7 ^$ y! J. M2 C
int main()
{
! t  l. ~, g4 B% h. z        int argc=4;
        char*argv[5] ={"-linkname", "F:\\Program Files\\Wolfram Research\\Mathematica\\8.0", "-linkmode", "launch", "NULL"};//指定连接方式
' j* u- n- G3 _3 |/ f        MLEnvironment env;//mathlink的环境变量
        MLINK lp;
. K# b! G, o. d        env=MLInitialize(NULL);//初始化
# }9 l) |/ \9 i, T8 A3 m        if(env==NULL)
# |% A( n, Y4 C! @) P" _        {
/ b& n, C( t- P% v3 _+ C0 a) _9 \* T                printf("initialization error!\n");
                return 1;
: c9 _  O, q2 S: ]  g8 K) y        }
  }& r2 F3 ]( ~2 t0 J8 v        lp=MLOpen(argc, argv);
$ X- ~; T9 E+ [6 j        if(lp==NULL)
  x& ^' A' W. q* G; p; m0 l        {
                printf("open error!\n");
                return 1;
        }
$ ]! W2 }6 X3 x9 `' h        for(int i=1;i<=10;i++)
                printf("cos(%d)=  %lf\n",i,f(lp,i));//计算f（1）至f（10）

1 ?  y2 q6 {# u! C7 {/ g9 x        MLClose(lp);
        MLDeinitialize(env);

        //printf("Hello World!\n");
. p0 i8 O! V1 Y6 ?0 e- p        return 0;
1 k! X" _2 Z) t5 V# C& ?}
: k  k6 a4 y7 ]6 p3 R# I" R

double f(MLINK lp, double x)
- ^2 ], a$ u+ r7 I# Y9 B9 l{
$ ?6 v4 w# N& D6 B9 w        double res;
. {9 d% L2 J9 `$ I        int pkt;
        static int done=0;
9 r' \7 F2 M; ~  O& H/ e/ @# _, ]9 V        if(!done)
, H2 T" R& G/ X# |/ @        {
                //char* expression="cos(x)^5";
                MLPutFunction(lp,"EvaluatrPacket",1);
7 h# T1 d* I: d$ _; P- h; x                MLPutFunction(lp,"ToExpression",1);
( {/ d) p0 ^. D+ v' c) C                MLPutString(lp,"cos(x)^5");
             MLEndPacket(lp);
                while((pkt=MLNextPacket(lp),pkt)&&pkt!=RETURNPKT);
: f; t1 @: t9 V( Z# P                MLNewPacket(lp);
                done=1;
        }
        MLPutFunction(lp,"EvaluatePacket",1);
        MLPutFunction(lp,"N",1);
        MLPutFunction(lp,"F",1);
5 [- a$ U; Z" S, {7 {; N$ {7 |        MLPutReal(lp,x);
% F. W  ]) |# J4 r        MLEndPacket(lp);
        while((pkt=MLNextPacket(lp))&&pkt!=RETURNPKT);
8 _$ D0 J# q$ r0 y( Z        MLNewPacket(lp);
; N& [4 i# ]. N        MLGetReal(lp, &res);
        return res;
}
2 B0 c7 K1 e0 {( g& m+ a

遥望科西嘉

这么大个论坛没人了么？

遥望科西嘉

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