沧浪之水

//remzf.h//实现remez算法#include#includevoid remz(a,b,p,n,eps,f)int n;double a,b,eps,p[],(*f)(double);{ int i,j,k,m; double x[21],g[21],d,t,u,s,xx,x0,h,yy; if(n>20)n=20; //逼近多项式的最高次数为

#include<stdio.h>//四阶龙格库塔void runge(float(*f)(float x,float y),float a,float b,float y0,int N){ float x=a,y=y0,K1,K2,K3,K4; float h=(b-a)/N; int i; printf("第一问\n"); pr...

#include<stdio.h>float f1(float x,float y) //第一问{ return y*y;//方程}float f2(float x,float y)//第二问{ return 0.1*(x*x*x+y*y);}//四阶龙格库塔提供出发值void runge(float(*f)(float x,float ...

#include<stdio.h>#include<math.h>float f1(float x){ return exp(-x*x);}float f2(float x){ return sin(2*cos(x))*sin(x)*sin(x);}float f3(float x){ return sin(x)/sq...

（1）Hu.mfunction a=hu(x,y,m)S = zeros(1,2*m+1);T = zeros(m+1,1);for k = 1:2*m+1 S(k) = sum(x.^(k-1));endfor k = 1:m+1 T(k) = sum(x.^(k-1).*y);endA = zeros(m+1,m+1);a= zeros(m+1,1...

A=[2 4 3 -5 10 -12 4 11];b=[10 -13 -26 25];x=zeros(2,1);m=A'*bn=A'*Ax=n\m结果：m = -179 692n = 129 -83 -83 306x = 0.08172.2836

题目： 算法： #include<stdio.h>#define N 3void difference(float y[],float f[4][4],int n){ int k,i; f[0][0]=y[0]; f[1][0]=y[1]; f[2][0]=y[2]; f[3][0]=y[3]; for(k=1;k&...

题目 拉格朗日插值法流程图 (1)二次插值：#include<stdio.h>float lagelangri(float x[],float y[],float xx,int n){ int i,j; float *a,yy=0; a=new float[n]; for(i=0;i<=n-1;i++) { ...

实验目的 (1) 通过上机编程，复习巩固以前所学程序设计语言及上机操作指令。 (2) 通过上机计算，了解舍入误差所引起的数值不稳定性。实验内容 用两种不同的顺序计算 ，分析其误差的变化。实验步骤 用两种不同的顺序计算，即从n=1加到n=10000和从n=10000加到n=1。程序设计#include<stdio.h>#include<math.h&gt...

1 . 实验目的 (1) 通过对二分法与牛顿迭代法作编程练习与上机运算，进一步体会二分法与牛顿迭代法的不同特点。 (2) 编写割线迭代法的程序，求非线性方程的解，并与牛顿迭代法作比较。 2 . 实验内容 (1) 用牛顿法求下列方程的根： (2) 编写一个割线法的程序，求解上述各方程。 3 . 实验步骤 牛顿迭代法流程图： 割线法：将上图所示的牛顿法流程图的迭...

用列主元消去法解下列方程组： #include<stdio.h>#include<math.h>void main(){void zhu(float *,int,float[]); //声明函数zhuint i; float x[4]; //x为解float c[4][5]={1,1,0,3,4, ...

用LU分解法解方程组 ，其中 #include<stdio.h>void main(){ float x[4]; int i; float a[4][5]={48,-24,0,-12,4, -24,24,12,12,4, 0,6,20,2,-2, ...

编写用改进平方根法解方程组 的程序，并解下列方程组： #include<stdio.h>main(){ float x[6],y[8]; int i; float a[6][7]={0.5,-0.5,0,0,0,0,-1, -0.5,1.5,-0.5,-0.25,0.25,0,0, ...

编写用追赶法解三对角线方程组的程序，并解下列方程组：#include<stdio.h>#include<math.h>void main(){ int i; int n=4; float x[4]; float c[4][5]={2,-1,0,0,5, -1,2,-1,0,-12, ...

分别用雅可比 迭代法与高斯塞德尔迭代法解下列方程组： 雅可比迭代法：#include<stdio.h>#include<math.h>#define eps 1e-6#define max 100//雅可比迭代法void jacobi(float *a,int n,float x[]){ int i,j,k=0; float eps...

#include<stdio.h>#define N 10void modeuler(float(*f)(float,float),float x0,float y0,float xn,int n){ int i; float yp,yc,x=x0,y=y0,h=(xn-x0)/n; printf("第一问\n"); printf("x[...