24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

poko宝贝

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1624.1
散金: 100
帖子: 455
在线: 29.9小时
虫号: 485268
注册: 2007-12-27
性别: MM
专业: 高分子物理与高分子物理化

[求助] 求助！！！使用 Numerical recipes in C++求解非线性方程组已有1人参与

Golbally Convergent  for Nonlinear Equation
调用的newt函数。没有学过指针函数什么的，照葫芦画瓢写的很烂。编译提示错误是：
newt.cpp(22): warning #159: declaration is incompatible with previous "fmin" (declared at line 342 of "/usr/include/bits/mathcalls.h"
float fmin(float x[]);
newt.cpp(38): internal error: assertion failed at: "shared/edgcpfe/lower_c99.c", line 2870
f=fmin(x);
^
code如下,其中void *vecfunc(int n,float x[],float fvec[])是很复杂的非线性方程组，这里为了简洁易见，简单列了几个方程，可能解不出来。
Numerical recipes in C++书里格式是 void (*vecfunc)(int n,float x[],float fvec[])，但是调用的时候老出错，就把括号去了。
整个调用的函数包已经打包在附件，急用，请大神帮忙。
//////////////////////////////////////////////////

#include <math.h>
#define NRANSI
#include "nrutil.h"
#include "nrutil.cpp"

#include "newt.cpp"
#include "fdjac.cpp"
#include "fmin.cpp"
#include "lnsrch.cpp"
#include "lubksb.cpp"
#include "ludcmp.cpp"

void *vecfunc(int n, float x[], float fvec[]);
int main()
{
int n;
n=5;
float x[n];
float *fvec;
int *check;

check=ivector(1,n);
fvec=vector(1,n);

newt(x, n, check,  vecfunc);

free_ivector(check,1,n);
free_vector(fvec,1,n);
return 1;
}

void *vecfunc(int n,float x[],float fvec[])
{
//float x[5];
//float fvec[5];

fvec[1]=x[1]*x[2]+x[5]-1.0;
fvec[2]=x[1]+x[2]-x[3];
fvec[3]=x[1]*x[4]+x[2]+x[3]-10.0;       
fvec[4]=x[1]+x[2]+x[3]-1.0;
fvec[5]=x[4]+x[5]-1.0;               
}

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : Golbally_Convergent__for_Nonlinear_Equation.zip

2017-02-09 13:21:52, 12.55 K