I want to plot Vcr as a function of m. In the solution of the program it is obtained that Vcr is
Vcr= 0 1 2 3 4
so something should appear, but I get this (nothing):
I don't know where and how I have to put the function plot in the code to obtain what I want. Here is the code:
close alls=input('Introduce s(1-40): ');n=input('Introduce n (1-100): ');t=1000;mf=zeros(1,s);Vf=zeros(s,n);cf=zeros(s,n,t);Vcr=zeros(1,s);Ccr=zeros(1,s);syms c V m;F=-m+V+c+i*(m-V-c);Re=real(F) Im=imag(F) for h=1:s mf(h)=h; m=mf(h); for j=1:n Vf(h,j)=j; V=Vf(h,j); PRe=subs(Re) ; d=0.01 ; %step size
D=20; Nvec=[1:d:D]; step=length(Nvec); for g=1:step cf(h,j,g)=g; c=cf(h,j,g); Prec=subs(PRe); if Prec>-0.0001 && Prec<0.0001 cr=c; break else cf(h,j,g)=cf(h,j,g); end end PImag=subs(Im) if PImag==-1900 Vcr(1,h)=V; Ccr(1,h)=c; else Vf(h,j)=Vf(h,j); end end
end
disp ('Las velocidades de migraciĆ³n son: ')Ccrdisp('las velocidades criticas son: ' )Vcrplot(m,Vcr)
Thank you for your help!
Best Answer