Help with plotting 2d and oscillation

5 views (last 30 days)
Jacob Andréasson
Jacob Andréasson on 9 Nov 2019
Skärmavbild 2019-11-09 kl. 17.33.06.png
Hey!
I need som help with the script for exercise C. I dont know if i have understood the exercise right and i would love som help if my scipt is correct.
Here is my script:
%% Exercise 2.a
clear
clf
clc
k=1
L=6
x=linspace(0,L);
fplot(@(x) (((2.*k)./L).*x), [0,(L/2)])
hold on
fplot(@(x) (((2.*k)./L).*(L-x)), [(L/2),L])
%% Exercise 2.b
clear
clf
clc
k=1;
L=6;
T=10;
x=linspace(0,L,100); t=linspace(0,T,100);
n=1:100;
for i=1 : length(x)
for j=1 : length(t)
u(i,j)=(8/(pi^2)).*sum((1./(n.^2)).*sin(n.*pi/2).*cos(n.*pi.*t(j)/6).*sin(n.*pi.*x(i)/6));
end
end
surf(t,x,u)
xlabel('t'); ylabel('x'); zlabel('u')
title('Exercise 2.b');
%% Exercise 2.c
clear
clf
clc
k=1;
L=6;
T=10;
x=linspace(0,L,100); t=linspace(0,T,100);
n=1:100;
for i=1 : length(x)
for j=1 : length(t)
u(i,j)=(8/(pi^2)).*sum((1./(n.^2)).*sin(n.*pi/2).*cos(n.*pi.*t(j)/6).*sin(n.*pi.*x(i)/6));
end
end
plot(x,u)
xlabel('x'); ylabel('u');
title('Exercise 2.c');

Answers (0)

Categories

Find more on Mathematics in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!