수치해석(경계치 조건)

소개글

수치해석 경계치 조건에 대한 수치해석 입니다.(A+)
프로그램 및 쉬운 주석
화면캡처 등등...

목차

프로그램 소스
화면캡처
오차분석

본문내용

#include <iostream> //입출력 헤더파일
#include <cmath> //수학적 함수표현 헤더파일
#include <iomanip> //출력문구에 공간을 할당 해줄수 있도록 하는 헤더파일
using namespace std; //네임스페이스 할당
#define e 2.7*************52*************266 //e값에 대한 정의
double analytic(double x); //아날라틱 솔루션에 대한 함수
int main()
{
///////////////////변수선언/////////////////////////////
double mat_A[101][101]={0,}; //좌측에 놓일 행렬
double mat_D[101]={0, }; //우변에 놓일 행렬
double s=0; //문제에서의 우변
double q=0; //문제에서의 T항의 계수
double h=0; //등간격
double r=0; //tridiagonal 해법에서 사용될 임시변수
double x=0; //아날리틱 솔루션에 사용할 변수
double t[101]={0, }; //실제 출력값에 대한 변수
double b[101]={0,}; //tridiagonal 해법에서 사용될 변수
double b_prime[101]={0,}; //tridiagonal 해법에서 사용될 임시변수
double d_prime[101]={0,}; //tridiagonal 해법에서 사용될 임시변수
double eps=0; //오차에 대한 변수
int i=0; //for문에서 사용할 변수
int j=0; //for문에서 사용할 변수
int n=0; //등간격 갯수
double x_0=0; //T(0)
double x_n=0; //T(10)
////////////////과학적 표기및 숫자표기법//////////////////
cout.setf(ios_base::left,ios_base::adjustfield);
cout.setf(ios_base::showpos);
cout.setf(ios_base::showpoint);
cout.precision(4);
ios_base::fmtflags old=cout.setf(ios_base::scientific,ios_base::floatfield);
//////////////////입력값 입력////////////////////
cout<<"초기값 첫구간 입력=>";
cin>>x_0;
cout<<"n초기값 끝구간 입력=>";
cin>>x_n;
cout<<"n등간격으로 나눌 구간 입력=>";
cin>>n; //각 입출력
/////////////////행렬 A를 정의////////////////////
h=(x_n-x_0)/n; //등간격에대한 정의
for(i=1;i<=n-2;i++)
{
mat_A[i][i+1]=1/pow(h,2);
mat_A[i+1][i]=1/pow(h,2);
//차분식에서 T(i-1),T(i+1)의 계수를 의미함
}
q=(-0.1); //q값정의
for(i=1;i<=n-1;i++)
{
mat_A[i][i]=2*(-1)/pow(h,2)+q;
//차분식에서 T(i)의 계수를 의미(행렬에서의 대각선 부분)
//나머지 부분은 변수초기화 하면서 0으로 다 초기화 하였음
}

참고 자료

없음