PDE.zip_C Builder

This is a 2D heat equation Partial differential equation.





󰇛













Here we are taking two cases

Case 1: when 

Case 2: when +T





























































󰇜  󰇛







󰇧









󰇨



󰇧





















󰇨

 







󰇜

The solution for Case 2 is

The difference is only in the factor of extra T in case 2

Because of this the value of T at different x and y values are increase which is show in the graph 2.

So the C++ code for Case 1 is given as

#include<iostream>

#include<fstream>

using namespace std;

int main()

{

std::cout<<"Enter the time ";

double A = dt/(dx*dx);

//double B = 1/dt;

ofstream outfile;

outfile.open("data1.dat");

//Inital condition

float x = 0;

float y;

for (int i=0 ; i<n ; i++)

{

y=0;

x = x+dx;

}

// Boundary

for (int k=0 ; k<=n_step ; k++ )

{

for(int i=0 ; i<=n ; i++)

{

T[0][i][k] = 1;