下载 >  开发技术 >  C++ > 偏微分中的古典隐格式的实现

偏微分中的古典隐格式的实现

2010-05-21 上传大小:467B
偏微分中的古典隐格式的VC++的实现:
#include<iostream.h>
#define N 20
void main()
{
 double u[22][20],x[20],alpha[20],beta[20],gama[20],w[20],g[20];
 double h=0.05,r=0.4,k=16.0;
 int i,n;
 for(i=1;i<N/2;i++)
 {
  x[i]=i*h;
  u[i][0]=x[i];
 }
 for(i=N/2;i<N;i++)
 {
  x[i]=i*h;
  u[i][0]=1-x[i];
 }
 for(i=2;i<N;i++)
  alpha[i]=r;
 for(i=1;i<N;i++)
	 beta[i]=1+2*r;
 for(i=1;i<N-1;i++)
	 gama[i]=r;
 for(n=0;n<6;n++)
 {
  w[1]=gama[1]/beta[1];
  g[1]=u[1][n]/beta[1];
  for(i=2;i<N-1;i++)
  {
   w[i]=gama[i]/(beta[i]-alpha[
                            i]*w[i-1]);
   g[i]=(alpha[i]*g[i-1]+u[i][n])/(beta[i]-alpha[i]*w[i-1]);
  }
  g[N-1]=(alpha[N-1]*g[N-2]+u[N-1][n])/(beta[N-1]-alpha[N-1]*w[N-2]);
  u[N-1][n+1]=g[N-1]; 
  for(i=N-2;i>0;i--)
	  u[i][n+1]=w[i]*u[i+1][n+1]+g[i];     
  for(i=1;i<N-1;i++)
	  cout<<u[i][n+1]<<endl;
 }
}                        
...展开收缩
综合评分:0
开通VIP 立即下载

评论共有0条

 
C语言程序,差分算法,显示格式隐式格式 立即下载
积分/C币:3
偏微分方程数值解法的MATLAB源码--古典显式格式求解抛物型偏微分方程等 立即下载
积分/C币:5
二维抛物线方程交替方向隐格式 matlab程序 立即下载
积分/C币:5
DuFort-Frankel格式求解椭圆-抛物型偏微分方程组 立即下载
积分/C币:3
mathematica 差分法解偏微分方程 立即下载
积分/C币:3
偏微分中的迎风格式的实现 立即下载
积分/C币:3
COMSOL求解偏微分方程模型 立即下载
积分/C币:5
古典加密算法之置换密码和代换密码 立即下载
积分/C币:3
置换密码 古典密码算法 立即下载
积分/C币:5
全部古典密码加密与解密的C语言程序 立即下载
积分/C币:5

热点文章

VIP会员动态

0 1 2
关闭
img

spring mvc+mybatis+mysql+maven+bootstrap 整合实现增删查改简单实例.zip

资源所需积分/C币 当前拥有积分 当前拥有C币
5 0 0
点击完成任务获取下载码
输入下载码
为了良好体验,不建议使用迅雷下载
img

偏微分中的古典隐格式的实现

会员到期时间: 剩余下载个数: 剩余C币: 剩余积分:0
为了良好体验,不建议使用迅雷下载
VIP下载
您今日下载次数已达上限(为了良好下载体验及使用,每位用户24小时之内最多可下载20个资源)

积分不足!

资源所需积分/C币 当前拥有积分
您可以选择
开通VIP
4000万
程序员的必选
600万
绿色安全资源
现在开通
立省522元
或者
购买C币兑换积分 C币抽奖
img

资源所需积分/C币 当前拥有积分 当前拥有C币
5 4 45
为了良好体验,不建议使用迅雷下载
确认下载
img

资源所需积分/C币 当前拥有积分 当前拥有C币
5 0 0
为了良好体验,不建议使用迅雷下载
开通VIP
img

资源所需积分/C币 当前拥有积分 当前拥有C币
5 4 45
您的积分不足,将扣除 10 C币
为了良好体验,不建议使用迅雷下载
确认下载
下载
无法举报自己的资源

兑换成功

你当前的下载分为234开始下载资源
你还不是VIP会员
开通VIP会员权限,免积分下载
立即开通

你下载资源过于频繁,请输入验证码

您因违反CSDN下载频道规则而被锁定帐户,如有疑问,请联络:webmaster@csdn.net!

举报

若举报审核通过,可返还被扣除的积分

  • 举报人:
  • 被举报人:
  • 举报的资源分:
  • *类型:
  • *详细原因: