Question: I am given this matlab code which solves the boundary derivatives for two point forward and two point backward accuracy schemes. I need to change

I am given this matlab code which solves the boundary derivatives for two point forward and two point backward accuracy schemes. I need to change the code to solve for 3 point forward and 3 point backward boundary schemes.

These are the equtions I am solving for:

change the code to solve for 3 point forward and 3 point

MATLAB R2018b-academic use PLOTS EDITOR PUBLISHMEW Find Files Compare Sa Print GoTo Comment %23 Breakpoints RunRun and Advance Run and Advance New Open Save Find Time UsrsKevin SterlingDesktopMatlabLabsKevin Sterling Lab 4Chapter8 Partial Current Folder Editor-CAUsers Kevin Sterling\Desktop\Matlab\Labs Kevin Sterling Lab 4\Chapter8 PartiaMEX4.m Name EX4.m 2-load ('reaults.mat') %Write the velocity components into array of 250x51 uxx- reshape (ux, lx, ly) 5-1 uyy = reshape (uy, 1 x, ly) ; 6-dudy (250,51)-0.0 7dvdx (250,51)-0.0 Calculate the two derivatives in the whole domain 9- for i-1 : 250 for 1:51 if (i1&51) 13 elseif (1 &j51) dvdx (i,j)-(uyy (i+1,j)-uyy(i,j))/(itl - i)i elseif (250&&-1) 18- dudy (i.j)-(uxx (i,j+1)-uxx(i,5))/5+1-j) elseif (250&&J51) elseif (i-1&& j > 1 j 1 && j 1 &&i250&j51) dudy (1,j)(uxx (1,j) - uxx(i,j-1))j-j-1)) dvdx (,j)(uyy (+1,j)-uyy (i-1,j))/(i+1 - (-1)) dudy (i,j)(uxx (i,j+1)ux (i,j-1))/(j+1-j-1)) dvdx (i 1+1-(1-1)) Command Window MATLAB R2018b-academic use PLOTS EDITOR PUBLISHMEW Find Files Compare Sa Print GoTo Comment %23 Breakpoints RunRun and Advance Run and Advance New Open Save Find Time UsrsKevin SterlingDesktopMatlabLabsKevin Sterling Lab 4Chapter8 Partial Current Folder Editor-CAUsers Kevin Sterling\Desktop\Matlab\Labs Kevin Sterling Lab 4\Chapter8 PartiaMEX4.m Name EX4.m 2-load ('reaults.mat') %Write the velocity components into array of 250x51 uxx- reshape (ux, lx, ly) 5-1 uyy = reshape (uy, 1 x, ly) ; 6-dudy (250,51)-0.0 7dvdx (250,51)-0.0 Calculate the two derivatives in the whole domain 9- for i-1 : 250 for 1:51 if (i1&51) 13 elseif (1 &j51) dvdx (i,j)-(uyy (i+1,j)-uyy(i,j))/(itl - i)i elseif (250&&-1) 18- dudy (i.j)-(uxx (i,j+1)-uxx(i,5))/5+1-j) elseif (250&&J51) elseif (i-1&& j > 1 j 1 && j 1 &&i250&j51) dudy (1,j)(uxx (1,j) - uxx(i,j-1))j-j-1)) dvdx (,j)(uyy (+1,j)-uyy (i-1,j))/(i+1 - (-1)) dudy (i,j)(uxx (i,j+1)ux (i,j-1))/(j+1-j-1)) dvdx (i 1+1-(1-1)) Command Window

Step by Step Solution

There are 3 Steps involved in it

1 Expert Approved Answer

Step: 1 Unlock blur-text-image

Question Has Been Solved by an Expert!

Get step-by-step solutions from verified subject matter experts

Step: 2 Unlock

Step: 3 Unlock

Students Have Also Explored These Related Databases Questions!

Solve with MATLAB! Part 1: Concentration profile using finite difference approximation Equation 4 shows a changing concentration with respect to radial position as well as to time, making it a...

Diffusion out of capillary Part 1: Concentration profile using finite difference approximation Equation 4 shows a changing concentration with respect to radial position as well as to time, making it...

function = explicit(mx) %% Problem Number: 49 L = 45; tmax = 80; U = 0.8; D = 1.2; Xs = 4; Cs = 10; XFlux = 36; % % % % % % % m s m/s m2/s m kg/m m % Calculating dx and dt based of mx input dx =...

Make a MATLAB code that will answer Task 1 and Task 2. Computational Problem/Code 1 (Related to Project 1.3h, not really h, but the spirit of h with less work.) Implement equations/algorithms 1.17...

ECE 2240: Numerical Methods for Engineers Winter 2017 Lab 5: Numerical Differentiation and Integration Apr 3-7, 2017 The purpose of this lab is to: 1) Implement numerical methods for approximating...

Upload MAILAB code using link below Show that if you discretise the equation shown in @1 using the differentiation schemes above and implementing the boundary conditions, you will get a set of...

Q1.1 1 Point As shown in lectures, the derivative of a function y() can be calculated using the second order central, and first order forward and backward difference schemes (x) = (yi+1 - yi-1)/(24)...

assignment is to implement a shooting method to solve this nonlinear boundary value problem. In doing so, we make use of methods developed to solve nonlinear equations and initial value problems....

We would like to solve the steady state heat transfer problem given by the drawing in Fig. 7.17. We have a composite solid material of area L x L y , with the left half having conductivity k 1 and...

MATLAB CODE for above mentioned task I am getting errors in the area marked in bold. clear all; clc; close all; %% Define all coefficients for the 4-element Windkessel Model h = 0.01; % time step...

In Problem, find each indicated derivative and simplify [10, log t] dt

The following costs were incurred by Griff Co., a manufacturer, during 20x1: Accounting and legal fees Freight-in Freight-out Officers' salaries Insurance Sales representatives salaries 25,000...

X will pay the swap bank annual payments on $ 1 0 , 0 0 0 , 0 0 0 with the coupon rate of LIBOR - 0 . 1 5 % ; in exchange the swap bank will pay to company X interest payments on $ 1 0 , 0 0 0 , 0 0...

Seved Help 14 Wisconsin Snowmobile Corp. is considering a switch to level production Cost efficiencies would occur under level production, and aftertax costs would decline by $31,500, but inventory...

=+1 What makes managing employee health and safety programs around the world so difficult?

=+3 Describe the kinds of health, safety, and security problems that international employees encounter and design a crisis management program to deal with them.

=+4 Why is designing and implementing a global human resource information system so difficult? What kinds of problems need to be overcome?