Facing problem about plot the velocity profile

Respected sir i have attched one picture here , likethis picture i want to plot velocity(u) vs y plot

my formula for velocity profile is below whatever I have done in Matlab code I am writing below the formulae

U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))

  % t_o=106;
  t_e=156;
 deltavals = 0.6 ;
alphavals = [-0.3,- 0.2 ,-0.1,0, 0.1,0.2,0.3 ];
yvals = 0.3:0.1:1;
[ alpha, y] = ndgrid( alphavals, yvals);
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
plot(y,U,'-r') end

[Image Analyst]

Format your code please

Give us all the code, especially all the variable values that go into making up U. We can't run your code without it!

[Debasis Roy]

Respected sir value of variable t_o=100 t_e=150 mu=0.3 delta=0.6 alpha values are 0.1,0.2,0.3,0.3,0.4,-0.2,-0.5

[Image Analyst]

OK, now, what is y?

t_o = 100;
t_e = 150
mu = 0.3
delta = 0.6
alpha = [0.1, 0.2, 0.3, 0.3, 0.4, -0.2, -0.5]
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
% t_o = 106;
t_e = 156;
deltavals = 0.6 ;
alphavals = [-0.3,- 0.2 ,-0.1,0, 0.1,0.2,0.3 ];
yvals = 0.3:0.1:1;
[ alpha, y] = ndgrid( alphavals, yvals);
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
plot(y,U,'-r')

If this is still a problem for you, give us everything we need to help you.

[Souvick roy]

sir due to problem in value i have non dimensionalise my final equation. i have attached the picture. in wquation one is valid from y =0 to y=0.3 and equation two is valid for y=0.3 to 1

I want to divide y in 0.01 interval

balue of delta 0.3

value of k=100

value of beta 0.01

please help me in wriying the programme and plotting