二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf

二维有限体积法计算热传导及源码

//#include "stdafx.h"

#include

#include

#include

#include

#include

using namespace std;

#define q 500

#define k 1000

void main ()

{ //input the value you want:

double L,dx,dy,T,Ax,Ay,d;

int m,n,i,j,kk,mm ;

//char str1[20];

char *str1="conf.dat";

//cout<>str1;

//string str1("data.");

ifstream fin(str1);

if(!fin)

{

cout<

exit(1);

}

cout<

fin>>L>>m>>n>>T;

cout<

cout<

cout<

//creat the maxtrix

mm=m*n;

dx=0.3/m;

dy=0.4/n;

vector ae(mm),ap(mm),aw(mm),an(mm),as(mm),su (mm),sp(mm) ;

Ax=L*dy;

Ay=L*dx;

//结¨¢点Ì?0,左Á ¨®下?角?顶£