clear all; close all;
syms T U t w
a=0.3; M=1; h0=1; h1=1; k0=1; k1=1; p=0.5; d=0.5; Z=10;
%t=0:1:10; %w=0.75;
H=@(s) (p*Z-d-exp(-s))*k1;
K=@(r,g) p*Z*k1-(d+exp(-r))*k1-2*exp(a*r).*g;
for n=1:10
y=int(H(U),0,10);
yy= h0 + ((1-a)./M)*H(t).*h1 + (a./M)*y;
h1=yy;
end
for n=1:10
z=int(K(U,T),0,10);
zz= k0 + ((1-w)./M).*K(t,yy).*k1 + (w./M)*z;
k1=zz;
end
zzfcn = matlabFunction(zz)
t=0:0.1:10;
w = 0.1:0.1:0.75;
[T,W] = ndgrid(t,w);
figure
surf(T, W, zzfcn(T,W))
xlabel(‘t’)
ylabel(‘omega’)
zlabel(‘k(t)’)
set(gca, ‘ZScale’,’log’)
colormap(turbo)
colorbarclear all; close all;
syms T U t w
a=0.3; M=1; h0=1; h1=1; k0=1; k1=1; p=0.5; d=0.5; Z=10;
%t=0:1:10; %w=0.75;
H=@(s) (p*Z-d-exp(-s))*k1;
K=@(r,g) p*Z*k1-(d+exp(-r))*k1-2*exp(a*r).*g;
for n=1:10
y=int(H(U),0,10);
yy= h0 + ((1-a)./M)*H(t).*h1 + (a./M)*y;
h1=yy;
end
for n=1:10
z=int(K(U,T),0,10);
zz= k0 + ((1-w)./M).*K(t,yy).*k1 + (w./M)*z;
k1=zz;
end
zzfcn = matlabFunction(zz)
t=0:0.1:10;
w = 0.1:0.1:0.75;
[T,W] = ndgrid(t,w);
figure
surf(T, W, zzfcn(T,W))
xlabel(‘t’)
ylabel(‘omega’)
zlabel(‘k(t)’)
set(gca, ‘ZScale’,’log’)
colormap(turbo)
colorbar clear all; close all;
syms T U t w
a=0.3; M=1; h0=1; h1=1; k0=1; k1=1; p=0.5; d=0.5; Z=10;
%t=0:1:10; %w=0.75;
H=@(s) (p*Z-d-exp(-s))*k1;
K=@(r,g) p*Z*k1-(d+exp(-r))*k1-2*exp(a*r).*g;
for n=1:10
y=int(H(U),0,10);
yy= h0 + ((1-a)./M)*H(t).*h1 + (a./M)*y;
h1=yy;
end
for n=1:10
z=int(K(U,T),0,10);
zz= k0 + ((1-w)./M).*K(t,yy).*k1 + (w./M)*z;
k1=zz;
end
zzfcn = matlabFunction(zz)
t=0:0.1:10;
w = 0.1:0.1:0.75;
[T,W] = ndgrid(t,w);
figure
surf(T, W, zzfcn(T,W))
xlabel(‘t’)
ylabel(‘omega’)
zlabel(‘k(t)’)
set(gca, ‘ZScale’,’log’)
colormap(turbo)
colorbar function MATLAB Answers — New Questions

​