Answered

Heat transfer differential equation

Call the function this way: [t,y] = ode45(@f,[0,5],1673)

Heat transfer differential equation

Call the function this way: [t,y] = ode45(@f,[0,5],1673)

3 Tage ago | 0

Answered

2D information to 1D colorplot

Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

2D information to 1D colorplot

Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

3 Tage ago | 0

Answered

Quiver with a system of ODEs.

For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

Quiver with a system of ODEs.

For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

4 Tage ago | 0

Answered

Solving nonLinear equation and Linearized Eqaution using ODE45

YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

Solving nonLinear equation and Linearized Eqaution using ODE45

YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

4 Tage ago | 0

| accepted

Answered

How to plot magnitude and directions between 2 points?

try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

How to plot magnitude and directions between 2 points?

try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

4 Tage ago | 0

| accepted

Answered

surface fixed point value Z

I want to know the value of Z from the surfc function Here is it

surface fixed point value Z

I want to know the value of Z from the surfc function Here is it

7 Tage ago | 0

Answered

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?

Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?

Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

11 Tage ago | 1

| accepted

Answered

I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.

Try @(V,T)

I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.

Try @(V,T)

13 Tage ago | 0

| accepted

Answered

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?

alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?

alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

14 Tage ago | 0

Answered

Euler 3D rotation between two vectors

Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

Euler 3D rotation between two vectors

Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

14 Tage ago | 0

| accepted

Answered

Simultaneous differential equations - derivative of y^3 wrt t

Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

Simultaneous differential equations - derivative of y^3 wrt t

Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

14 Tage ago | 0

Answered

Numerically derive a continous, non-symbolic function

Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

Numerically derive a continous, non-symbolic function

Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

16 Tage ago | 0

Answered

3D Point Cloud - Gaussian Curvature

one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

3D Point Cloud - Gaussian Curvature

one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

17 Tage ago | 0

Answered

Runge Kutta for system of eqs

I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

Runge Kutta for system of eqs

I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

18 Tage ago | 0

Answered

The problem of finding intersection point of triangle mesh

To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

The problem of finding intersection point of triangle mesh

To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

19 Tage ago | 0

| accepted

Answered

How can i solve a system which contains matrices

Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

How can i solve a system which contains matrices

Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

21 Tage ago | 1

| accepted

Answered

Using a for-loop to score values in different intervals

Correct form if 0 <= x(i) && x(i) < 1

Using a for-loop to score values in different intervals

Correct form if 0 <= x(i) && x(i) < 1

22 Tage ago | 0

| accepted

Answered

Intersection between line and circle when line ends inside circle (using geom2d)

Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

Intersection between line and circle when line ends inside circle (using geom2d)

Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

22 Tage ago | 1

| accepted

Answered

Removing points from meshgrid (mask)

Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

Removing points from meshgrid (mask)

Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

23 Tage ago | 0

Answered

How can I find the surface normal of a 8 by 3 matrix?

Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

How can I find the surface normal of a 8 by 3 matrix?

Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

23 Tage ago | 0

Answered

How do I increase the circle size and animate a stick figure?

Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

How do I increase the circle size and animate a stick figure?

Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

23 Tage ago | 0

Answered

Finding specific data point x for given y

I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

Finding specific data point x for given y

I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

23 Tage ago | 0

Answered

Surface between 3D lines

You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

Surface between 3D lines

You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

23 Tage ago | 0

| accepted

Answered

two body problem using ode45

- equation of equilibrium Is this condition fullfield?

two body problem using ode45

- equation of equilibrium Is this condition fullfield?

24 Tage ago | 0

Answered

Numerically Integrating the differential equation below?

Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

Numerically Integrating the differential equation below?

Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

25 Tage ago | 0

Answered

Streamline plot not complete

Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

Streamline plot not complete

Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

25 Tage ago | 0

| accepted

Answered

2D Matrix revolution

Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

2D Matrix revolution

Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

25 Tage ago | 1

Answered

How do I align the vertices and the adjacency graph of a Buckyball?

Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

How do I align the vertices and the adjacency graph of a Buckyball?

Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

27 Tage ago | 0

Answered

Find the coordinates of a point chosen on a plot

Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

Find the coordinates of a point chosen on a plot

Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

29 Tage ago | 1

| accepted

Answered

How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?

Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?

Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

etwa ein Monat ago | 0