problem but nobody seems to be able to explain or solve it in plain English. Could somebody please explain what this error is and how to fix it?
I have a simple function:
function [r]=Mec134function[w,theta_deg]
t2=10000;
theta_rad=[theta_deg./180].*pi;
t1=55090./[10*sin[theta_rad]];
rx=[t1.*cos[theta_rad]]-t2;
ry=w-[t1.*sin[theta_rad]];
r=[[rx].^2+[ry].^2].^0.5;
end
It seems to give this error for line 3 but I'm not sure why.
Thanks.
Accepted Answer
Your function defines 2 input arguments [w and theta_deg]. When you run Mec134function, you must specify exactly two inputs, otherwise you will get the error "Not enough input arguments".
For example, if you run the Mec134function in the command window without specifying any arguments:
You get this error:
Not enough input arguments.
Error in Mec134function [line 3]
theta_rad=[theta_deg./180].*pi;
If you run the Mec134function and specify two input arguments, "w" and "theta_deg" [assuming "w" and "theta_deg" are defined], you do not get the error message:
\>> Mec134function[w,theta_deg]
If you have the file "Mec134function.m" open in the Editor and you try to run the function by pressing the "Run" button or F5, MATLAB runs the Mec134function without any input arguments, and you get the error "Not enough input arguments". The "Run" button dropdown menu then opens prompting you to enter values for the missing input arguments.
Add the desired values and press enter. The values you enter are set as the default inputs when you click the "Run" button or F5 in the future.
To change the values, press the down arrow below the "Run" button and enter new values.
More Answers [19]
sorry the function is
function [r]=Mec134function[w,theta_deg]
t2=10000;
theta_rad=[theta_deg./180].*pi;
t1=55090./[10*sin[theta_rad]];
rx=[t1.*cos[theta_rad]]-t2;
ry=w-[t1.*sin[theta_rad]];
r=[[rx].^2+[ry].^2].^0.5;
end
if that's clearer :-]
Thanks :-]
will have a look at the "getting started"
I have a simple function:
function [r]=Mec134function[w,theta_deg]
t2=10000;
theta_rad=[theta_deg./180].*pi;
t1=55090./[10*sin[theta_rad]];
rx=[t1.*cos[theta_rad]]-t2;
ry=w-[t1.*sin[theta_rad]];
r=[[rx].^2+[ry].^2].^0.5;
end
that seems to give this error for line 2 but I'm not sure why.
I too am very new to Matlab, and tried to run the code above in .m [JP Donlon on Nov. 7th]. However, I keep getting an error stating "Not enough input arguments." I'm not sure what this means, because I have attempted to run other code by professors which works on other computers. Is it something with my preference settings?
Also, when I run the Code Analyzer, there are no issues...not sure what is going on.
How to rectify this error
Error using DetectFace [line 68]
Not enough input arguments.
the code is given below
I=imread['lena.jpg'];
minFace = 20;
maxFace = 4000;
overlappingThreshold = 0.5;
numThreads = 24;
if nargin > 2 && ~isempty[options]
if isfield[options, 'minFace'] && ~isempty[options.minFace]
minFace = options.minFace;
end
if isfield[options, 'maxFace'] && ~isempty[options.maxFace]
maxFace = options.maxFace;
end
if isfield[options, 'overlappingThreshold'] && ~isempty[options.overlappingThreshold]
overlappingThreshold = options.overlappingThreshold;
end
if isfield[options, 'numThreads'] && ~isempty[options.numThreads]
numThreads = options.numThreads;
end
end
if ~ismatrix[I]
I = rgb2gray[I];
end
candi_rects = NPDScan[model, I, minFace, maxFace, numThreads];
if isempty[candi_rects]
rects = [];
return;
end
clc; clear; close all;
xo=0.4;
A=[];
b=[];
Aeq=[];
beq=[];
Q=100;
R=1;
N = 50;
U0= zeros[100,1];
x=xo; h = 0.1;
xo=[-0.5,0.5];
options = optimoptions[@fmincon,'Algorithm','sqp'];
U = fmincon[@cost1,U0,[],[],[],[],[],[],@confuneq,options];
for k =1:N
S1= F[x[k],U[k]];
S2=F[x[k]+0.5*h*S1,U[k]];
S3=F[x[k]+0.5*h*S2,U[k]];
S4=F[x[k]+h*S3,U[k]];
x[k+1] = x[k] + [1/6]* [S1+ 2*S2+ 2*S3 + S4]*h;
k = k + 1 ;
end
plot[U];
grid on
figure[]; plot[x]
grid on
I new to matlab.Can anyone help me to fix the iisue here.
Its showing not enough input arguments.
Hello sir,I am newer to matlab. I am getting error in this code like "preprocessing requries more input arugument to run". can you please to run this program. my project topic is recognition and matching fake logos using filters.
function img=preprocessing[I]
[x y o]=size[I];
if o==3
I=rgb2gray[I];
end
I=im2double[I];
med=medfilt2[I];
figure,imshow[med]
title['MEDIAN FILTERED IMAGE']
[psnr_med,mse_med]=psnr[I,med]
out7= imfilter[I, fspecial['average']];
figure,imshow[out7]
title['MEAN FILTERED IMAGE']
[psnr_avg,mse_avg]=psnr[I,out7]
gau=imfilter[I,fspecial['gaussian']];
figure,imshow[gau]
title['GAUSSIAN FILTER IMAGE']
[psnr_gau,mse_avg]=psnr[I,gau]
psf=fspecial['gaussian',7,10];
image1=imfilter[I,psf,'conv','circular'];
var1=[1/256]^2/12;
var2=var[I[:]];
wei=deconvwnr[image1,psf,[var1/var2]];
figure,imshow[wei];title['WEINER FILTERED IMAGE'];
[psnr_wei,mse_wei]=psnr[I,wei]
psnr_all=[psnr_med,psnr_avg,psnr_gau,psnr_wei];
psnr_max=max[psnr_all];
val=find[psnr_all==psnr_max];
if val==1
img=med;
disp['median have high psnr'];
elseif val==2
img=out7;
disp['mean have high psnr'];
elseif val==3
img=gau;
disp['gaussian have high psnr'];
else
img=wei;
disp['weiner have high psnr'];
end
I am getting error for below function as "Not enough input arguments. "
delayed_signal = mtapped_delay_fcn[input];
Hi I'm trying to run Dr. John Stockie's matlab code but I am getting a "Not enough input argument" error. I'm not very well verse with Matlab, so I would appreciate any help...Thank you. I am pasting the code:
function C = ermak[ x, y, z, H, Q, U, Wset, Wdep ]
Umin = 0.0;
ay = 0.34; by = 0.82; az = 0.275; bz = 0.82;
sigmay = ay*abs[x].^by .* [x > 0];
sigmaz = az*abs[x].^bz .* [x > 0];
Kz = 0.5*az*bz*U*abs[x].^[bz-1] .* [x > 0];
if U < Umin,
C = 0 * z;
else
Wo = Wdep - 0.5*Wset;
C = Q ./ [2*pi*U*sigmay.*sigmaz] .* exp[ -0.5*y.^2./sigmay.^2 ] .* ...
exp[ -0.5*Wset*[z-H]./Kz - Wset^2*sigmaz.^2/8./Kz.^2 ] .* ...
[ exp[ -0.5*[z-H].^2./sigmaz.^2 ] + ...
exp[ -0.5*[z+H].^2./sigmaz.^2 ] - sqrt[2*pi]*Wo*sigmaz./Kz .* ...
exp[ Wo*[z+H]./Kz + 0.5*Wo^2*sigmaz.^2./Kz.^2 ] .* ...
erfc[ Wo*sigmaz/sqrt[2]./Kz + [z+H]./sqrt[2]./sigmaz ] ];
ii = find[isnan[C] | isinf[C]];
C[ii] = 0;
end
and the error message refers to "sigmay" in line 31
function test[num1, num2,small,s]
load [small, num1, num2]
s = sum[num1, num2]
end
this the code for this i'm getting these errors
Not enough input arguments.
Error in test [line 3]
load [small, num1, num2]
Hello. Am new in Matlab and I want to do my assignment with this function refraction_2layers and Matlab is saying error using refraction_2layers [line 18 and 25] Here is the whole program
function refraction_2layers[v1, v2, z, FIRST_ARRIVALS_ONLY];
if nargin < 4 FIRST_ARRIVALS_ONLY = 0; end
x = [0:5:300];
t1 = x./v1;
t2 = [2*z*sqrt[v2^2-v1^2]/[v1*v2]]+x./v2;
xcrit = 2*z*v1/[sqrt[v2^2-v1^2]];
if isreal[xcrit]
a = min[find[x>xcrit]];
end
crossover = [[2*z*sqrt[v2^2-v1^2]]/[v1*v2]]/[1/v1-1/v2];
b = max[find[x> project2
Error using fzero [line 306]
FZERO cannot continue because user-supplied function_handle ==> temp failed with the error below.
Not enough input arguments.
Error in project2>Moles [line 30]
T = fzero[@temp, T0, options];
Error in odearguments [line 90]
f0 = feval[ode,t0,y0,args{:}];
Error in ode45 [line 115]
odearguments[FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin];
Error in project2 [line 7]
[xt, L] = ode45[@Moles, xtspan, L0];
hi
I am new to matlab. i am getting error message "extract_features" requires more input arguments to run.
anderror in command window :
\>> Extract_Features
Not enough input arguments.
Error in Extract_Features [line 2]
img1 = imread[filename];
code is written below:
function Extract_Features[filename,flag]
if ndims[img1] == 3; img1 = rgb2gray[img1]; end
disp[['Extracting features from ' filename ' ...']];
plot[fir[fir1,1],fir[fir1,2],'r+'];
plot[fir[fir3,1],fir[fir3,2],'bo'];
filename2=filename; filename2[end-1]='x'; filename2[end]='t';
save[filename2,'fir','-ascii'];
Can someone help me please
envelope = sqrt[movmean[rec_EMG.^2], 'window'];
I was trying to do get the RMS but it says:
Error using movmean Not enough input arguments.
I didn't understand that, as I already use two arguments there.
Thanks fo ryour help
Can any one help me this please
function u = Anti_Tv[g,my,gamma]
gHS = uint8[imadjust [g]];
gGC = uint8[255.*[[double[g]./255].^[gamma]]];
g = double[g[:]];
n = length[g];
b = zeros[2*n,1];
d = b;
u = g;
eer = 1;k = 1;
tol = 1e-3;
Lambda = 0.05
[B, Bt, BtB] = DiffOper[sqrt[n]];
Not enough input arguments
Dear Matlab experts, If anyone one of you would like to assist me running the below code i would be really greatfule to you.
The error i am getting.
qardlecm
Not enough input arguments.
Error in qardlecm [line 24]
nn = size[data,1];
The code i want to run
....................................................................................................................................................................................................
function[bigphia,bigpia,thett,distthett] = qardlecm[data,ppp,qqq,tau]
pd = makedist['normal','mu',0,'sigma',1];
hb[jj,1] = [4.5*normpdf[icdf[pd,tau[jj,1]]]^4/[nn*[2*icdf[pd,tau[jj,1]]^2+1]^2]]^0.2;
hs[jj,1] = za^[2/3]*[1.5*normpdf[icdf[pd,tau[jj,1]]]2/[nn*[2*icdf[pd,tau[jj,1]]^2+1]]][1/3];
xx = data[:,2:size[data,2]];
ee = xx[2:nn,:] - xx[1:[nn-1],:];
eei = zeros[nn-qqq,qqq*k0];
eei[:, ii+1+[jj-1]*qqq] = ee[[qqq+1-ii]:[nn-ii],jj];
yyi[:,ii] = yy[[1+ppp-ii]:[nn-ii],1];
X = [eei[[size[eei,1]+1-size[yyi,1]]:size[eei,1],:], xxi[[size[xxi,1]+1-size[yyi,1]]:size[xxi,1],:], yyi];
X = [eei, xxi, yyi[[size[yyi,1]+1-size[xxi,1]]:size[yyi,1],:]];
ONEX = [ones[size[X,1],1],X];
Y = yy[[nn-size[X,1]+1]:nn,1];
bt = zeros[size[ONEX,2],ss];
[bt1] = qregressMatlab[Y,ONEX,tau[jj,1]];
fh[jj,1] = mean[normpdf[-uu[:,jj]/hb[jj,1]]]/hb[jj,1];
barw = zeros[nn-1,qqq*k0];
barw[jj:[nn-1],[k0*[jj-1]+1]:k0*jj] = ee[2:[nn-jj+1],:];
tw = [ones[nn-1,1], barw];
mm = [xx[[qqq+1]:nn,:]'*xx[[qqq+1]:nn,:] - xx[[qqq+1]:nn,:]'*tw[qqq:[nn-1],:]*inv[tw[qqq:[nn-1],:]'*tw[qqq:[nn-1],:]]*tw[qqq:[nn-1],:]'*xx[[qqq+1]:nn,:]]/[nn-qqq]^2;
bb[jj,1] = 1/[[1-sum[bt[2+[qqq+1]*k0:1+[qqq+1]*k0+ppp,jj],1]']*fh[jj,1]];
qq[jj,ii] = [min[psu,[],1]' - tau[jj,1]*tau[ii,1]]*bb[jj,1]*bb[ii,1];
midbt[:,jj] = bt[2+qqq*k0:1+[qqq+1]*k0,jj]/[1-sum[bt[2+[qqq+1]*k0:1+[qqq+1]*k0+ppp,jj],1]'];
bigbt = reshape[midbt,[],1];
bigbtmm = kron[qq,inv[mm]];
wwj = zeros[nn-ppp,qqq*k0];
yyj[:,jj] = yy[[ppp+1-jj]:[nn-jj],1];
wwj[:,jj+[ii-1]*qqq] = ee[[ppp-jj+2]:[nn-jj+1],ii];
kk = zeros[nn-ppp,ss*ppp];
ONEX = [ones[nn-ppp,1],xxj,wwj];
[bbt] = qregressMatlab[Y,ONEX,tau[ii,1]];
kk[:,jj+[ii-1]*ppp] = kkk;
llla = [kka'*kka - kka'*tilw*inv[tilw'*tilw]*tilw'*kka]/[nn-ppp];
wwj = zeros[nn-qqq,qqq*k0];
yyj[:,jj] = yy[[qqq+1-jj]:[nn-jj],1];
wwj[:,jj+[ii-1]*qqq] = ee[[qqq-jj+2]:[nn-jj+1],ii];
kk = zeros[nn-qqq,ss*ppp];
ONEX = [ones[nn-qqq,1], xxj, wwj];
[bbt] = qregressMatlab[Y,ONEX,tau[jj,1]];
kk[:,jj+[ii-1]*ppp] = kkk;
llla = [kka'*kka - kka'*tilw*inv[tilw'*tilw]*tilw'*kka]/[nn-qqq];
cc[jj,ii] = [min[psu,[],1]' - tau[jj,1]*tau[ii,1]]/[fh[ii,1]*fh[jj,1]];
bigpia = zeros[ss*[ppp-1],ss*[ppp-1]];
psu = inv[llla[[jj-1]*[ppp-1]+1:jj*[ppp-1],[jj-1]*[ppp-1]+1:jj*[ppp-1]]]*llla[[jj-1]*[ppp-1]+1:jj*[ppp-1],[ii-1]*[ppp-1]+1:ii*[ppp-1]]*inv[llla[[ii-1]*[ppp-1]+1:ii*[ppp-1],[ii-1]*[ppp-1]+1:ii*[ppp-1]]];
bigpia[[jj-1]*[ppp-1]+1:jj*[ppp-1],[ii-1]*[ppp-1]+1:ii*[ppp-1]] = cc[jj,ii]*psu;
midphi[:,jj] = bt[2+[qqq+1]*k0:1+[qqq+1]*k0+ppp,jj];
bigphi = reshape[midphi,[],1];
bigphia = [bigphia1; bigphia2; bigphia3];
dg = [nn^[1/2],0,0; 0,nn^[1/2],0; 0,0,nn];
r2 = sum[tilwb,1]*[nn-2]^[-1];
r3 = sum[xx[3:nn,1],1]*[nn-2]^[-3/2];
rh9 = xx[3:nn,1]'*xx[3:nn,1];
QQQ = [r1, r2, r3 ; r4, r5, r6; r7, r8, r9];
psiu[rr,jj] = tau[jj,1]-1;
sigmma = psiu'*psiu*[1/[nn-2]];
sigma1 = mean[psiu1.^[2]];
sigma2 = mean[psiu2.^[2]];
sigma3 = mean[psiu3.^[2]];
distmt1 = nn*fh[1,1]^[-2]*sigmma[1,1]*inv[dg]*inv[QQQ]*inv[dg];
distmt2 = nn*fh[1,1][-1]*fh[2,1][-1]*sigmma[1,2]*inv[dg]*inv[QQQ]*inv[dg];
distmt3 = nn*fh[1,1][-1]*fh[3,1][-1]*sigmma[1,3]*inv[dg]*inv[QQQ]*inv[dg];
distmt4 = nn*fh[2,1][-1]*fh[1,1][-1]*sigmma[2,1]*inv[dg]*inv[QQQ]*inv[dg];
distmt5 = nn*fh[2,1]^[-2]*sigmma[2,2]*inv[dg]*inv[QQQ]*inv[dg];
distmt6 = nn*fh[2,1][-1]*fh[3,1][-1]*sigmma[2,3]*inv[dg]*inv[QQQ]*inv[dg];
distmt7 = nn*fh[3,1][-1]*fh[1,1][-1]*sigmma[3,1]*inv[dg]*inv[QQQ]*inv[dg];
distmt8 = nn*fh[3,1][-1]*fh[2,1][-1]*sigmma[3,2]*inv[dg]*inv[QQQ]*inv[dg];
distmt9 = nn*fh[3,1]^[-2]*sigmma[3,3]*inv[dg]*inv[QQQ]*inv[dg];
distcon1 = A11 + A12 + A13;
distcon2 = A21 + A22 + A23;
distcon3 = A31 + A32 + A33;
distcon4 = A41 + A42 + A43;
distcon5 = A51 + A52 + A53;
distcon6 = A61 + A62 + A63;
distcon7 = A71 + A72 + A73;
distcon8 = A81 + A82 + A83;
distcon9 = A91 + A92 + A93;
distthett = [distcon1, distcon2, distcon3 ; distcon4, distcon5, distcon6 ; distcon7, distcon8, distcon9];
thett1 = bt[2,1] + bt[3,1];
thett2 = bt[2,2] + bt[3,2];
thett3 = bt[2,3] + bt[3,3];
thett = [thett1 ; thett2 ; thett3];
Edited: DGM on 22 Feb 2023
I get Unrecognized function or variable 'theta'.
how to fix it
function hpol =polar_dB[theta,rho,rmin,rmax,rticks,line_style]
error['Requires 5 or 6 input arguments.']
theta = th[:,ones[1,nr]];
theta = th[:,ones[1,nr]];
if isstr[theta] || isstr[rho]
error['Input arguments must be numeric.'];
if any[size[theta] ~= size[rho]]
error['THETA and RHO must be the same size.'];
next = lower[get[cax,'NextPlot']];
fAngle = get[cax, 'DefaultTextFontAngle'];
fName = get[cax, 'DefaultTextFontName'];
fSize = get[cax, 'DefaultTextFontSize'];
fWeight = get[cax, 'DefaultTextFontWeight'];
set[cax, 'DefaultTextFontAngle', get[cax, 'FontAngle'], ...
'DefaultTextFontName', font_name, ...
'DefaultTextFontSize', font_size, ...
'DefaultTextFontWeight', get[cax, 'FontWeight'] ]
hhh=plot[[0 max[theta[:]]],[min[rho[:]] max[rho[:]]]];
v = [get[cax,'xlim'] get[cax,'ylim']];
ticks = length[get[cax,'ytick']];
elseif rem[rticks,3] == 0
inds = 1:[length[th]-1]/4:length[th];
xunits[inds[2:2:4]] = zeros[2,1];
yunits[inds[1:2:5]] = zeros[3,1];
rinc = [rmax-rmin]/rticks;
for i=[rmin+rinc]:rinc:rmax
plot[xunit*is,yunit*is,'-','color',tc,'linewidth',0.5];
text[0,is+rinc/20,[' ' num2str[i]],'verticalalignment','bottom' ];
cst = cos[th]; snt = sin[th];
plot[[rmax-rmin]*cs,[rmax-rmin]*sn,'-','color',tc,'linewidth',0.5];
cst2 = cos[th2]; snt2 = sin[th2];
cs2 = [[rmax-rmin-george]*cst2; [rmax-rmin]*cst2];
sn2 = [[rmax-rmin-george]*snt2; [rmax-rmin]*snt2];
plot[cs2,sn2,'-','color',tc,'linewidth',0.15];
plot[-cs2,-sn2,'-','color',tc,'linewidth',0.15];
text[rt*cst[i],rt*snt[i],int2str[abs[i*30-90]],'horizontalalignment','center' ];
loc = int2str[180-[i*30+90-180]];
text[-rt*cst[i],-rt*snt[i],loc,'horizontalalignment','center' ];
axis[[rmax-rmin]*[-1 1 -1.1 1.1]];
set[cax, 'DefaultTextFontAngle', fAngle , ...
'DefaultTextFontName', font_name, ...
'DefaultTextFontSize', fSize, ...
'DefaultTextFontWeight', fWeight ];
if theta[i]*180/pi >=0 && theta[i]*180/pi =90
xx[i] = [rho[i]-rmin]*cos[-theta[i]+pi/2];
yy[i] = [rho[i]-rmin]*sin[-theta[i]+pi/2];
elseif theta[i]*180/pi < 0
xx[i] = [rho[i]-rmin]*cos[abs[theta[i]]+pi/2];
yy[i] = [rho[i]-rmin]*sin[abs[theta[i]]+pi/2];
if strcmp[line_style,'auto']
q = plot[xx,yy,line_style];
axis['equal'];axis['off'];
if ~hold_state, set[cax,'NextPlot',next]; end
Dear Matlab expert, please help me. My Matlab is 2016a, I try to run a code but there is an error "Not enough input arguments."
Error in astar [line 3]
ssNode = startNode;
function [ClosedList,cost,heuristic,func,iteration] = astar[source,target,weights,heuristics,startNode,goalNode]
[s,t,n,sNode,gNode] = refactor[source,target,weights,sNode,ggNode];
ClosedList = struct['Path' ,sNode,'Cost',0,'Heuristic',heuristics[sNode],'F',heuristics[sNode]];
OpenList = [OpenList ClosedList];
while[isGoalReached[OpenList,gNode]==0 && ~isempty[OpenList]] [minI,minP] = minPath[OpenList];
newPaths = getNewPaths[s,t,weights,heuristics,minP]; OpenList = [OpenList newPaths];
[~,minP] = minPath[OpenList];
ClosedList = n[minP.Path];
heuristic = minP.Heuristic;
function [minIndex,ClosedList] = minPath[paths]
ClosedList = paths[minIndex];
if[paths[i].F < ClosedList.F]
ClosedList = paths[minIndex];
function isGoal = isGoalReached[paths,goalNode]
[~,minP] = minPath[paths];
if[minP.Path[length[minP.Path]] == goalNode]
function weight = getWeight[s,t,weights,nodeA,nodeB]
if[s[i]==nodeA && t[i]==nodeB]
function paths = getNewPaths[s,t,w,h,path]
uniqueNodes = getNodes[s,t];
currentNode = path.Path[length[path.Path]]; childs = getChilds[s,t,currentNode];
if[isempty[find[path.Path==childs[i], 1]]]
c = path.Cost + getWeight[s,t,w,currentNode,childs[i]];
heur = h[uniqueNodes==childs[i]];
p = struct['Path',[path.Path childs[i]],'Cost',c,'Heuristic',heur,'F',f];
function childs = getChilds[source,target,node]
childs = sort[target[source==node]];
function nodes = getNodes[s,t]
nodes = unique[horzcat[s,t]];
function [s,t,n,sn,gn] = refactor[source,target,~,startNode,goalNode]
uNodes = unique[horzcat[source,target]]; n = uNodes;
uNodes = unique[horzcat[source,target]];
[~,sIndex] = ismember[source[i],uNodes];
[~,tIndex] = ismember[target[i],uNodes];
Sol this programm
function I = Int[rho, z, t, fdphi]
global D lam alpha0 L0 Ep tau gamma M2 f w;
gamma = n2 * lam ./ [2 * pi];
w0 = [1e-3 * M2 * lam * f] ./ [pi * w];
wz = w0 * sqrt[1 + [z ./ z0]^2];
Rz = z * [1 + [z0 / z]^2];
r = linspace[0, wz, 100];
Leff = L * [1 - exp[-alpha0 * L0]] / alpha0;
dphi = k * Leff * gamma * Ep * 1e21 * exp[-2 * r.^2 ./ wz.^2] .* exp[-4 * log[2 * t^2 / tau^2]] ./ [1 + [z / z0]^2];
Iint = Iint - 1i * k * r.^2 / [2 * Rz];
Isec = k * rho[m] * r / D;
Isec = besselj[0, Isec] .* Iint;
Isec = trapz[r, Isec] * dr;
Iprim = [4 * pi^2 * Ep] / [lam * D]^2 * [w0 / wz]^2 * exp[-alpha0 * L0];
Iprim = Iprim * exp[-4 * log[2] * t^2 / tau^2];
xlabel['Radial Coordinate [mm]'];
title['Far Field Intensity'];
rho = linspace[-15/2, 15/2, 100];
xlabel['Radial Coordinate [mm]'];
title['Far Field Intensity'];
Np=100;a=le-10;X=a*[1:1:Np];t0=[hbar^2]/[2*m*[a^2]]/q;L=a*[Np+1];
T=[2*t0*diag[ones[1,Np]]]-[t0*diag[ones[1,Np-1],1]]-[t0*diag[ones[1,Np-1],-1]];
[V,D]=eig[T];D=diag[D];[Enum,ind]=sort[D];
El=D[ind[1]];psil=abs[V[:,ind[1]]];P1=psil.*conj[psil];
E2=D[ind[50]];psi2=abs[V[:,ind[50]]];P2=psi2.*conj[psi2];
Ean=[[[hbar*pi]^2]/[2*m*[L^2]]/q]*[1:Np].*[1:Np];
xlabel[' Eigenvalue Number,alpha-->'];
See Also
Categories
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
An Error Occurred
Unable to complete the action because of changes made to the page. Reload the page to see its updated state.