Question

The second problem, the one I'm asking about, is only making changes to the first problem, which I have included the correct code for. Thank you!!

een be modeled by diding the plate io adiscrese grid of cels and sulting the change in tenspentureof ell over 90 2.2 3) dared es dlebel css (saggestion:use5 both ves w fint, ut increme this 1o 20asee you get yeur pragranı warking. Your progrars shcald do the Sllasing lain the cs alus of the prid erptures and enay tesp wll be used to ongnethe te vel Poript te Eser·nd efter 'erperature valaes Rr the top, boom, let, teht sides. Also prompt and enter an initiel tepentare T fer the inneriar cels. Fo this leh, you may iniialie ell the inaioe Initielie tenp usingthese values and disply the istial conts of the terp amay on the cesole 2. Copy temp ald 3 Leepova all the ieor cels (bar sot the odge cellst s,aceady nca true afer all cella are checood This is alenper popam thas useal se lere re a number ofhiets reik ahaat theae qu stigul ฝีd dnu ob wish your purt.vr bvtare arv.hy wrinng งเ'y codo. Lay Idouble print α.ต tīt array ary tite, and so te fitxtke, will be useful in debugging your prog om. stake แถเ you can csplan ta fleticg bstween old array th tmparay, and why your pro- stake แถเ yo1tdemnd how to have your propan kxp through only wvrlorof amay Make sare you end what the cenverpesee eterin dest and what role it plays in yoproge #include-iostream> #include-ioma nip> #include using namespace std; const int NROWS-20, NCOLS = 20; void display (double temp NCOLS]); void update(double temp [) [NCOLS], double old[ INCOLS], double tol, bool& steady); int main) ofstream fout; string filenametemp.dat"; fout.open (filename.c str)); double top, botton, left, right, initial, old[NROWS][NCOLS), temp [NROWS] [NCOLS], tol; int i, j; bool steadyfalse; cout top bottom >> left >> right; cout initial; cout tol; for (i=0; ǐ

Answer 1

#include <iostream>

#include <cmath>

#include <string>

#include <iomanip>

#include <fstream>

using namespace std;

const int NROWS= 10, NCOLS= 10;

void display(double temp[][NCOLS], int rows, int cols);

double computeDiff(double oldTemp[][NCOLS],double newTemp[][NCOLS],int rows,int cols);

void update(double temp[][NCOLS], double T[][NCOLS], double tol, int rows, int cols,bool &steady);

int main() {

               ofstream fout;

               string filename = "temp.dat";

               fout.open(filename.c_str());

               double temp[NROWS][NCOLS], old[NROWS][NCOLS], tol;

               int i, j;

               bool steady = false;

               string infile;

               int rows,cols;

               cout<<"Enter the filename : ";

               getline(cin,infile);

               ifstream fin(infile.substr(0,infile.length()-1).c_str());

               if(fin.is_open()) // if file can be opened

               {

                              fin>>rows>>cols; // read the rows and cols

                              i=0,j=0;

                              // read till end of file, reading the values

                              while(!fin.eof())

                              {

                                             fin>>temp[i][j];

                                             if(j == cols-1)

                                             {

                                                            i++;

                                                            j=0;

                                             }else

                                                            j++;

                              }

                              fin.close(); //close the file

                              for(i=0; i<rows; i++){

                                             for(j=0; j<cols; j++)

                                                            old[i][j]=temp[i][j];

                              }

                              display(old,rows,cols);

                              cout<<"Input the tolerance : ";

                              cin>>tol;

                              if(tol <= 0)

                              {

                                             cout<<"Tolerance must be positive"<<endl;

                                             return 1;

                              }

                              while(steady == false)

                              {

                                             update(temp,old,tol,rows,cols,steady);

                                             for(i=0; i<rows; i++){

                                                            for(j=0; j<cols; j++){

                                                                           old[i][j]=temp[i][j];

                                                            }

                                             }

                              }

                              cout<<endl;

                              display(old,rows,cols);

                              for(i=0;i<rows;i++)

                              {

                                             for(j=0;j<cols;j++){

                                                            fout<<left<<setw(9)<<fixed<<setprecision(3)<<temp[i][j]<<" ";

                                             }

                                             fout<<endl;

                              }

                              fout.close();

               }else

                              cout<<"Cannot open file"<<endl;

               return 0;

}

void display(double temp[][NCOLS], int rows, int cols)

{

               int i, j;

               for(i=0; i<rows; i++){

                              for(j=0; j<cols; j++)

                                             cout<<left<<setw(9)<<fixed<<setprecision(3)<<temp[i][j];

                              cout<< endl;

               }

}

double computeDiff(double oldTemp[][NCOLS],double newTemp[][NCOLS],int rows,int cols)

{

               int i,j;

               double sum=0;

               for(i=0;i<rows;i++)

               {

                              for(j=0;j<cols;j++)

                                             sum += pow(newTemp[i][j]-oldTemp[i][j],2);

               }

               return sqrt(sum);

}

void update(double temp[][NCOLS], double T[][NCOLS], double tol, int rows, int cols,bool &steady)

{

               int i,j;

               steady = true;

               for(i=1; i<rows-1;i++){

                              for(j=1;j<cols-1;j++){

                                             temp[i][j]=(T[i-1][j]+T[i+1][j]+T[i][j-1]+T[i][j+1])/4.0;

                                             if(computeDiff(T,temp,rows,cols) >= tol ){

                                                            steady = false;

                                             }

                              }

               }

}

//end of program

Output:

Input file:

3 4 1.0 1.0 1.0 2.0 1.1 1.2 1.2 2.0 1.2 1.2 1.2 2.5 2

Console:

Enter the filename hw5btest1.dat 1.000 1.000 1.000 2.000 1.100 1.200 1.200 2.000 1.200 1.200 1.200 2.500 Input the tolerance 0.1 1.000 1.000 1.000 2.000 1.100 1.163 1.331 2.000 1.200 1.200 1.200 2.500

Output file:

1.000 1.000 1.000 2.000 1.100 1.163 1.331 2.000 1.200 1.200 1.200 2.500