MATLAB: How to calcul equation using function Matlab

functionif statementMATLABmatlab function

Hi,
I try to write a function for calcul eqution u(x,y t)=exp(-t)*sin(pi*x/L)*sin(pi*y/L) on domaine L-shaped
function z=U(x,y,time,L)
x=[0;0;0;0;0;0;0;0;0;0;0;1;2;3;4;5;6;7;8;9;10;10;10;10;10;10;9;8;7;6;5;5;5;5;5;5;4;3;2;1;1;2;3;4;5;6;7;8;9;1;2;3;4;5;6;7;8;9;1;2;3;4;5;6;7;8;9;1;2;3;4;5;6;7;8;9;1;2;3;4;1;2;3;4;1;2;3;4;1;2;3;4;1;2;3;4];
y=[0;1;2;3;4;5;6;7;8;9;10;0;0;0;0;0;0;0;0;0;0;1;2;3;4;5;5;5;5;5;5;6;7;8;9;10;10;10;10;10;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;5;5;5;5;6;6;6;6;7;7;7;7;8;8;8;8;9;9;9;9];
z=zeros(size(x));
time=0.15; L=10; E=exp(-time)
%gamma1
 if(x==0 & y>=0&y<=L)
     z=E.*sin(pi*y(x==0&y>=0&y<=L)/L).*sin(pi*0/L);
 
 %gamma2 
elseif( y==L & x>0&x<=L/2)
 z=E.*sin(pi*x(y==L&x>0&x<=L/2)/L).*sin(pi*L/L); 
%gamma3 
elseif( x==L/2 & y>=L/2&y<L)
 z=E.*sin(pi*y(x==L/2&y>=L/2&y<L)/L).*sin(pi*5/L);
 
%gamma4
 elseif( y==L/2 & x>L/2&x<=L)
 z=E.*sin(pi*x(y==L/2&x>L/2&x<=L)/L).*sin(pi*5/L)
 
%gamma5
 elseif (x==L & y>=0&y<L/2) 
z=E.*sin(pi*y(x==L&y>=0&y<L/2)/L).*sin(pi*10/L);
 
 %gamma6 
elseif( y==0 & x>0&x<L)
 z=E.*sin(pi*x(y==0&x>0&x<L)).*sin(pi*0/L);
 
%gamma7
 % elseif ??? I don't how to  write in this case !
z=E.*sin(pi*x/L).*sin(pi*y/L); 
end
end
Please Help me.
Thanks

Best Answer

function z=U(x,y,time,L)
x=[0;0;0;0;0;0;0;0;0;0;0;1;2;3;4;5;6;7;8;9;10;10;10;10;10;10;9;8;7;6;5;5;5;5;5;5;4;3;2;1;1;2;3;4;5;6;7;8;9;1;2;3;4;5;6;7;8;9;1;2;3;4;5;6;7;8;9;1;2;3;4;5;6;7;8;9;1;2;3;4;1;2;3;4;1;2;3;4;1;2;3;4;1;2;3;4];
y=[0;1;2;3;4;5;6;7;8;9;10;0;0;0;0;0;0;0;0;0;0;1;2;3;4;5;5;5;5;5;5;6;7;8;9;10;10;10;10;10;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;5;5;5;5;6;6;6;6;7;7;7;7;8;8;8;8;9;9;9;9];
z=zeros(size(x));
time=0.15; L=10; E=exp(-time)
%gamma1
 mask = (x==0 & y>=0&y<=L)
 z(mask) = E.*sin(pi*y(mask)/L).*sin(pi*0/L);
 
 %gamma2 
mask = ( y==L & x>0&x<=L/2)
 z(mask) = E.*sin(pi*x(mask)/L).*sin(pi*L/L); 
%gamma3 
mask = ( x==L/2 & y>=L/2&y<L)
z = E.*sin(pi*y(mask)/L).*sin(pi*5/L);
and so on.
Related Question