lab5q2exp.m

algrithm for optimization matlab source code for heat conducting plate 3 methods are included

clc;
close all;
clear all;

n = 19;
m = 400;
beta = 1;

dx = (2*pi)/(n+1);
x=0;
xmax = 2*pi;

dt = 10/m;
t=0;
tmax = 10;

lamda = (beta*dt)/(dx^2);
T = zeros((tmax/dt),(xmax/dx)+1);
i = 1;
j = 1;

for x = 0:dx:xmax
    T(1,j) = sin(x);
    j = j+1;
end

for i = 1:tmax/dt
    for j = 2:xmax/dx
        T(i+1,j)= lamda.*(T(i,j-1)+T(i,j+1)+(1-2*lamda).*T(i,j));
    end
end

X = 0:dx:xmax;
Y = 0:dt:tmax;
%mesh (X:Y:T);
figure
subplot(3,1,1)
plot(Y,T(:,2));
%AXIS([0 10 0 1]);

subplot(3,1,2)
plot(Y,T(:,5));
%AXIS([0 10 0 1]);

subplot(3,1,3)
plot(Y,T(:,11));
%AXIS([0 10 0 1]);


%Tdeltax1    = B(:,2);
%Tdeltax4    = B(:,5);
%Tdeltax10   = B(:,11);