I have boundary conditions u(0,t)x = 0 and u(1,t) = 0. and initial conditions u(x,0) = u0(x) = −cos 2pix on the interval 0 ≤ x ≤ 1
So far I have created for the heat equation
I am now trying to create
function [pl,ql,pr,qr] = heatbc(xl,ul,xr,ur,t)
for the initial conditions but I dont understand how to plug in the boundary conditions for pl ql pr and qr.
Any help would be really appreciated.