Hi, I need to write a function U like this :
if x=0 then U(0,y)=0
if y=0 then U(x,0)=0
if y=10 then U(x,10)=sin(pi*x/L).*sin(pi*10/L) for x=[0,L/2]
if x=10 then U(10,y)=sin(pi*y/L).*sin(pi*10/L) for y=[0,L/2]
if y=5 then U(x,5)=sin(pi*x/L).*sin(pi*5/L) for x=[L/2,L]
if x=5 then U(5,y)=sin(pi*y/L).*sin(pi*5/L) for y=[L/2,L]
So, i have a vector of x and y, for example x=[0 0 0 0 1 ….5….10..0..1..5….0] and y=[0 0 0 0 1 ….5….10..0..1..5….0], and npt=size(x)=size(y). I did code about it but it not work
function z= U(npt,x,y,L) for m=1:npt if( x(m)==0) z=0; end end for m=1:npt if x(m)==L z= sin(pi*y/L).*sin(pi*10/L); end end for m=1:npt if y(m)==0 z=0; end end for m=1:npt if y(m)==L z = sin(pi*x/L).*sin(pi*10/L); end end for m=1:npt if x(m)==L/2 for y(m)=L/2:L z= sin(pi*5/L).*sin(pi*y/L); end end end for m=1:npt if y(m)==L/2 for x(m)=L/2:L z = sin(pi*5/L).*sin(pi*x/L); end end end end
Best Answer