how to calculate a difference of function with any parameter

function y=f(x)

for i=10

y=i*sin(x)

end

end

calculate the change in y/change in x vs i

Answer 1

[Torsten]

function main
  h = 1e-6;
  xleft = -pi;
  xright = pi;
  n = 200;
  dx = (xright-xleft)/n;
  x = xleft:dx:xright;
  change_in_y_by_change_in_x = (f(x+h)-f(x))/h;
  plot(x,change_in_y_by_change_in_x)
  hold on
  plot(x,cos(x))    % for comparison, plot (sin(x))' = cos(x)
end

function y = f(x)
  y = sin(x)
end

[shiv gaur]

pl plot x vs dy/dx

[Torsten]

Done.

[shiv gaur]

Dr torsten you help me for the earlier program i stuck at some point pl help me to solve the same

[shiv gaur]

Dr torsten this the program you help me for solving equation I am little bit stuck in program pl help

program is below

function shiv7

%K0= (2*pi/1530)*1e6:(2*pi/10)*1e6:(2*pi/1630)*1e6;

%K0= linspace((2*pi/1530)*1e6,(2*pi/1630)*1e6,5);

%D1=1e-9:1e-9:3e-7;

D1= linspace(1e-9,3e-7,15);

for j=1:numel(D1)

da = D1(j);

p0 = 1;

p1 = 1.5;

p2 = 2;

TOL = 10^-8;

N0 = 100; format long

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=3;

while i <= N0

b = DELTA2 + h2*d;

D = (b^2 - 4*f(p2,da )*d)^(1/2);

if abs(b-D) < abs(b+D)

E = b + D;

else

E = b - D;

end

h = -2*f(p2,da)/E;

p = p2 + h;

if abs(h) < TOL

disp(p)

break

end

p0 = p1;

p1 = p2;

p2 = p;

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da ) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=i+1;

end

%pv(j)=p;

P(j)=real(p)

end

plot(D1,P)

end

function y=f(x,da )

l=6330*(10)^(-10);

k0=2*pi/l;

nc=1.33;

ec=nc^2;

na=1.63;

ea=na^2;

nf=1.59;

ef=nf^2;

nm=0.064+1i*4;

em=nm^2;

ns=1.495;

es=ns^2;

df=0.5e-6;

dm=0.03e-6;

kc=k0*sqrt(x^2-ec);

ka=k0*sqrt(ea-x^2);

kf=k0*sqrt(ef-x^2);

km=k0*sqrt(em-x^2);

ks=k0*sqrt(x^2-es);

rfm=(kf-km)/(kf+km);

rms=(km-ks)/(ks+km);

rfa=(kf-ka)/(kf+ka);

rac=(ka-kc)/(ka+kc);

a=(1-rfm)/(1+rfm);

b=(1-rms*exp(2*1i*dm*km))/(1+rms*exp(2*1i*dm*km));

c=(1-rfa)/(1+rfa);

f=(1-rac*exp(2*1i*da*ka))/(1+rac*exp(2*1i*da*ka));

pfms=2*atan(1i*a*b);

pfac=2*atan(1i*c*f);

y=2*kf*df-pfms-pfac;

end

we have to same find out the delta(y)/delta(na) vs da plot thanks

[shiv gaur]

plot is like this

[Torsten]

Please include the function f in which you calculate "N" dependent on "da" callable as

N = f(da)

[shiv gaur]

Dr torsten N is the root ie p function is dependent on da also the change in the value of p for different value of da for some value of na that is given this is running program so Dr torsten pl help the by varying the value of da the value of changing p ie N so we have to calculate the DP/Dda vs da plot

this is running script so minor mistake in program so pl help for three value of na 1.45 1.49 1.50

DP/Dda vs da thanks fig like att .file