MATLAB: How to plot using all the for loop values

integrationMATLAB

Hi,

Here I am pasting a simple code which does a simple integration from – infinity to + infinity:

syms x
C=1;D=1;E=1;N=1;W=1;
alp=1;
z=1;
for P=1:1:20
a=(sqrt(-1))*(C-P);
b=(D-(sqrt(-1))*(P+x));
c=(E+(sqrt(-1))*(N-P));
d=(W^2+x.^2);
e=b*c*d;
f=a/e;
S=alp*(int(f,-inf,inf));
out = vpa(S);
k=double(out)
end

******** The problem is to plot a graph (S Vs P) using all the values of P. Can anyone help me in this regard?

Thanking You!

Best Answer

You have to use: k(P) instead of k

syms x
C=1;D=1;E=1;N=1;W=1;
alp=1;
z=1;
for P=1:1:20
  aa=(sqrt(-1))*(C-P);
  b=(D-(sqrt(-1))*(P+x));
  c=(E+(sqrt(-1))*(N-P));
  d=(W^2+x.^2);
  e=b*c*d;
  f=a/e;
  S=alp*(int(f,-inf,inf));
    out = vpa(S);
    k(P)=double(out)
end
plot(abs(k))

Related Solutions

MATLAB: How to plot a 2D graph using all the for loop values

Your errors:

you can not use a negative index
when you write for w=-20:20, at the end of the loop w=20, is just one value
you should use plot (w,E) instead of plot(E,w)

alp=1;
K=1;
E0=1;
E=[];
for w=-20:1:20;
  a=K*(cos(pi*alp/2)*w^alp+w^(2*alp));
  b=1+2*(cos(pi*alp/2)*w^alp+w^(2*alp));
  c=a/b;
    E(end+1)=E0*(1+c);
end
figure(1)
w=-20:20
plot(w,E)
xlabel('w');
ylabel('E');
title('Plot of E vs w')t(E,w)

MATLAB: How to apply ‘parfor’ in this for-loop code

You need to rewrite the code in approximately the form:

      mu_11=zeros(Z_N,N_c-1,Z_N,N_c-1);
      nu_11=zeros(Z_N,N_c-1,Z_N,N_c-1);
      syms s_1 s_2
      for n=1:Z_N
          mu_11n = zeros(N_c-1,Z_N,N_c-1);
          nu_11n = zeros(N_c-1,Z_N,N_c-1);
          z_s1_cavitysiden = z_s1_cavityside(n,:);
          for k=1:N_c-1
              for m=0:Z_N-1 % Defining Functions
                  f=cot(pi*(z_s1_cavitysiden(k)-l(m+1)*s_1-zeta(m+1))*exp(-i*(pi/2-beta*pi/180))/Z_N/h-m*pi/Z_N);
                  g=s_1*cot(pi*(z_s1_cavitysiden(k)-l(m+1)*s_1-zeta(m+1))*exp(-i*(pi/2-beta*pi/180))/Z_N/h-m*pi/Z_N);
                  for p=1:N_c-1 % Perform Numerical Integration
                      mu_11n(k,(m+1),p)=vpaintegral(f,s_1,S_1(p),S_1(p+1),'RelTol',1e-4,'AbsTol',0,'MaxFunctionCalls',Inf);
                      nu_11n(k,(m+1),p)=vpaintegral(g,s_1,S_1(p),S_1(p+1),'RelTol',1e-4,'AbsTol',0,'MaxFunctionCalls',Inf);
                  end
              end
          end
          mu_11(n,:,:,:) = mu_11n;
          nu_11(n,:,:,:) = nu_11n;
          progress_11=100*n/Z_N
      end

It is also recommended to move calculations out of loops when possible. For example, -i*(pi/2-beta*pi/180))/Z_N/h does not change with n, k, or m, and could be pre-computed. Likewise, z_s1_cavitysiden(k) does not depend upon m, so the value could be computed inside the k loop instead of inside the m loop.

Best Answer

Related Solutions

MATLAB: How to plot a 2D graph using all the for loop values

MATLAB: How to apply ‘parfor’ in this for-loop code

Related Question